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Preface 

This  thesis  demonstrated  the  feasibility  of  analyzing  the  output  intensity  pat¬ 
terns  of  an  optical  fiber  to  determine  the  angle  of  arrival(AOA)  of  incichuit  l;!s<'r 
radiation.  Adequate  information  was  contained  in  the  low  fretpiency  Fourier  com¬ 
ponents  to  determine  .‘\OA  within  1°. 
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Ahhitract 

riiis  tlicsis  f’xainincd  thr  foasioilit y  of  using  oi)tical  filx'rs  for  unglc  of  .a- 
ri\al(AOA)  icju.  S])(H'ifically.  a  unifonii  ainplit  ucU'  plain'  wave  was  t  ransniil  1  cil 

tlu’ough  a  2tniu  diaiiK'tc'r,  tnull iinod*',  stt'p-iiidox,  p'aslir  filx'r.  Eihcr  l('ngtlis  id  Mn. 
2/11.  and  ‘,Mi/  wore  invostigatc<l.  I  lu' out  put  inti'iisiU’  pattern  Iroin  each  optical  Idier 
was  detected  l)y  a  512  x  512  rest)lut.ion,  cliarge  injt'ction  d('\’ic('(('ll))  camera,  i  lu' 
data  array  from  the  canK'ia  was  [)reprocessed  with  a  nu'dian  filter,  and  normalizi'd 
to  the  peak  int('nsity.  .\  fast  Eouri<'r  traiisform  coinx'itf'd  tlu'  rlata  to  the  frequency 
domain  wlu're  tin'  fundann'iital  and  first  thr('<'  harmonic  fr('qu('nci<'s  w('r('  (-xtracti'd. 
'1  h('se  1!)  numhers  represented  comj)onents  of  a  featur<'  vector  re|)r('s<'nting  tlu'  |)al- 
terii  related  to  a  spc'cific  AOA.  'I’wenty-  six  t.emj)lat<'  vectors,  (-ach  representing  a 
1°  incrc'ment  from  0°  to  25°.  were  created,  'lest  vt'ctors  wc'ix'  conquiicd  to  template 
Vi'ciors  using  ihi<'lid<'an  distance.  I  lie  /ninimnm  distance  classified  the  lest  \'ecloi' 
.\()A. 
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ANGLE  OF  ARRIVAL  DETECTION 
TFIROUGII  ANALYSIS  OF 
OPTICAL  FIBER  INTENSITY  PATTERNS 


/.  I II  trod  uc  lion 

I  h<‘  Air  I''c)rc'<'  is  currcut  ly  in\<‘st  it^nl  iii.i;  iiK'tluxIs  lo  (Ictcriiinic  tin'  aiiulc  (j1 
arrival!  AO  A )  of  diroct  od  cncrtry  wc'apoiis,  primarils'  lascis.  An  AOA  di'icci  or  lias 
s('\'('ral  applications.  It  conid  Ix'  ns('d  on  an  aircraft  to  alert  tlic  pilot  to  incident 
laser  radiation  whicli  conid  blind  crc'wiiMnnbers  jx'rniain'nt  ly  or  for  as  little  as  a  t('w 
seconds,  tlu'rohy  initting  their  lives  in  j<H)|)ardv.  Atiollier  application  is  loi'  use  on 
I'econnaissancr'  satellites  to  i\'arn  of  )as('r  radiation  which  conid  daniaxe  \  ital  coinpo- 
nents.  .An  AOA  detector  can  Ix' ns('d  in  conjunct  ion  with  pi)W('r  met  ers  ainl  sj.ectrntn 
analvzc'rs  to  determine  pow('r  ami  fre<m<'nc\'  of  th<'  sonrc(\  which  conid  be  se\'i'ral 
hundred  kilom('t<'rs  away.  I'his  information  is  vital  to  tin'  intelligence  eoinmnnity. 
riiis  system  must  be  com|)act,  li, gill  weight  sliock  rc'sistant.  highly  reli.dile.  and  sim¬ 
ple  in  constrnct ion,  with  an  angtdar  resohilion  of  1  dc'gia'c  or  less  (tt:7).  I  his  thesis 
will  in\'estigate  a  solution  to  this  [H'oblem  (hrongh  the  use  ol  a  passix'e  (aei<'d  np.m) 
detector,  riie  output  intensit}'  pattern  from  a  large  diametei'  optic. d  liber,  specili- 
cally  a  ■]injn  diameter,  mnltimode.  plastic  lilxM'  is  analyzed  to  determine  the  .\(),\. 
Since  this  detc’ctor  will  nse  short  libersfi.e.  Sent  or  less)  and  will  contain  no  moving 
parts,  it  will  nu'et  the  previonslv  iiK'ntioned  nxjnirement s. 

/./  Approach  and  Scopr 

riie  initial  efhrri  in  this  thesis  is  to  research  the  leasibility  ol  .\().\  systems 
using  o|)tical  fibers.  I'his  includes  a  t  In'oret  i<'al  determination  ol  the  intensity  pioiile 
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at  t  1k'  out  put  of  I  he  u|)t  ical  lilx'r  using  g<s)in('t  rical  opt  i<s.  1  lie  sc( oud  pa  1 1 
ol  this  rc'sc'arcli  is  ('xporiiiiciit al.  A  uniform  amplitude  plane  \\a\'<’  will  dlnmiiiati- 
a  short  s('gm(  nt  of  ■\iinn  diameter  plastic  (ilx-r  with  the  output  impinging  on  a 
old  X  old  cl!arg('  inj('ction  (l(‘viee(( '1 1))  canu'ra.  Since  this  detector  :  to  determine 
AOA  of  laser  input  at  a  suhstantial  distance  from  the  sourci',  the  in|)ut  is  as>umed 
to  l)('  a  uniform  amj)litude  jilane  wave.  I'lie  data  is  then  rixluced  to  a  'dot)  x  deli 
array,  smootlu'd  through  a  medi<m  hltc'r  ojX'ration.  and  then  normalized  based  (jn 
the  highest  intensity  \alue.  following  this,  tin'  ix'sults  are  f'ourier  t  ranslormed.  1  his 
l!e<|uency  domain  jirocessing  iiwoUcs  ext  ract  ion  of  the  lundaini’iit  al  fii'(|uenc\'  and 
the  Inst  three  harmonics  to  torin  a  lea  t  lire  vi'ctor  for  each  aiiule  ol  nicidencc,  f  nia  lly 
ih'  .'C  ji'alure  \'eciojs  are  compared,  using  .Minkowski  di>tance.  to  randi)in  angle  test 
\'e(  tors  ill  order  to  determiiK'  .\().\  of  the  test  samples.  1  his  proci'ss  is  (/iitliiied  in 
figure  1 . 

/.d  ( )r(/(uii:(i/ion 

1  his  thesis  hegnis  with  a  look  at  two  ,\().\  detectors  utilizing  optical  lihers. 
It  then  t  heoix't  i(  ally  in\'es|  igat  es  tin'  iirojiagat  ion  of  light  through  .iii  optical  liher 
to  explain  tin  caialioid  output  pattern  si'iui  in  short  lihers.  Chapter  d  ilescrihes  the 
lahoratory  setu|).  data  acquisition,  ami  iirocessing.  Results  ate  discussed  in  ('ha|)t('r 
1  with  conclusions  and  re<-ontnieti<lat  ions  in  chapter  o. 


II.  Swrnniary  of  Current  Kriowledyt 


\arious  systems  ha\('  tx'c'ii  dc'Vi'lopcd  in  r('reiit  yc'ars  to  determine  AOA  of 
liyht  energy.  One  such  system  uses  a  (|uadraii<  detector  to  gi\'e  enmigli  information 
to  allow  armonsi  vehicle  crews  to  home  in  on  tlie  .soni’ce  using  the  vi'l  icle's  night 
si^ht  (■'):1  ).  Another  system  utilizes  interferometry  tt'chniques  to  determiiK'  ikjI  only 
.\().\  of  the  incident  laser  radiation  l)ut  also  its  wavelength  (5:1).  .A  third  d('\ic(' 
imes  th('  g('ometry  of  a  ti'iangular  culx'  rellector  with  jdiotodc-tectors  aligned  with 
the  thr('esid('s  to  determine  AOA  (o:l-IO).  I  here  aif'only  two  .AO.A  systi-ms  using 
o|)t  ical  lilx'i's  (t).  Id). 

J.  I  Optical  |■  lb(  V  Di  !<  ctors 

I  he  first  system  was  inv('st  igat('d  in  198-5  at  Wright  -  Pat  t  (uson  Ah'H.  01!.  This 
system  utilized  an  array  of  l<'ns('s  moun(c>d  on  a  hemis[)herical  dome.  Tach  lens  was 
separated  l>y  'J‘J.5  degia'es  allowing  ovc'rlapping  fi('lds  of  vi('w(  !•  0\’).  In  t)t  her  words, 
eaeh  h’lis  saw  a  small  poi  t  ion  of  what  its  adjaccuit  lens  saw.  1  his  resnlteil  in  a  total 
t)l  ^1  lens  apertures.  ITich  lens  focuscxl  light  etiergy  into  a  colh'ction  of  optical  filx'is 
ol  diflei'ent  lengths.  Pacli  Inmdle  of  libers  t('rminated  onto  a  photodetector.  Two 
dilleicnt  t.N'Jx's  (d  photodeteclors  were  used.  .Silicon  [)hot odet eel ois  writ'  used  to 
detect  laser  wa\('lengths  from  l).()////i  to  l.()//n)  whiU*  germanium  photodelectors  op¬ 
erate  .1  m  the  1.1/////  to  1.9/;///  s|x’ct  ral  region,  commonl}'  known  as  the  infrared(  I H ). 

I  he  phot  odet  ect  oi  sigintls  were  used  to  e.xcite  mnh'rdamix'd  circuits  (M);',!).  Due 
to  the  overlap  ol  the  i'()\  ol  adjac/'iil  lensv's.  a  mathematical  rehal  ionship  I/elw/'en 
tie  /e/'o  ci(/ssings  ol  liinh'i'dainped  circint  oscillations  and  the  optic;d  cenli(/id  of 
/ipertnre  illnmiiiiit  ion  was  ex|)loited.  .Signal-to-nois('  was  alt/ned  ;is  necessary  to  v;u  y 
(h'lei  tion  |/rob;diility  ;tinl  fidse  alarm  rate' ( H);()).  desis  were  cenidncted  with  v.irying 
<h"ji'<'<'s  anel  type's  of  baekgrounel  lighting.  IT'sulls  showe'd  the’  baikgre/unel  lighting 
liael  negligible' e'lle'cl  een  .\ ( ) .A  e let e 'ct  ie/ii  of  t he' hise'i'  1  igh t  (  1  ():55-.'55 ) ,  ,\ngle'e)l  arriv;d 
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accuracy  of  ±5  (i('grccs  was  obtaiucd  'I’liis  is  good.  bu(  it  is  desired  (o  )/a\e 

an  accuracy  of  1  degree  or  less 

riic  second  system  is  the  s\d)jecl  of  current  researcli  by  a  small  linn  based  in 
.\l('ido  I’ark,  California  and  hope's  to  obtain  the  needee!  a.'cnracy  (fill'd),  d  ins  syste'in 
is  composed  of  several  hundred  apertures  on  the  surface  of  a  ln'mispherical  dome. 
Ifach  ape'rture  has  a  20  degre'e'  I’OV  and  does  not  reepiire  lense's  like  the'  pr('\ious 
system  (0:5).  This  system  uses  optical  fibers  to  transmit  the  intc'iisity  incident  on 
th('  ape'rt tires  to  individual  silicon  photodetectors.  A  processor  using  an  algorilhm 
similar  to  that  used  in  Radar  Warning  Receive'rs  performs  the  necessary  comptirisons 
of  iutc'usities  from  th<'  irradiatt'd  filx'rs  to  dt'lermiiH'  .\0.\  (f):7).  It  accom|)lishes  this 
by  utilizing  a  look-up  table,  d  his  t.abh'  was  creatt'd  by  illuminating  the  hemis])h<'re 
frtun  different  dirt'ctions,  pt'rhirming  the  aforemi'iit iont'd  signal  anal\'sis.  and  storing 
the  rt'sults  in  memory.  In  act  ual  us('.  th<'  system  will  corrt'latt' t  In'  int ensity  of  incident 
last'!'  energy  to  the  test  patterns  ston'd  in  memory.  A(),\  of  tin'  last'r  radiation  can 
by  determiiK'd  from  tliis  corrt'lation  (9:5:5-15).  In  order  to  achit'vt'  an  ovt'rall  l''(/\’ 
of  2-  st('radians(hemisphere),  these  photod<'t<'ctors  art'  spact'd  ('vt'iy  !>  ih'gn'es.  once 
again  establishing  o\'('rlap|)ing  bOV  (f):8).  'rids  ovt'i  lap  is  ni'ct'ssary  so  t  hat  tin'  ent  irt' 
2r  sif'radian  fit'ld  is  deteett'd.  It  also  insures  that  each  st'iisor  is  not  iinh'pi'inh'iit 
but  interconnected  so  that  the  data  rt'ct'ivt'd  can  be  combined  <md  corrt'lated  m  an 
ellecti\'e  mamn'r  (9:18-25).  1  he  ojttical  (iln'rs  in  tln'st'  s\stems  an'  strictK'  usi'd  to 

transmit  tin'  iiitt'iisity  from  ;m  apt'iture  with  a  snudl  l'()\  . 

I.ifjhi  Propdijalinii  in  an  Opticdl  Pihiv 

1  his  research  e.xphdts  the  amilysis  of  the  intensity  proliles  uf  individual  non 
imaging  filx'is  to  gain  additional  iniornuit ion  at  the  output  t)l  the  liber,  ('aptaiii 
(leiio  Coh'  was  the  first  |)erson  to  invt'stigate  this  t  heoi  \'.  lids  thesis  will  perlorm 
the  experimental  portion  of  (’a|)tain  (’oh''s  initial  invest  igiit  ion  (1).  .\ou  imaging 
im|>lies  that  there'  is  no  lens  syste'm  used  to  focus  ;iu  image  jdaue  on  the  deti'ctor 


array  in  tlio  (’ID  camera.  In  afldition,  multimode  step-iiule.x  fibers  do  not  transmit 
image's  well,  due  to  inle'rmodal  dispersion.  A  look  at  how  light  pro])agat('s  down 
a  fil)er  will  help  to  better  understand  this  phenomenon.  Since  the  diameter  of  the' 
lib('r  used  was  significantly  larger  than  the  wavelength  of  the  incident  laser  radiation, 
geometrical  o[)tics  (i.('.  ray  tracing)  can  be  n.sed  to  examine  wave'  propagation  within 
the'  fibe'r  (7:b5).  I’or  the-  analysis,  incielent  light  rays  will  be'  se])arate'd  into  two 
classes.  The  first  class,  meridieinal  rays,  are  easy  to  track  since'  t  he'y  stay  wit  hin 
the'  meridional  plane'.  A  meridional  plane  is  any  plane  which  contains  tlu'  axis  of 
s\nimetry  of  the'  liber  (core'  axis).  Ideally,  the'ie'  is  only  eine'  me'iielieinal  ray  fear  an 
ine  ieh'iit  plane'  waN'e.  'flu'  se'ee.nel  elass,  skew  rays,  re'pre'se'iit s  all  eiff-axis  rays  (7;2.'i). 
f  igure'  2  displays  the'  merielie)nal  ray  passing  threnigh  a  fibe'r.  ,Sne'irs  law 

II i  sin  0,  —  II (  sin  0, 


is  n.se'el  to  ele'le'rmine'  its  [jatli  as  it  pre)j)agale's  eleiwn  the'  fibe'r  (7:17).  As  the-  ray 
ceintae  ts  the  e'litrance  of  the'  (ibe'r,  seime'  of  the'  e'lie'rgy  is  re'lle'ct e'el  anel  the'  re'maineler 
is  transmitted.  As  the'  angle’  eef  inciele’iice  incre'ase's,  h'ss  e'lu'rgy  is  t  ransmit  t  e'el  into 
the’  filje-r  (I;  1  Oti).  The’  transmit  t  e'el  ray  is  bent  towarels  the'  core'  axis  si  ne'e’ t  he  ineh'X  eif 
re’fraction  erf  the  ceire’  is  gre'ate’r  t  han  that  erf  air.  I'his  is  Du’  case'  of  e'xternal  re’lie'e  t  iern. 
In  the  case  erf  e’xte’rnal  refle’ct  iern,  some  erf  the-  e'liergy  is  always  t  reinsmit  te’el.  As  the’ 
ray  contacts  the  core-claeleling  inte'iface  it  is  once  again  re’lle'e'te’el.  heiwe’ve'r,  this  time 
the'  ray  is  cerntacting  the  snrfae-e'  erf  a  h'ss  erptie-ally  ele'iise  me'elinm  (i.e’.  hrwe'r  ineh'X 
of  re’fraet  iern)  which  repre'se’iits  the'  case  erf  inte’rnal  re'fh'cl  iern.  dims  the'  ray  is  be’iit 
away  frerm  the  snrfaee’  nerrmal.  The'  |rerint  erf  cerntaet  e  an  be'  merele’le'd  as  a  plane' 
tange'iit  ter  the’  |rerint  erf  e'erntae'l  whieh  ean  Ire'  ealle'il  the'  re'lle'et  ion  plane'.  .As  the’ 
angle’  d  de’e-re’ase's.  er  ine're'ase’s  imlil  it  e’e(uals  bt)  ele’gre'e's  and  the'  I  ransmit  le'd  ray 
is  partdh'l  ter  the’  eerre’-e'laehling  inte’ilae'e’.  This  \alue'  erf  o  is  knerwn  as  the’  e  ritie'al 
angle’.  0,rii-  for  any  value’  erf  o  h'ss  than  Omt-  the'  ray  will  e’xire'rie'iie'e'  tertal  inte'rind 


re'lle’et  iern.  lele’ally  this  me’ans  that  neriie'  erf  the'  ray  is  t  I'ansmit  t  e'd  inter  the’  elaihling. 


Figurr  2.  Path  of  Mrri<liotial  Hay 


but  in  reality  some  minimal  loss  always  occurs  in  the  form  of  a  surface  wave  or 
evuiiescent  ,vave  (4:100,107).  The  ray  continues  down  the  fiber  experiencing  intcunal 
reflection  at  each  core-cladding  interface  until  it  reaches  the  end  f'f  the  fiber.  I'lie 
ray  contacts  the  end  of  the  fiber  at  an  angle  of  (j)  degrees  and  once  again  experiences 
internal  reflection  since  n.a,>  <  ncore-  Thus  the  ray  upon  exiting  the  fiber  will  be  bent 
away  from  the  core  axis.  Figure  3  presents  a  graphical  representation  of  this  angular 
dependence  (7:17-24). 

The  maximum  entrance  angle  which  allows  total  internal  reflection  can  be 
calculated  by  Snell’s  law.  Using  figure  3  (b)  as  a  reference: 

iicoTf  sin  (90°  -  (pcrxi )  =  y>chd  a  ( 1 ) 

where  a  =  90°.  Thus: 

sin  (90°  -  dcrit)  =  =  fos  (2) 

^^corc 

for  the  core-cladding  interface. 

At  the  fiber  entrance: 


Sin  OjixQj-  —  ^Irorf  1^111 


and  since: 


eriuation  3  becomes: 


sin 


A,  ril  —  ^1  COS*  Ocrit^ 


(;5) 


(1) 


Our  Sill  —  ^^I'orr  ^1  (OS 
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refleciion 

piano 


(a)  End  View 


skew  ray  path 


ray  path  projected 
on  fiber  end  face 


skew  ray  path 


core  axis 


(b)  3-D  View 


Figure  4.  Path  of  Skew  Ray 


Substit  uting  equation  2  into  equation  5  yields: 


^^air  Sin  0 —  ^^coTf^  1 


Since  n„,>  =  1,  eciuation  6  reduces  to: 


=  siif 


l  or  the  fiber  used  in  this  research.  =  1.492  and  =  1 .  1 1  (i  n'sulting  in  a 

maximum  AOA  of  approximately  28°  to  achieve  the  critical  angh'  at  t  he  core-cladding 
inlerface(7:24.2.'5).  Figure  4  illus(rat<'s  the  |>ath  of  a  ski'w  ray.  Fhis  ray  ('xpeiieiices 


meridional  rays 


skew  ravs 


kigure  5.  Critical  Angle  ('oni|)arisoii  of  Meridional  and  Skew  Kays 

the  sanu’  reflection  and  refraction  as  tlie  meridional  ray  as  it.  contacts  the  I'lit  ranee  to 
tlu'  fiber  since  it.  is  part  of  the  same  uniform  amplitude  jilane  wavi-.  It  continues  to 
tra\-el  parallel  to  the  meridional  ray  until  it  reaches  the  first  core-cladding  interfaci'. 
At  this  interface  it  reflects  diffi'i'ently  tlian  the  nuuddional  ray  since  its  ndlection 
plane  is  oriented  differently,  as  slunvn  in  figure  4  (a).  I’iguri'  4  (b)  shows  how  this 
ra\'  continues  its  projiagation  along  the  fiber. 

Figure  b  is  a  magnifii'd  end  viiwv  of  a  small  sample  of  adjaccuit  rays  composing 
the  incident  uniform  amplitude  plane  wave.  I'he  c('nter  ray  is  tlu'  meridional  ray. 


('ardioid  Output  Pattrni 


Figure  6.  Cardioid  Output  Intensity  Pattern  of  an  0])ti(al  I'iher 

All  other  rays  arc  skew  rays.  These  views  only  show  the  rays  as  tlnyv  enti'r  the  filu'r 
and  contact  the  core-cladding  interface  for  the  first  tinu'.  (’ontinuing  propagation 
down  the  fiber  is  not  shown.  Figure  •'>  (.a)  shows  that  if  the  incoming  nu'riditinal 
ray  is  at  the  critical  angle,  the  skew  rays  are  less  than  the  critical  angle  since'  6  .  in 
figure  3,  is  now  divided  into  a  tiHiisversc  as  well  as  a  longitudinal  coiupoiK'iit .  Thus 
the  only  ray  at  the  critical  angle  6cTit  is  tlie  meridional  ray.  If  the  angle  of  incidf-nce 
of  the  plane  wave  is  increased  slightly,  the  meridional  ray  as  we'll  as  the'  imme'eliat ely 
aeljacent  skew  rays  are  >  (f>cri/-  1  his  is  shown  in  figure'  5  (b).  f  igure'  •")  (c)  she)ws 
that  as  the  angle  of  incidence  e-ontinue's  to  incre'ase.  jnergression  eeent  inue's  outward 
Ireem  the  eenter  with  more  and  me>re’  skew  rays  reflecting  at  angle's  >  I'hgure'  G 

elisplays  the  enitput  ed  a  large’  eliamete'r,  multime)ele’,  ste'|)  inde'X  fibe'r.  1  he'  angle’ 
e)f  ineiele'iice’  is  much  less  than  the’  critical  angle’  anel  the’  fibe'r  is  e»nl\'  le)ng  e'Uongh 
te)  allow  one’  re’fh’e’t ie)n  at  the’  ee)re’-  elaeleling  inte'ifaee',  ,\s  e'N|)e'e  te-d.  the'  me'riiliiinal 
ray  passes  threjugh  the’  e-nel  ejf  the’  fibe’r  at  the’  highe'st  le)ngit  nelin.t!  angle'  while  the’ 
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skew  rays  sliare  their  longitudinal  angle  with  a  transverse  angle,  thus  reducing  the 
longitudinal  angle.  This  results  in  the  cardioid  pattern  seen  in  figure  6.  This  cardioid 
pattern  was  first  predicted  by  a  ray  tracing  program,  sensor. pas.  written  in  rurbo 
Pascal  by  J.  Scott  Hunt,  Cadet,  United  States  Air  Force  Academy.  Cadet  Hunt 
wrote  this  program  earlier  this  year  during  a  temporary  duty  assignment  at  the 
Air  Force  Weajions  I,ab.  As  the  angle  of  incidence  increases,  eventually  a  second 
refl('ction  is  allowed  to  take  place  in  the  fiber.  During  this  time  the  cardioid  pattern 
moves  through  the  center  of  the  output  to  the  other  side.  As  the  angle  of  incidence 
increases  beyond  the  critical  angle,  the  cardioid  |)attern  continues  to  move  as  before 
but  e.xperiences  a  loss  of  intensity  due  to  le.ss  transmission  (higher  refl('ction)  at  the 
entrance  of  the  fiber  and  loss  of  total  internal  reflection  within  the  fiber. 

2. 3  Prcproccssin  g 

do  remove  noise  spikes  from  images  and  smooth  th('  data,  a  iiK'ilian  filtc'iing 
opc'iation  can  be  performed.  I  he  process  is  simple.  For  an  .\’  x  nu'dian  fill<'r  the 
intensities  in  the  first  ,'V  rows  and  A’  columns  of  the  image  array  are  read.  I'liese 
.V*  numlters  are  then  sorted  lowest  to  highest.  The  middle  number  r('])rescnts  the 
median  value  and  it  is  placed  in  the  center  of  the  N  x  A'  array.  .A  new  image  array 
is  created  and  begins  with  the  .V  x  N  array.  Going  back  to  the  original  array,  the 
next  set  of  A’^  numbers  is  still  taken  from  the  first  N  rows  but  shifted  one  column 
to  the  right,  thereby  using  intensities  from  columns  one  to  .V  +  1.  The  process 
continues,  storing  the  median  values  in  the  new  array.  This  type  of  filter  will  not 
affect  step  changes  in  intensity  but  will  adjust  a  noisy  pixel  to  a  value  rt'presuntat  ive 
of  its  neighbors.  An  (A'  +  M)  x  (A'  +  .1/)  array  will  smooth  more  than  an  \  x  .V 
array,  howev('r.  it  may  take  significantly  longer  to  proc('ss  (fi). 


2.-j  Data  Reduction 

In  order  to  increase  processing  s|)eed,  sonu'  form  of  data  recliiction  must  Ke 
employed.  Information  about  tin'  shape  and  magnitude  of  tlu'  output  of  an  optical 
fiber  contained  in  a  large  array  of  pixels(i.e.  250  x  250)  can  be  reduced  by  I'onrier 
dVansform(F"r)  methods.  A  F'F  decomposes  a  waveshape'  into  an  inlinite-ly  countable' 
set  of  sines  and  ce)sine's.  The  f'd’  e'einatie)n  is  (2:31  il): 

/•’(u,c)  =  ^{/(.r,//)}  =  r  r  /(.r,,v)r-J''^‘^'''+'^')e/.re///  (S) 

J J  -  'X' 

I  he  ge'iie'ral  shape'  e)f  the'  i>atte'rn  is  ce)ntaine'ei  in  the  lowe-r  lree|ue'ncie's  e)l  the-  f  1 
while'  more'  eh'taileel  infe)rmat  ie)n,  such  as  sharp  e'e)rne'rs.  is  containe-ei  in  the-  liighe-r 
fie'(|ue'ne'y  ce)m()one'nts  of  the'  !•' I'  (3:202.203).  As  an  e'xample',  the:  sepiare'  wa\'e  in 
(igure  7  (a)  cati  be*  re'|)rese'nte'el  by  a  sum  of  sine's  and  cosines.  In  this  s])e'e-ific  e’xtmiple' 
the  coe'flicie'iits  of  the  cosine  te'rms  as  well  as  the  e've'ii  harmeniic  sine'  te'rms  are'  ze-ro. 
Thus  only  oeid  harmonic  sine'  wave's  are  ne'cele'el  te>  re’f)re'se’n(  tiiis  eeaee’forttt.  In  the' 
ge’iie'r  d  e  ase'.  he)weve'r.  all  ceee'fllcie'nt s  will  have  ne)n  ze'ie)  value-s.  As  sheewn  in  figure'  7 
(b).  the'  gene'ral  shape'  e)f  the'  sepiare'  wave'  can  be'  ae'hie've'e!  using  e)nly  the'  lirst  li'W 
lowe'r  freepiency  conijieene'nt  s.  In  the'  e'ase'  e)f  a  ste)re’el  eligital  ;uray  e)(  inte'iisit  ie's. 
the'  array  ean  be  converle'd  te)  the'  fre'e|ne'ncv  ele)maiii  by  me'ans  eel  a  last  hourie'r 
t  ransfe)rm(I'T’T).  Altheeugh  the' input  array  may  be'  re'al  niimbe'rs  e)nly.  the'euitpiit  e)f 
the'  I  f'  1  will  be'  a  re'al  array  anel  an  imaginary  arra\'.  l  lie'se'  arrays  are'  re'ielaee'd  with 
a  single'  magnitnele'  array  by  e'eembining  e'ae'h  array  e'le'iue'ut  thre)Ugh  the'  hermula: 

M aijuit ad<  ~  \J{  R(  +  [  fnt )-’ 

1  his  magnitnele'  array  ee)ntains  the'  le)we'r  fre'epie'ue  v  ce)mpe)ne'nl s  in  e'aeh  e)f  the-  feeur 
eeerne'is.  In  oreh'r  te)  place  the'  ee’iite'r  fre'e(ue'ne  v  at  the'  ee'iite'r  ol  the'  array  surreeuiieh'il 
by  the'  le)we'r  harmeenies  with  fre'epie'iicy  ine-reasing  euit warel  freem  the'  e  e'iite'i  .  the'  lirst 
anel  thirel  as  we'll  as  the'  se'e'eenel  anel  fe)nrth  epiaelrants  must  be'  swa|)pe'ei,  eeunmouly 
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railed  ciuadrant  swapping,  din'  fuiidaim-iital  freciuenry  is  ncav  in  the  rentei'.  1  he 
second  harmonic  is  composed  of  the  eight  mimhers  snironnding  tin'  cenler.  I  he 
third  liarmonic  consists  of  the  sixt('en  mimhers  surrounding  the  second  harmonic. 
I'his  process  contitiues  outward  to  the  edge  of  tin'  array.  1  lu'ri'fore.  iidormalion 
regarding  tlie  general  shape  of  a  large  array  can  he  stored  in  tin'  fiist  h'w  harmonics 
of  a  tnagnitnde  array  of  spatial  frecpiencies. 

d.  j  Compari.^on  of  lma(j(s 

Since  informat  ion  concerning  t  he  shapi'  of  an  image  can  he  stort'd  in  thelowi'r 
harmotiics  of  spatial  freipiencies  these  components  can  repri'sent  the  components  of 
a  s'ectcjf  cotnmonly  calh'd  a  fetitnri'  ci'ctor.  1  his  results  in  a  \'ector  of  .\  compoiu'iits 
for  ('ach  image.  C’omptiritig  the  images  is  Just  a  matter  ol  calculating  distances 
het  wei'ii  .\ -ditnetisional  letiture  vectors.  I  his  is  done  using  tin'  .Minkowski  distanci': 


where  A  is  an  iiilegei'.  I'oi'  A  :  1.  \t \  is  calh'd  hi.n  distance  and  h  i'  A  —  2  it  is 

known  as  i.iicliilf  an  distance.  .\  2  1)  e\am]>l<'  is  shown  in  ligure  S, 

I  he  tips  of  the  two  vectors  lie  at  (2.2)  atid  ('>.1).  laxi  ilistance  ec|U.ds: 

-  (.2  -  2)  +  (  1  -  2)  =  1 

and  the  I'.uchdean  distance  cipials; 

M,  ^  y/j.'.  -  21-  f  (  1  -  2  I-  v'Vtl 

11  the  distance  in  .\  Space  hel  ween  M  vectors  is  large  enough  to  remain  sepai.ited 
when  noise  is  added  to  the  svstem.  1  iie  Minkowski  distance  can  he  used  lo|-  patli-rn 
leeomnt  ion.  1  his  |s  accomplished  hv  assigning  a  template  \ecl(jr  to  ea<  li  know  n 


I'igur('  8.  Taxi  Distaiux'  and  I'urlidoaii  Distance’ 

slia|)('  and  calculat  ing  the  distance  to  a  test  vx’ctor  r('i)r(>s(’nt  ing  an  unknown  slia|)(’. 
1  he  closest  K'niplate’  ve'Ctoi'  classili<‘s  t  h<'  t('st  ve'ctor. 

J.(i  SuiimKtrji 

1  his  cha|)t(’r  examined  two  o[)tical  fiber  syst<'ms  which  use'  the’  jtowe'r  rie)m 
'-e  vcral  libers  to  eletermine  AOA.  Pro|)agation  of  light  radiatieni  through  a  libe-r 
was  |)re’sente’ei  in  order  to  nneh'rstanel  the  cardioid  out|)nt  fiom  a  she)rt  fibe-r.  l  lu'n 
siiiDot  hing.  data  re'elne't  iejii.  anel  a  pat te’rn  recogidt  ion  t e'chniepie'  using  fe-at  lire’  \  e’e  t ens 
wi'ie’  re’viewe’el  in  an  e'fre)rt  te)  eleme)nst rat e’  l<'asibilit \’  ot  ele'te'rmining  AOA  thiejiigh 
an.dysis  of  t  he  inte'nsity  pat  tern  lre)m  '  he  e)ut  ()nt  of  an  opt  ievd  fibe-r.  Chapter  d  elel  ails 
the  labeaatory  setup.  e'e)ui pme-nt  used,  anel  proce'ssing  seaftwai'e'  ns-'e!. 


III.  Expcrrmc.nUil  ApproacJi 


1  lu'  ('xpc'rinK'ut al  pcirtion  of  lliis  tliosis  list'd  a  Siiiricon  fraiiu'  <>,ial)l)('i'.  o|)1i- 
cal  Ix'iich,  atid  softwarr  written  in  '  I’lirlio  ("  and  'ADA'  to  pt'rfonn  tlie  necessary 
piocessing  of  tlu'  data.  All  '  I’nrlio  C  programs  were  written  by  the  author.  1  he 
■.\1).\'  program  list'd  for  the  fast  Fourier  t ransformf FF  1)  was  snp|)iif'd  liy  (It'orge 
l.etonrneau.  Cilv8!ll)  (8).  'I'his  chapter  examines  the  design  of  tin-  tijitical  bench, 
data  collt'ction.  and  t  lit'  [irtrct'ssing  of  that  tlala.  I'lit'  reader  is  assiimetl  to  be  fa¬ 
miliar  with  use  of  an  IR.M  or  compat  iblt' perstrnal  ct)mi)Ut<'r(  l’(')  using  the  MS-DOS 
operat  ing  systt'in. 

/  Spn-icoD  l-'nniK  (Ifahbfr 

1  his  ('xpt'riment  rt'tinirt'd  the  nt'ct'ssary  lianlwart'  and  software  to  analyze  the 
output  intt'iisity  pattern  of  an  optical  fiber  as  it  was  rotatt'd  through  seeeral  angles 
ol  incidi'iict' of  incoming  laser  ratliation.  'I'hc  Sjiiricon  22-”)()  Last'r  Heani  Diagnostics 
Sy'iein  was  chost'ii  since  it  iirovitletl  high  image  resolution.  This  systt'm  was  used 
to  grab  fraitu's  of  tlata,  store  that  data  in  a  file,  and  display  a  3-D  jilot  of  thi'  iniagt.'. 
I  lie  system  consists  t)f  a  (uocessing  carti  which  ii.st's  two  slots  in  an  IB.M  conijuitible 
personal  com|)nter(  I’O).  x't'rsion  u.t)  leaser  Bt'am  Diagnostic  software,  a  ( '1 1 )  b  Ft '1 1 
.'12  .X  •')  1 2  Charge’  Injection  De\ice(( 'll))  came'ra,  a  real-t  ime  display  monitor,  and 
the  necessary  connect  iiig  cables.  .As  just  meiit  ioiic'd.  t  Ih'  22-')I)  syst  em  is  ad\'ert  ised  to 
Work  within  an  IB.M  .A  !  .  IB.M  or  IBM  compatible.  .Alter  much  expe-riment  at  ioii 
with  \arious  BC's.  it  was  d’s  -'ve’ed  that  the22'>()  apjiarently  only  works  on  an  IB.M 
I  .  Ibis  may  be  due  to  an  addi'essing  pidbh'iii  within  the  /,-21.^'s  and,  IB.M  'isti 
ill  which  the  22.''()  system  failed  to  operate'.  De-tails  re'gareliug  the-  ope-ratiou  of  the- 
.8|i)ricon  22.u()  base-r  Be-am  Diagiiostie-  Syste-m  eaii  be-  found  in  appeinlix  .A. 
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Figure  9.  ()|)li(-al  Hencli  Layout 


J  Design  of  Oj)tiral  lUnrh 

Since  this  was  a  new  area  of  rc'sc'arcli,  the  optical  Ix'nch  had  to  he  designed 
and  l)uilt.  d’his  r<’(iuired  the  fabrication  of  several  lixlurc's.  ’Mu'  end  rc'sult  was  a 
sini|)le.  reliable,  and  easy  to  us('  system  for  gat  Inn  ing  data.  Figiiin'  t)  details  tlu' 
ojjt  ical  bench  design. 

Since  this  detector  will  Ix'  used  at  large  distanen's  from  tlx'  souicc'  of  lasc-r 
radiation,  the  input  to  the  detector  was  assumed  to  Ix’  a  imiform  amplitude  plain' 
wav<'.  A  •')0m\V  .\ewport  hniuiprin’iit  C'orporation(.\F(’)  II('-.\('  continuous  wave 
laser  [M'ovich'd  input  to  the  system.  !  In'  Ix'am  passed  through  a  \ariable  attenuator 
and  onto  a  beam  steering  instrument.  This  reduced  Ix-am  intensity  and  made-  the 
job  ol  beam  alignment  signilicant  ly  (-asier.  Once  the  beam  was  aligned  parallel  with 
the  top  and  one  edge  of  the  optical  Ix-nch  using  a  ri'ference  poini  on  a  inii\'al)|e 
scieen.  an  .\F,('  beam  collimator  using  a  .a//m  pinhole  was  placc'd  in  tin'  beam  path. 


To  simplify  the  job  of  aligimu'nt,  the  collimator  was  bollc'ci  to  a  transverse',  line 
resolution,  movable  stage  which  was  mounted  on  a  laboratory  jack.  The  output  of 
the  collimator  was  a  2.5cm  diameter  Gaussian  [>rorile.  This  was  passed  through  a 
Wurin  aperture  aligned  on  the  reference  point.  'I’his  aperture'  matched  the  size'  e)r  the' 
fibe'f  and  was  used  strictly  to  align  the  camc'ra  and  fiber  a.sse'inbly.  d'ei  facilitate'  the' 
mounting  of  the  fibei  within  the  camera,  special  C-mount  fixture's  we're  fabrie-ated 
with  '5mm  hole's  through  the  center  allowing  a  snug  fit  of  the  fil)er  within  the'  fixture', 
'riiree  fixture's  measuring  1  -^in,  and  were  fabricated  to  facilitate  thre'e'  fibe'is 
in  h'ligths  of  In;,  2?’n,  and  'Mu  re'spe'ctively.  The  fixture's  were  thre'aeled  outside'  to 
se  ie'W  into  the  C-mount  lens  casing  of  the  camera.  This  f)re)\  ide'd  rigiel  e-enipling 
of  the  fiber  and  camera.  Each  fiber  protrueU'd  Im???  from  the  e'lid  of  the  C-me)unt 
fixture  inside  the  camera.  This  insureel  hum  se])aration  between  the  enel  e;)f  the 
fibe'r  and  the  glass  plate  covering  the  CID  detector  array.  Since  the  glass  jilate  was 
1  .'uum  thick,  tiie  distance  between  the  end  of  each  fibc'r  and  t  lu'  detector  was  2.5mm. 
lo  simulate  \‘ariovis  angh's  of  arrival(AOA),  the  camera  ass('m!)Iy  was  placed  on  a 
rotating  stage.  This  rerpiired  another  fixture  to  be  designc'd  and  fabricatt'd  to  mount 
tli('  camera  to  the  rotation  stag*'.  This  fixture  allowed  the  canu'ra  to  be  translated 
in  a  radial  direction  to  facilitate  alignment  of  tin'  <'nd  of  tiu'  fibi-r  with  the  rotational 
axis  while  kec'ping  lateral  rium'inent  to  a  miniimmi.  Tlu'  camera  and  rotation  stagi' 
was  mountx'd  to  a  t  rans\'ers<'.  fine  ix'solution.  movable'  stage'.  Spaea'i  s  we're  use'd 
to  aeljust  height  since  another  laleoratory  jack  was  not  available’.  The'  came'ra  was 
aligned  by  cejitering  the  filx’r  on  the  '5tum  Ix'am  passing  threengh  the  aix-rture.  To 
ce’nte’r  the  end  of  the  fibe’r  e)n  the’  axis  of  rotation,  the’  eame'ia  was  Inrne'el  eui(live' 
mexh’)  and  then  rotate’el  api)re)ximale’ly  25  eh'gre'e's  fre)m  beue'sight.  I  he’  came’ra  was 
the'u  translated  radially  fe>r  maximum  inte'iisity  and  le)eke'el  elown.  This  is  shown  in 
figure’  ID. 

The’  -imru  ape’rture’  was  t  he-n  re'()lae-('d  l>y  2  ne'iitral  ele-nsity  (ilte-rs  with  a  eeem- 
biiie'el  e)i)tie'al  eh'ilsit  y  ( O I ) )  eef  1.2(i.e’.  I(l'^^).  In  aeielitieen  the'  variable-  atle'iiuator  hael 
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Figure  10.  Centering  End  of  Fiber  on  Axis  of  Rotation 

to  be  adjusted  to  maximum  attenuation(OD  =  2)  to  keep  the  camera  saturation 
to  a  minimum.  Since  the  ^inni  optical  fiber  sampled  a  small  portion  of  the  2. arm 
Gau.ssian  beam,  a  uniform  amplitude  plane  wave  was  simulated.  Since  the  camera 
array  was  sensitive  to  visible  light,  and  no  filtering  was  ix-adily  available,  data  was 
takc'ii  under  conditions  of  very  low  ambient  light. 

Preparation  of  Fibers 

Before  the  fibers  could  be  used  for  (his  research,  th('  ends  had  to  be  i)olisheil  as 
flat  and  smooth  as  possible.  I'l;-'  following  |)rocedure  di-tai's  this  work.  Tlu'  ])olish('r 
us<'d  was  a  Buf^hler  Fibrmet  dual  |)laten  poiislier.  The  fiber  was  cut  ap|)roximat('ly 
oversize  by  etching  oiu'  side  of  the  fiber  with  a  razor  blade.  Tlu'  filx'r  was 
tin'll  grasped  in  both  hands  and  snaivpi'd  to  achieve'  a  rather  clean  break,  'fin'  filx'r 
was  then  inserted  into  an  adapter  |)late  mounted  on  tin'  polisher.  I’his  adapter 
plate  was  sjiecifically  fabricated  for  (he  '.Unin  filx'r  used  in  (his  research.  After  tlu' 
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nixT  was  propc'rly  inountc'cl  on  the  polisher,  a  Vlfnn  grit  alumiiiuni  oxide  polishing 
disc(yellow)  was  placed  on  the  right  platen.  These  discs  are  color  coch'd  and  range 
from  0.3/<m  to  12/071  in  surface  texture.  The  polisher  was  started  and  a  bead  of 
fibrinet  polishing  extender  was  placed  on  the  disc.  After  15  minutes,  the  course  disc 
had  removed  approximately  \inin  of  material  froin  the  end  of  the  fiber.  The  end  of 
th('  fiber  was  now  flat  but  rougli.  At  this  time  the  I2fitn  disc  was  replaced  with  a 
d//;n  silicon  carbide  disc(gray)  and  a  aluminum  oxide  disc( white)  was  placed 

on  the  left  platen  of  the  polisher.  A  small  bead  of  polishing  extender  was  placed  on 
each  disc.  Once  again  the  fiber  was  placed  against  the  right  platen,  this  time  for  only 
2  to  3  minutes.  'I'he  fiber  was  then  smooth  enough  to  be  finished  on  the  left  platen. 
This  fiinshing  step  took  anywhere  from  6  to  10  minutes.  The  adapt c'r,  with  the  liber 
still  securely  clamped,  was  then  removed  and  the  end  of  the  fiber  examined  with  a 
small  magnifying  lens.  Removing  the  adapter  and  fiber  rather  than  the  fiber  alone 
h('lped  to  insure  that  the  end  of  the  fiber  would  be  reinstalled  i)arallel  to  the  surface 
of  the  polishing  disc  if  furtiier  poiisliing  was  necessary.  Two  fibers  were  polished  for 
each  of  the  three  lengths. 

d.  /  Gathering  Template  Data 

With  the  camera  in  live  mode,  the  camera  fiber  assembly  was  rotatc.'d  to  find 
boresight(i.e.  zero  degree  .'\0A)  on  the  fiber.  'I'he  stag<?  was  then  locked  down  and 
the  angular  measurement  on  the  rotation  stage  noted.  In  this  case,  zero  di'gree  AOA 
on  the  fiber  was  0°  +  25'  on  the  rotation  stage. 

Data  from  the  3d;  fiber  was  taken  first  since  the  longt'r  (ilx'r  ('xhibit.ed  tlu' 
greatest  amount  of  change  per  degn'e  of  motion.  This  was  ('xpected  due  to  t  he  larger 
number  of  reflections  imposed  by  the  increa.sed  length.  'I'his  frame  was  exc'cuted  and 
stored  in  file  0DEG3.  J'his  r('|7resented  a  zero  degree  AOA  on  tin'  3;;;  fiber.  I  he 
stag('  was  rotated  oik'  degree  to  the  1°  i  25'  mark  on  th('  lotational  stage,  d  his 
frame  was  cxccitted  atid  .stor(’d  in  file  lDF/13.  ft  is  iiiiporlani  (lud  the  files  are 
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named  in  this  manner  with  no  extension  since  these  are  the  filenames  recognized 
b\'  the  preprocessing  software.  This  procedure  was  repeated  through  25°  resulting 
in  a  total  of  26  template  files.  Data  was  then  gathered  for  the  2m  and  lin  fibers. 
Those  files  were  labeled  *DEG2  and  '^'DEGl  respectively  where  the  *  rcpresenis 
some  arbitrary  AOA.  The  angle  of  incidence  was  the  only  parameter  allowed  to  vary 
during  the  initial  data  gathering  portion  of  this  research. 

3.5  Preprocessing 

To  help  determine  the  optimum  preprocessing,  a  program  was  written  which 
performed  various  operations  on  the  data.  These  included  median  filtering,  normal¬ 
ization  based  on  highest  intensity,  log  compression,  location  of  peak  intensity,  and 
thresholding  to  name  a  few.  This  program  called  logSxS.c  allowed  direct  manipu¬ 
lation  of  full  size  Spiricon  files.  Further  details  on  this  program  can  be  found  in 
appendix  B.  Through  the  use  of  this  program  various  preprocessing  schemes  were 
tested  and  the  results  displayed  on  the  monitors  to  observe  the  effects.  These  results 
were  used  to  design  the  final  preprocessing  program  preproc.c. 

The  program  preproc.c  first  reduced  the  512  x  512  Spiricon  data  to  a  256  x  256 
array.  This  was  done  for  two  reasons.  First,  because  of  the  small  size  of  the  fiber 
relative  to  the  sensor  array,  the  fiber  output  intensity  pattern  fit  within  the  center  256 
X  256  pixels.  Second,  and  most  importantly,  this  reduced  the  263,052  byte  Spiricon 
file  to  65,536  bytes.  This  reduced  file  was  strictly  data  witli  no  headc'r  or  paranu'tcr 
information.  Next,  the  data  was  median  filtered.  This  median  filter  differed  from  the 
filter  described  in  chapter  2  in  that  it  placed  the  median  value  into  the  center  pixel 
of  the  array  currently  being  filtered  instead  of  creating  a  new  array.  Therefore  each 
median  sort  was  performed  using  a  previously  sorted  median  value.  This  was  named 
a  recursive  median  filter.  Through  the  u.se  of  the  previously  mentioned  log3.r3.c 
program,  this  was  discovered  to  be  more  efficient  than  the  mi'dian  filti-r  described  in 
chapter  2  .  This  proce.ss  helped  remove  the  significantly  higher  intensitii's  int  roduci'd 


by  overly  sensitive  pixels  in  the  CID  camera  array.  This  op('ration  thus  removed  a 
great  deal  of  unwanted  noise  while  preserving  the  shape  of  the  pattern. 

Following  the  filtering  operation,  the  data  may  or  may  not  be  converted  to  a 
compressed  logarithmic  scale.  This  helps  increa.sc  the  low  intensities  with  little  effect 
on  the  higher  value  intensities.  This  option  was  not  .selected  during  this  research 
effort.  The  program  next  entered  a  normalization  routine  which  scaled  the  data 
such  that  the  highest  intensity  in  the  file  was  set  to  2.55.  This  took  advantage 
of  the  full  range  of  data  values  allowed  by  Spiricon.  Following  this,  Fie  data  was 
thresholded.  This  option  allowed  the  user  to  not  only  choose  a  threshold  value  but 
also  select  whether  the  data  shoidd  be  binarized.  In  other  words,  evc-rything  below 
the  threshold  value  is  set  to  zero,  but  if  the  binary  option  is  selected,  then  everything 
above  the  threshold  is  set  to  255.  Otherwise,  tho.se  values  abov('  the  threshold  are 
left  untouched. 

To  help  select  a  threshold  value,  a  histogram  option  can  be  chosen  in  the  be¬ 
ginning  of  the  program.  This  option  places  a  histogram  of  the  normalizi'd  data  for 
each  file  into  the  file  TEMPH1ST.DOC  or  TESTHIST.DOC  dej)ending  on  whetlu'r 
template  or  test  data  is  being  preprocessed.  The.se  files  also  contain  the  sum  of 
int('nsities  and  the  peak  value  before  and  after  normalization,  as  well  as  the  coor¬ 
dinates  of  the  peak  intensity.  Figure  11  illustrates  t  his  histogram,  ddu'  oul])ut  files 
are  appended  with  a  .256  extension  and  placed  in  the  PHOCDA'fA  subdirectory  of 
SPIR.  This  program  will  process  multiple  files  without  user  inter\enl  ion.  I'nrther 
information  on  this  program  can  be  found  in  a[)pendix  C. 

."7.6’  r  IULt!<OllUj 

Information  concerning  the  general  shape  of  an  image  is  cariT'd  in  the  low 
fix'quency  Fourier  components  (3:262,263).  Thus,  an  attcmipt  was  made  to  write 
an  FFT  program  to  run  on  the  PC  that  could  proc<'ss  a  256  x  256  unsignc'd  byte 
array.  Two  problems  were  encountered.  'I'he  first  involvc'd  tlu'  conversion  of  each 
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HISTOGRAM  OF  DATA 

(Value  in  parentheses  is  value  on  spiricon  plot) 
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Fillin’  11.  Histogram  of  Normalizrd  Data 


intf'nsity  value  from  one  byte  (unsigned)  to  eight  bytes(.signed),  thereby  allowing 
double  precision  floating  point  calculations  in  computing  the  hFd'.  d’his  process 
increased  the  size  of  each  file  to  524,288  bytes.  This  added  to  the  second  prol)lem  of 
limited  memory.  The  IBM  PC  allocates  Random  Access  Memory(RA.M)  in  GlKi)yte 
segments.  Significant  time  was  spent  trying  to  perform  the  FFT  on  this  large  data 
file  within  the  memory  constraints  imposc'd  by  the  PC.  Fventually  this  was  overcome, 
lint  the  results  from  the  FFT  were  erroneous.  During  the  troubleshooting,  it  was 
discovered  that  George  Letourncau  had  a  working  version  of  an  FFT  written  in  ADA 
which  ran  on  the  VAX  11/780  using  the  V^MS  operating  s\'st('in  (8).  d'his  [irogram 
accepted  a  256  x  256  array,  unsigned  byte  format,  and  n'tunu'd  a  256  x  256  unsigiu'd 
byte  array  representing  the  magnitude  of  the  FFT.  This  iirogram  performed  all  the 
necessary  quadrant  swapping  to  center  the  frequency  com/ionenis. 

To  use  this  program,  transfer  the  data  files  from  the  PC  to  the  mainframe  in 
the  following  manner.  First  copy  the  data  files  to  floppy  disks.  Take  these  disks  to 
room  243  and  install  in  drive  B  of  a  PC  containing  the  File  I'ransfer  Program(  F'f  P). 
Drive  B  will  handle  double  density  and  high  density  disks,  fype  F  f P  followed  by 
the  name  of  the  computer  which  has  the  FF'F  programfi.e.  k'd'P  i780a).  .\t  tlu'  FTP 
prompt  type  BINAlTl'  since  data  and  not  text  is  being  t ransfern'd.  .Now  type  .MPF'f 
b:’'  .256  for  a  multiple  transfer,  and  all  files  will  be  transferred  to  the  mainfraiiK'. 
Once  transfer  is  complete,  type  BYE  to  return  to  DOS.  Logon  to  tin'  mainframe 
from  an  appropriate  terminal  and  type  RUN  2DFT.  Type  in  the  filename  of  the  file 
to  be  processed  at  the  prompt.  Upon  comjiletion  of  the  f  F'!'.  the  program  will  ask 
for  the  name  of  an  output  file.  'J’y|)e  the  original  filename  with  an  f'd’  predix.  for 
example,  if  the  input  file  is  0DEG3.256,  the  out])ut  file  will  be  F'fODFt  13.256.  After 
tin'  filename  is  selected,  the  |)rogram  will  continue  to  iK'iform  a  Hough  transform, 
so  type  CTRL  Y  to  exit  the  program  and  restart  it  for  the  next  file.  .After  all  files 
have-  been  Fourier  transformed,  transfer  the  files  back  to  drive  B  on  tin-  PC  by  using 
MGF/f  instead  of  .MPUT  at  the  F'fP  prompt. 


When  this  program  was  first  used,  some  unexpected  results  w('re  obtained. 
Most  of  the  files  had  outputs  of  zero  magnitude.  It  appeared  that  the  FFT  program 
r('(iuired  binary  input  since  the  only  files  which  yielded  expected  result.s  were  tho.se 
with  intensities  of  255.  All  files  should  have  been  normalized  to  255,  but  a  roumling 
error  resulted  in  a  normalized  peak  value  of  254.  However,  five  of  the  26  template 
files  had  peak  intensities  of  255  when  the  frame  was  executed  and  therefore  did 
not  require  normalization.  The  rounding  error  was  corrected,  and  the  files  wc-re 
re|)rocessed.  This  time  the  files  were  binarized  to  insure  the  entire  image  was  Fourier 
transformed.  The  threshold  was  set  at  64  to  reduce  noise.  Tlu'se  files  yi('lded  good 
results.  After  further  iiux'stigation  it  was  confirmed  that  the  Fk' F  program  did  not 
r('cognize  values  below  255.  'I'lie  FF'F  program  was  correct<'d  and  bettc'r  rc'sults  were 
obtained.  T’lese  results,  as  well  as  problems  with  the  noise  reduction  threshold,  are 
discussed  in  chapters  4  and  5. 

To  extract  the  fundamental  and  first  three  harmonics  of  the  F1''T.  a  program 
was  written  called  template. c.  This  program  accepted  input  files  with  the  FT  prefix, 
riie  files  were  read  in  one  at  a  time  and  values  in  the  center  7x7  array,  which 
represented  the  desired  low  frequency  components,  were  written  into  each  row  of  t  he 
file  d’FMPLATF.DAT  in  the  TE.MF^LATE  directory.  This  result etl  in  a  temi)late  file 
of  26  rows  and  49  columns.  Each  row  represented  a  49  com])onent  feature  v('ctor 
containing  the  shape  information  of  the  original  image.  A  program  called  dist(inc(  .c 
was  written  to  calculate  the  separation  between  each  of  the  26  vectors  in  49- D  sjmee 
to  insure  adequate  distance  betw('en  the  vectors,  thus  allowing  K'ss  chance  of  error 
when  classifying  a  noi.sy  test  vector,  ddie  distance  between  each  template  vector  and 
its  nearest  neighbor  was  placed  itt  the  file  rE.MFDlS  IfiDOC.  luirtlier  information 
on  these  two  programs  can  be  found  in  ap|)endix  D  and  ai)pendix  E  res])ect iv('ly. 


7  Classification  of  Test  Vectors 

Gathering  lest  data  was  the  next  step.  The  camera  was  rotated  to  an  arhilraiy 
angle  and  locked  down.  The  AOA  was  taken  from  the  rotational  stage  atul  written 
down  for  future  reference.  This  frame  was  executed  and  stored  in  lih-  0  I  S  F*.  I  Ik'  0 
in  this  case  represented  the  0  number  test  vector,  not  the  number  of  degrees,  since 
our  objective  was  to  classify  this  file  using  the  previously  created  template  file  to 
determine  the  AOA.  The  *  represented  the  fiber  length,  similar  to  tlu'  pri'vioiisly 
mentioned  template  file.s.  This  proce.ss  was  repeated  8  more  times  resulting  in  t)  test 
vectors  numbered  0  through  8.  d'hese  files  were  processed  in  the  same  manner  as 
th('  template  files  except  for  a  minor  ilifferi'nce.  This  time  the  tfinplatr.c  program 
created  a  TEST.D.AF  file  containing  the  !)  test  vectors  and  th('  di.stann.r  |)rogram 
calculated  the  distance's  between  the  template  vectors  and  ('ach  te'st  vector,  flx'se' 
distances  were  stored  in  TF^STDIST.DOC.  The  minimum  distance  classilie'd  the  test 
v('ctor  ;\O.A.  I'he  results  are  pix'sented  in  Chapter  4. 

J.S  Software  X'erification  Procedures 

A  significant  portion  of  the  time  invested  in  this  researcli  was  dedicalcfl  to 
software.  This  inclndefl  not  only  programming  but  also  verifying  correct  results. 
‘C  is  not  a  very  restrictive  language,  and  many  incorrectly  written  programs  will 
still  compile,  producing  incorrect  results.  Thus,  verifying  correct  output  at  various 
stages  of  programming  was  necessary.  I’he  first  jirogram  written,  lo<)S.r3.c^  contained 
se\-eral  routines  which  were  individually  verified.  Verification  also  took  place  anytinu' 
a  change  was  made  to  the  jirogram. 

Programming  the  median  filti'r  was  a  little  tricky  since'  the  data  arrays  nse'd 
in  this  tlu'sis  elfort  we're  to  large  te)  be  sterre'el  in  the  PC  me-meirv.  Fo  ove'ree)m(' 
this  [)re)blem,  the  elata  was  se'gme'iite'd  anel  ru)  more  than  U)t)  reiws  wvvr  reael  in 
at  enie’  time.  In  aider  to  correct}}'  [X’tform  the  me'elian  filtering  ope'ratiern,  the'se' 
segments  had  to  overlaj).  A  file'  naiiu’el  PLAY  cerntaining  e-ernse'cnt  i\  e'ly  incre'asing 


integers  was  created  to  verify  tlie  correct  operation  of  the  median  filt('r.  es|)('cially 
at  tfie  seams  where  tlie  segments  overlapped.  Since  a  median  filter  will  not  affect 
consecutively  increasing  numbers,  a  correctly  operating  median  filtc'r  would  produce 
no  change  in  this  file.  Once  this  small  16  x  16  fde  demonstrated  correct,  opeuation  of 
the  median  fdter,  the  routine  was  changed  to  handle  the  .512  x  -512  Spiricon  arrays. 
'The  output  of  this  larger  array  was  then  displayed  on  the  Spiricon  real-time  monitor 
and  compared  to  the  original  image  to  verify  correct  filtc-ring  of  the'  largt'r  array. 
For  further  verification,  a  binary  image  of  a  circle  and  a  scpiare  were  also  cn-ated 
inside  a  .512  x  512  array.  A  median  filter  will  n<<t  afr{'ct  st('j)  change's  in  intensity. 
Comparison  before  and  after  processing  showed  no  change  in  either  of  these  two 
arrays.  This  routine  also  calculated  the*  sum  of  the  intensities.  This  sum  was  also 
used  as  a  reference  to  verify  any  change's  late'r  maeh'  in  the'  pre>gram. 

The  next  part  of  the  jrrogram  was  the  ne)rmalizat ion  |ne)ce'ehire'.  To  ve'iifv 
oi)eration  of  this  routine,  scaling  factors  and  rande)m  pixe'l  int e'lisil  ie's  befe)re  and 
afte'r  normalization  were  printe'el  to  the  scre'e'ii  and  comi)areel  to  hanel  calculate'el 
re'sults.  The  log  compre'ssion  anel  the'  thresholding  routine  we're'  hanelh'el  in  the'  same 
manner. 

,\11  routine's  in  lofiJx-i.c.  including  one's  ne)t  me'ntione'el  abe)Ve'.  we're'  adelit  ie)nally 
ve'iifie'd  using  .Norton’s  IIF.X  ('elite)r.  Re'sidts  from  we'ie'  use'ei  te)  he'l])  ve'iify 

preegrams  written  later. 

d'he  2dft  program  was  verifie-el  by  creating  .'5  files  ce)ntaining  a  re'ct  functie)n, 
eh'ha  function,  anel  a  constant  funetie)n.  'I'lie'  outj)ut  ])roeluce'd  the  e'\|)e'cte'el  sine 
functie)n,  eonstant  function,  and  elelta  function  respect  ive'lv  (2;.'UT).  Ne)rte)n‘s  llF.X 
e'elileu'  was  once'  again  use'ei  to  furthe-r  ve-rify  re'sults. 

Se'veral  pre)grams  we're'  writte'ii  in  ‘C  te)  ve'iify  the'  ejpe'iat  ieen  e)f  the'  Spiricon 
seiltware'.  for  example’,  using  a  pre)gram  e'alh'el  oddrow.c  whieh  se-t  odd  row  pixe'ls 
to  255  and  even  row  pixels  ('epial  to  zero,  it  was  disco\e're'el  that  Si)iricun  only  elis- 
plays  e'\e'ry  e)ther  re)W  of  eame'ia  eiata  ein  the'  re'al-time'  meanite)!'.  .Aneelhe'r  program 


calk'd  (lot.c  set  a  single  to  255  and  all  otlu'rs  to  zero.  Spiricon  is  suppcjsed 

access  a  us<’r  created  file  named  BADPIX  to  ignoi'f'  tliis  pi.xel  while  it  processes  tin' 
data,  riirough  the  use  of  dot.(\  it  was  discovered  that  this  k'ature  was  not  working 
correctly.  See  appendix  A  for  further  ii\forniation. 

This  chapter  examin('<l  the  (h'sign  of  theo[)tical  heiu  li.  I  his  inclutied  a  discus¬ 
sion  of  l)('am  aiignuK'ut  and  ojitical  comi)on('nts  used.  This  was  followed  hy  remarks 
on  filx'r  prc'parat ion.  data  gathering,  |)rocessing,  and  software'  ve'iificat ion  |)roce- 
duK's.  Idle  next  chapter  pix'sents  the  results  of  this  research  and  a  brief  analysis  of 
tliose'  re'sults. 


IV.  Data  Analysis 


The  previous  chapter  discussed  the  procedures  used  to  gattu'r  and  i)ro( css  the 
data.  I'his  chaoter  analyzes  the  results  ohtaiiu'd  from  this  rc'search.  1  lu'  non 
symmetry  observed  in  tlie  output  pattern  for  similar  angles  on  opposite  sides  of 
horesight  is  first  examined.  'Fhis  is  followed  by  a  presentation  of  tlu'  various  data 
gathered  and  a  discussion  of  tlu'  results. 

/f.l  Di.^cussion  of  Non  syntnu  trtj 

During  tlu'  data  gat  Inning  i)ortion  of  this  reseain  h,  a  non  symmetry  was  ol)- 
s('rved  between  opposite'  sides  of  boresight.  This  contradicts  the  belief  that  tin* 
cult  put  [)attei  n  of  a  circularly  symmc't  rical  fiber  should  also  be  symmc't  rical.  I  his 
anomaly  can  bc'  ('asily  explainc'd.  Whc'u  tin'  c'uds  of  the  libc'r  arc'  |)olish('d.  extrc'ine 
care  must  be  taken  to  polish  the  ends  orthogonal  to  tin'  core-  axis  of  the  fiber.  I  he 
fibc'rs  used  in  this  research  were  clampc'd  in  a  fixture  labricated  spc'cilically  for  a 
3;nn;  fibc'r.  '1  In'  ])olishing  |>rocc'ss  was  completc'd  for  c'ach  fiber  In'fore  it  was  re¬ 
moved  from  tin'  fixture.  Due'  to  the  inability  to  consistently  clam])  the  fibc'r  in  tin' 
fixture,  tolerances  in  the  |)olishing  machine  com])Oin'nt s.  such  as  bc'arings,  and  tin' 
inability  to  pc^lisfi  tin'  c'lids  pc'ifeci  ly  flat,  tin'  c'ntrancc'  and  exit  of  the  fiber  may 
neither  bc'  jiarallel  to  c'ach  other  nor  normal  to  tin'  core'  axis,  d  her''for('.  wln'ii  tin' 
libc'r  is  mountc'd  in  tin'  camera  and  the'  core'  axis  is  ccdlinc'ar  with  tin'  laser  bc'am. 
tin'  beam  may  contact  the  entrance  to  tin'  fiber  at  an  angle  other  than  0°.  1  he 

output  t  hc'rc'lorc'  apjn'ars  non  symmc'trical  as  tin'  camc'ra  is  rotatc'd  from  cun'  side'  of 
boresight  to  the'  other.  I  his  can  bc'  clemonst ratc'd  by  rotating  tin'  camc'ra  while'  in 
li\c  meulf'  and  watching  the  output  i)at.tc'rn  change'.  .Also  tolerances  in  the  threads 
ol  the  ('  mount  fixture'  and  the-  lac't  that  the-  hole'  may  ha\e'  bc'eii  drilled  at  a  slight 
angle'  within  the'  ('  meuinl  fixture'  add  to  tin-  non  s\nnne'tr\’.  riiiis  it  is  im|)c)rtant 


tliat  t()l('raiic'('s  hr  kc'pt  to  a  luinimuin  with  n'six-ci  to  tlu'  |)C)lislii!i,';  and  inonut  ing 
ot  the  fiber  in  a  practical  apj)licat ion. 

J  Pri srnlntion  of  Rrsulls 

Data  was  j^atlicrcd  for  all  tlircc  knigtlis  of  fiber,  d'his  consisted  of  grabbing 
a  Iraine  of  data  for  each  degree  of  incidence  from  0°  throngh  '2o° .  bach  franu'  was 
|)lac('d  in  a  file  and  (-acli  fih'  i)rocesse(l  to  create  a  If)  com[)on('nt  temi)lal  e  \-ect or  w  hich 
represented  that  part  icnlar  angle  of  arrivalf  AOA).  All  templatf'  vectors  w('re  i)laced  in 
the  fih'  1  E.MPLA  ri’..l)A  r.  As  stated  earlic'r.  the  fbn  fiber  was  ])ro<'C'ssefl  first  because 
of  th('mor('  ra|)id  ratc'of  change  of  output  pattern  p('r  degre<' of  change  in  .\0.\  when 
compaia'd  to  tin'  shorter  fibens.  Fliis  was  due  to  the  incn'ast'd  number  of  la'flections 
encount('r('d  by  tin-  laser  light  as  it  propagati'd  throngh  the  libc-r.  'rhr('('  diffc-ia'iit 
processing  l('chid<pies  weic  us<'d  for  t  his  filx'r.  Sinc(' t  he  fast  I'onriei' t  I'ansformf  1  'b  I' ) 
routine  initially  repnirt'd  binary  iiijnit.  a  threshold  of  tilfrangi'  b  --  was  set  in  the 

l)reproc('ssing  |)rogram  /naprec.cand  the  binary  option  chosc-n.  d'his  tlirc'sliold  was 
chosc'ii  bas('d  on  comparison  of  histograms  of  all  t('mplat<'  files.  I'est  data  was  then 
gathered  for  various  \OA  and  process<'d  in  the  same  maniK'r.  d'hese  M).\  are  shown 
in  colnmn  1  of  tabh'  1.  The  two  7°  ineasurenK’iit s  w('r('  takc-n  separately  to  check 
re|)eat ability,  liie  bnclidean  distance  Ix-twec'ii  each  t('st  V('ctor  and  all  template 
\eciors  was  caicnialed.  I  he  minimum  distanc('  classified  each  t('st  vector.  I  he 
residts  ar<’  shown  in  the  two  columns  labeh'd  Hl.X.A Rd’  ((>  1 )  of  t abh'  1 .  I'lu'  nnmber  in 
jiarent  liesc's  rf'prt'sent  s  t  (le  t  lire's  hold  vahu'  cliose'ii.  ih('  entry  in  (  he  drsi  a  nee  column 
repi'cse'iits  the  distance'  b('tw('('n  th(' tc'st  ve'ctor  anel  its  e  lose'st  te'inplate'  xa'e  teir.  It  is 
pi  i'se'iite'il  he're'  as  a  re'lat  i\'e'  e'omparise)!!  rat  he'r  t  han  abseilnte'  me'asnre'iiu'nt .  1  he'  ’  in 
the  estimate'd  ,\(),A  e'olnnin  iinlie'ale's  cori'e'ct  e  lassifie’al  ion  eif  tli  e  ti'st  \'('clor.  lour 
te’^l  \('cl()rs  we're  classilie'el  cotre'i’tly  but  e'rrors  as  high  as  (i  de'grees  wi'ii'  ap|iare'nl. 

It  was  bi'li<’\'e'd  that  senile'  inleermal leni  was  lost  wlii'ii  the  da'a  was  binarize'd 
so  the'  h  h  I  ineigram  was  mealifie'il  te»  ae-e'e'pt  iiein  binary  data.  I  bis  nein  binar\' 


Tal)le  1.  Classification  of  3/7?  Fiber 


Distance  Between  Test  and  Template  \  ectors 
3m  Filjer  at  0°  Reference 


Binary (64) 

Non- Binary 

Non-Binary(64) 

Act. 

AOA 

Est. 

AOA 

Distance 

Est. 

AOA 

Distance 

17 

*17 

29.02 

16 

12.96 

■0 

22 

20 

9.60 

21 

24 

21 

22 

6.78 

22 

10.49 

15 

12 

68.09 

*15 

■a 

21 

5.48 

20 

3 

158.93 

*3 

44.36 

7 

*7 

154.40 

5 

40.80 

*7 

160.70 

1 

13 

91.86 

/ 

35.95 

13 

135.36 

11 

1 

144.35 

*11 

21.12 

1 

110.28 

44A  corixH't 
W’CE  =  6" 

44%  correct 
WOE  =  2° 

33%  corrc'ct 
W'CE  =  6° 

33 


'( 'orrert  Classiticat  ion 
W'Cl'i  z=  Worst  Case  I'irror 


Table  2.  Classification  of  Im  and  2in  Fii)ers 


Distance  Between  Test  and  Template  Vectors 

2m  Fiber 

Izn  Fiber 

0°  Ref. 

0°  -  15'  Ref. 

0°  +  15'  Ref. 

Act. 

auA 

Est. 

AOA 

Distance 

Est. 

AuA 

Distance 

Est. 

AOA 

Distance 

17 

10 

23.41 

16 

27.53 

18 

26.19 

22 

17 

21.82 

*22 

22.58 

23 

17.66 

24 

*24 

16.97 

*24 

17.55 

25 

15.49 

15 

*15 

28.71 

*15 

33.23 

*15 

35.76 

21 

22 

21.86 

*21 

30.03 

22 

21.07 

3 

*3 

101.54 

*3 

38.65 

*3 

37.95 

7 

*7 

51.00 

6 

72.90 

*7 

56.66 

7 

*7 

46.95 

*7 

75.17 

*7 

59.00 

11 

*11 

43.47 

*11 

43.17 

12 

38.94 

67%  correct 
WCE  =  r 

78%  correct 
WCE  =  r 

44%  correct 
WCE  =  1° 

'orrecl  Classificat  ion 
WCK  =  Worst  Case  Error 

l)rocessin<:f  was  ])('rforined  with  and  without  a  tliresliold.  'I’Ik'  results  are  shown  in 
the  last  4  columns  of  table  L.  Fhe  number  (Gl)  rc'presents  the  thrc'shold  value.  It 
can  b('  s('('n  that  the  non  thresholded  processing  yif'lded  the  Ix-st  results  with  4  I'X 
(■<M  i<'ct  classification  and  a  worst  case  errorf W('l'y)  of  only  2°. 

d'he  ne.xt  research  ('ffort  examined  the  effect  of  fiber  length  on  AOA  defection. 
All  data  was  f)rocessed  non-binary  with  no  threshold,  based  on  tin'  prexious  results, 
lest  Data  was  first  gathered  for  the  2n/  fiber.  The  0°  (borc'sight )  reference  was  found 
to  be  0°  —  15'  on  the  rotational  stage,  'rids  was  40'  h'ss  than  wlnui  the  teunplatc' 
data  was  takfui  t  lu’  month  belore.  This  can  easily  lx-  account (xl  lor.  Whenever 
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the  camera  is  loosened  from  its  movmt  and  allowed  to  mov'e  radially  to  adjust  for 
different  fiber  lengths,  it  may  also  move  very  slightly  in  a  lateral  direction.  Another, 
more  significant,  factor  concerns  the  alignment  of  the  laser  beam.  In  the  month 
between  when  the  template  data  was  taken  and  the  test  data  was  taken,  the  beam 
steering  instrument  was  accidentally  adjusted.  This  required  complete  realignmc'ut 
of  the  beam.  In  addition,  the  0°  reference  is  determined  by  rotating  the  camera 
in  live  mode  until  a  full  intensity  circular  pattern  appears  on  the  screen.  This  is 
a  subjective  measurement  and  not  accurate  within  minutes  of  a  degree.  Taking  all 
tlic'se  factors  into  consideration,  it  was  not  surprising  that  boresight  had  changed 
from  when  the  template  data  had  been  taken.  Measurements  were  taken  at  the 
same  .AOA  increments  as  the  fiber.  Results  were  favorable  for  the  ‘lin  fiber  with 
GT'T  correct  classifications,  how('vcr.  errors  were  significant  wit  h  a  W'CIC  of  7°.  These 
results  are  shown  in  table  2.  All  angular  references  in  the  tables  arc'  relative  to  actual 
l)oresight  and  have  no  reference  to  the  rotational  stage. 

4-2. 1  Rrlinbilify  Tcftting  Immediately  following  this,  the  camera  was  realigned 
and  measurements  wore  made  on  the  hn  fiber.  Three  tests  were  made  with  the  \in 
filx'r  to  assess  reliability  in  determination  of  AOA.  The  first  two  tests  placed  the  test 
vectors  15'  on  either  side  of  the  0°  boresight  reference,  d'his  was  done  to  verify  the 
ability  of  the  system  to  continue  to  correctly  determine  AOA  with  a  small  angular 
deviation.  'I  hese  results  are  also  presented  in  table  2.  d’hese  two  tests  were  suc¬ 
cessful.  both  resulting  in  a  WCM  of  1°,  however  the  first  test  pc'rfoiined  31%  better, 
achie''ing  78%  correct  classification.  The  data  shows  that  all  c'rrors  in  the  first  tc'st 
were  on  the  low  sirle  and  all  errors  in  the  .second  test  were  high.  Thc'se  results  were 
e.xpectc'd  since'  the  data  was  first  takc'ii  to  the  negative  side  of  zero  and  then  shifted 
to  the'  positive  side  of  0°.  d'hese  rc'sults  are  also  presentc'd  in  table  2. 

4  2.2  of  1/  aud  Noo  Syminrhi/  1'he  third  relialulity 

test  with  the  If??  fiber  was  combinc'd  with  two  other  tc'sts.  d'hc'  first  tc'st  examiiu'd 


Table  3.  Classification  of  1?72  Fiber  With  Rotation  and  Variable  Intensity 


Distance  Between  Test  and  Template  Vectors 

Im  Fiber  at  0°  Reference 

Std.  Int. 

2  X  Std.  Int. 

Rotated  90° 
Std.  Int. 

Est. 

AOA 

m 

Dist. 

17 

*17 

40.02 

0 

47.75 

16 

3 

*22 

*22 

55.48 

*22 

l^i■ 

mam 

24 

*24 

16.12 

*24 

40.27 

25 

36.72 

22 

43.24 

15 

*15 

56.73 

*15 

80.51 

21 

*21 

36.22 

*21 

55.46 

10 

57.43 

3 

4 

60.23 

4 

104.39 

88.84 

7 

*7 

91.31 

8 

59.93 

6 

146.80 

7 

*7 

75.19 

*7 

94.48 

8 

63.07 

6 

147.08 

11 

12 

44.25 

12 

61.85 

23 

46.84 

12 

96.21 

78%  correct 
WCE  =  r 

67%  correct 
WCE  =  17° 

22%  correct 
WCE  =  12° 

11%  correct 
WCE  -  14° 

*  Correct  Classification 
V\  (’E  =  Worst  Case  Error 
Std.  Int.  =  Standard  Intensity 


the  affect  that  a  variation  in  incoming  intensity  would  have  on  AOA  classification. 
The  second  measurement  involved  rotating  the  fiber  90°  to  verify  that  non-  symmetry 
due  to  the  end  face  being  polished  non-normal  to  the  core  axis  would  affect  correct 
AOA  determination.  The  benefits  of  this  approach  were  two-fold.  First,  this  al¬ 
lowed  one  more  check  on  reliability.  More  importantly,  however,  tlu-  reliability  test 
provided  a  reference  as  to  how  good  the  result*^  of  the  intensity  and  non-symmetry 
test  really  were.  The  data  was  gathered  in  the  following  manner.  .Measurements 
were  taken  with  a  0  degree  borsight  reference  such  that  there  was  no  angular  devi¬ 
ation.  The  camera  was  rotated  to  the  first  test  angle  of  17°  and  the  data  gathered 
as  before.  Then,  the  .40  OD  neutral  density  filter  was  removed,  more  than  doubling 
the  intensity(i.e.  10'*°  =  2.51),  and  another  frame  was  grabbed.  The  .40  OD  filter 
was  replaced  and  another  .40  OD  filter  added,  thus  decreasing  the  intensity,  and 
data  gathered  again.  Finally,  the  added  .40  OD  filter  was  removed,  the  fiber  and 
C-rnount  threaded  fixture  rotated  90°,  data  gathered,  and  the  fixture  rotated  back 
to  its  original  setting.  This  procedure  was  repeated  for  each  of  the  nine  test  angles. 

The  results  are  presented  in  table  3.  The  first  column  once  again  shows  actual 
AOA.  The  next  two  columns  represent  the  third  reliability  test.  This  test  achieved  a 
coi  rect  classification  of  78%  with  a  WCE  of  only  1°.  Thus  the  reliability  of  measure¬ 
ment  of  AOA  using  an  optical  fiber  was  verified.  Using  this  as  a  comparison,  it  can 
be  seen  from  the  other  columns  in  the  table  that  both  a  change  in  intensity  and  the 
non-orthogonality  of  the  fiber  entrance  with  the  core  axis  affect  AOA  determination, 
riie  adverse  effect  of  intensity  in  AOA  detection  is  due  to  normalization  based  on 
p('ak  value.  This  effect  could  be  eliminated  through  normalization  based  on  total 
eiK’igy.  See  Chapter  5  for  further  comments. 

1  he  {in  fiber  api)eared  to  give  the  best  restilts  but  this  may  l)e  misleading. 
'File  fibers  may  have  gotten  rotated  inside  their  fixtures  between  tcunplate  aiid  t('st 
v('etor  data  gathering.  In  addition,  the  measurements  performed  in  this  research 
were  subject  to  human  error  each  time  the  beam  was  realigned  or  the  camera  ro- 
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Table  4.  Total  Intensities  and  Peak  Values  for  Ternplah'  Files 


Template  Intensities  Before  Normalization 


AOA 


0 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


3m  Fiber 


Total 

Intensity 


1282365 


1340016 


1490947 


1510842 


1520684 


1507261 


1658921 


1642212 


1608812 


1569449 


1513587 


1515596 


1468427 


1456887 


1439063 


1467542 


1409270 


1492634 


1411670 


1425437 


1285504 


1342986 

1130085 

1150165 

1076821 


2m  Fiber 

Total 

Intensity 

Peak 

1358228 

255 

1331934 

255 

1365646 

255 

1354307 

255 

1378817 

255 

1357861 

255 

1362819 

255 

1393817 

255 

1371979 

255 

1414075 

99 

1320108 


1301680 


1286716 


1280034 


1254365 


1247458 


1226093 


1259666 


1240201 

1216688 

1181791 

1134245 


u 


Im  Fiber 


Total 

Intensity  Peak 


1359795 


1363492 

255 

1344498 

1368844 

255 

1301121 

255 

1298955 

255 

1345798 

255 

BfeWAfdiM 

255 

1261141 

251 

1222484 

255 

1282444 

255 

1199557 

255 

1190632 

255 

1160464 

1157413 

1105218 

143 

1124846 

98 

1055813 

61 

1013414 

51 

Average  Noise  Intensity  =  1t)2,8-l() 
A\erage  Peak  Noise  =  12 
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tated.  The  longer  fibers  may  be  more  sensitive  to  these  slight  difFerences  which  can 
occur  between  tests.  It  is  possible  that  the  longer  fibers  may  outj)erform  the  \/n 
filler  in  more  controlled,  less  subjective  conditions.  This  may  also  hold  true  in  higln'r 
resolution  applications  where  template  data  is  taken  in  increments  of  less  than  1°. 
However,  a  look  at  table  4  does  indicate  that  the  lin  fiber  maintains  a  substantial 
peak  value  for  higher  AOA’s  which  may  play  a  part  in  its  apparent  better  perfor¬ 
mance  in  AOA  classification.  Further  research  will  have  to  be  conducted  in  this  arc'a 
before  this  can  he  validated. 

4.2.3  Discussion  of  TliiTsholding  Scherne  booking  at  table  1.  AOA  detection 
was  more  subject  to  error  when  thresholding  was  set  at  G4.  I  his  number  appc'ars 
to  be  sufficiently  low  to  eliminate  noise  while  still  maintaining  enough  signal  to 
determine  AOA.  The  reason  for  the  poorer  performance  can  be  attributed  to  two 
factors.  As  the  AOA  increases,  more  light  is  reflected  off  the  entrance  of  the  filx’r. 
d  ims  signal  intensity  is  reduced  while  noise  remains  relatively  constant,  ddiis  can 
b('  seen  in  table  4.  The  first  column  of  each  section  shows  the  total  intensity  and 
the  second  column  shows  the  peak  value  before  normalization.  Fach  total  intensity 
value  represents  the  sum  of  the  unsigned  bytes(0  —  255)  in  the  reduced  256  .\  25G 
array,  dhis  summation  was  performed  in  the  program  pnpioc.c.  do  give  some 
meaning  to  the.se  numbers,  five  reference  frames  of  data  were  taken  with  no  light 
on  the  camera  to  get  some  idea  of  the  noi.se  level  jiri'sent.  d  hese  frames  yiidded  an 
aierage  intensity  of  402,840  and  an  average  peak  value  of  12  bi'fore  normalization, 
for  AOA’s  above  20°  the  signal-to-noi.se  ratio  was  low  (uiough  to  cause  probh-ms  in 
.\()A  clas.sification.  i-dgures  12  and  13  show  how  the  higher  peak  vahu'  patterns  at 
lower  AOA's  a[)pear  before  and  after  normalization  as  compart'd  to  the  lower  ja'ak 
value  higher  .AOA’s  presented  in  figures  M  and  15.  d  he  second  factoi'  is  program 
rt'lated.  ddie  thresholding  routine  in  the  program  prrpi.jc.r  was  originally  placed 
in  IIk'  [uogram  to  binarize  the  data  or  to  .select  the  flata  abine  the  thi'eshold  for 
viewing  or  processing.  It  was  not  placed  in  tlu'  program  to  helj)  eliminate  noise. 


Figure  14.  83  Peak  Value  and  16°  AOA  Before  Normalization 


1  li('refoie,  it  was  placed  at  the  end  of  the  program,  after  the  normalization  routine. 
A  thresholding  routine  to  reduce  noise  should  he  added  to  the  (rrogram  aln'ad  ol  tlu- 
normalization  routine.  This  problem  is  addressed  in  the  recommendations  section  of 
cliapter  5. 

This  chapter  presented  the  data,  gathered  in  this  research  and  discussc'd  the 
results,  rhe  next  chapter  presents  conclusions  of  this  research  and  reccMnnu'ndations 
lor  future  research  in  this  area. 


V.  Recommendations  For  Future  Research 


I  Conclusions 

riiis  research  was  conducted  to  test  the  feasibility  of  using  tlie  output  patterns 
from  ojjtical  fibers  to  determine  angle  of  arrival( AOA ).  The  experimental  optical 
system  was  designed  and  assembled,  complete  with  camera  and  framegrabber.  Parts 
wei'e  lubricated  specifically  for  this  research  effort  to  allow  for  c'ase  and  reliai)ility 
during  the  data  gathering  portion  of  this  thesis,  d  hese  parts  included  threaded  ('- 
mount  camera  fixtures  designed  to  hold  the  fibers  firmly  in  place  without  damage, 
a  fixture  mounted  between  the  camera  and  rotating  stage  which  allowcal  only  radial 
movement  of  the  camera,  and  a  3mm  aperture  used  to  align  the  end  of  the  fiber 
with  the  axis  of  rotation.  Much  time  was  spent  learning  the  ‘lurbo  C"  language' 
and  writing  the  code  which  performed  the  signal  processing.  This  was  original  code' 
written  to  be'  compatible  with  the  Spiricon  system.  'Phis  eoele*  is  presente'el  in  the' 
appendices. 

.Several  are'as  affecting  AOA  detection  were  investigated  in  this  the'sis.  'I’lie'se' 
ine  lude'd  various  fiber  lengths,  laser  intensities,  and  processing  scheme-s.  J'he'  l/n 
fiber  gaw'  t)i<'  best  re'sult.s  althe)ugh  variatieuis  in  pre'processing  sehe'me's  or  higlie'r 
re'se)lut iejii  applications  may  affe'ct  this.  Variatie)n  in  sonree'  inlensily  le-snlte'd  in 
e'l  i eene'enis  .AO.A  classificat ieni,  ele'monst rating  the  m^e'e:!  to  moelifv  the-  pre)gram  //o - 
pror.c  to  cornpe'nsale  for  this  variation.  All  data  was  prese-nte'e!  as  an  aiel  to  lulluw  eni 
re'se'a  reli. 

.\  fe'W  e)f  tin'  sucex'.s.ses  acliie'Ve'el  in  this  tiu'sis  were'  ])rese'nte'd  abene'.  but  these' 
aie-  mine)!'  ce)mi)ared  te)  the'  major  sue-ce'ss  of  showing  that  angular  re'seiliit  iej!i  of  |° 
eer  le'ss  is  possible’  thre)ngh  the  analysis  e)f  the'  e)!itput  patte'rns  freem  an  optieml  fibe'r. 
Hee'e)mme'nelat  ie)ns  fe)r  future'  re'se’are'h  are'  |>rese’nte'el  in  the'  ne'Xt  si’cliou. 
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5 .  J  R(  CO  m  m  e  n  dation  s 

riiis  section  discusses  reconinuMidal ions  for  future  r('S(‘arcli  in  AOA  delect  iou. 
Hecoinnieudat ions  iucludc  upgrade's  to  the  optical  Ix'uch  se'tup.  correcting  for  \ari- 
ahle  intensity,  as  well  as  implementing  correlation  and  neural  network  processing. 

A  .Newport  Equipment  Corporation  (NEC)  model  l81  manual  rotational  stage 
was  used  for  this  research.  'Phis  rotational  stage  has  a  vernier  scaled  in  minutes  of 
degree  of  arc.  The  resolution  of  this  stage  was  adequate  for  this  initial  If-asihility 
study,  but  this  stage  is  small  and  can  support  the  weight  of  the  (’ll)  l  liCll  camera 
and  cables  onl\'  if  t  hat  weight  is  ev(>id\'  (list  ribu!('d  across  1  he  st  age.  .Since  t  he  camera 
and  cabh's  are  off  center  from  •?  to  6  inclu'S.  this  stage'  allows  tlu'  camera  assembly 
to  drop  at  a  slight  angle  such  that  it  is  no  longe'r  jraralh-l  with  the  top  ul  (he  opiical 
bench.  Since  the  laser  beam  is  aligiu'd  parallel  to  the  toj)  of  the*  Ix'iicli.  llu'  beam 
encounte'is  the  entrance  to  tlu'  fiber  at  a  slight  angle',  'flie'  use  of  the  NEC  model 
■19o  rotational  stage  is  suggested,  d'his  stage  is  larger  and  has  almost  thice'  time's  the 
load  capacity  of  the'  mode'l  -181.  d'his  is  also  a  motorize'd  stage  controlled  by  an  IHM 
or  compatible'  PC  using  the  .MS-DOS  operating  system,  d'his  added  feature'  allows 
stage'  rotation  in  incre'inents  as  small  as  .001  de'gre'e  while  automatically  ael justing 
for  backlash  in  the  mechanism,  d’he  NEC  495  should  make  (he  job  of  gatliering  data 
less  tedieuis  and  also  allow  rejx'atable  pe'rformance.  1  his  stage'  was  on  orde'r  at  the' 
time'  this  re'se'areh  was  jie'rforme'd  and  shexilel  arrive  in  De'cember  1989. 

Limited  time  evas  spe'ul  exi  pr('pre)e'e'ssing  sche'me's,  .Many  approai  lies  can  be 
take'll  to  pre'proce'ss  the'  data,  e'ach  giving  diffe'rent  re'sults.  Continuing  research 
should  include  changing  (he  [taramelers  u.se'd  in  the'  pre'|)re>e’e'ssing  in  hopes  ol  dis- 
cewering  t  he'  Ix'st  pre'proe-essing  se-he-me.  One  e'.xample'  weeuld  be'  te)  analyze  t  he  la-sult  s 
of  log  compre'ssing  the'  data  siiiea'  limite'd  time'  diel  not  allow  a  chance  to  iiiNmke  this 
o|)lion.  In  addition,  the'  |)re'])roe-('.ssing  prergram  should  be'  moelilie'd  to  add  a  noise- 
tliie'shold  ahe-ad  eif  the  nerrmalizatiexi  routine'. 

Proce'ssing  was  rnaele-  more  elilliemlt  with  the-  Ef'd'  program  ein  the  mainlrami- 
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aiul  the  necessary  transh'r  of  fih's.  The  files  had  to  l)e  processt'd  one  at  a  time  hy 
tfi(’  ns<'r.  Idle  {)rograni  could  be  modified  to  run  mnltiple  programs  \sithoiil  iim'I’ 
intervention,  but  the  files  would  still  have  to  be  transferred  to  the  mainframe  and 
back.  Ffy  complet  ing  the  work  necessary  toenalrh'  the  fft.c  program  to  run  on  the  PC', 
all  data  gathering  and  processing  can  be  performed  on  tlu'  same  comjniti'r.  1  he  [fl.c 
program  can  he  combined  with  the  prcproc.c,  tcinj>latf .r,  and  distaprr.r  programs 
to  process  multiple  files  without  user  intervention.  Add  to  this  the  possibility  of 
automatic  data  gatliering  with  the  computer  controlled  rotating  stage  triggering  tin' 
S|nricon  frame  grabber  to  grab  a  franu'  at  each  incri'iiu'iit  of  rotatioi..  and  further 
resi'arch  can  he  done  more  efncii-ntly.  Onc<'  the  tenpilate  files  are  created  the  user 
can  move  the  camera  to  an  arbitrary  angh',  call  a  program  on  the  PC  winch  would 
(h)  all  the  necessary  processing  and  comparisons  to  the  tcmplati'  file,  and  hud  the 
.\().\  of  the  ti'st  \ecfor  in  the  output  file. 

In  practici'  this  detector  will  lx*  snbj<'ct  to  various  intensities  of  laser  radiation, 
huriher  ri'search  should  involve  using  the  variable  attenuator  to  changt'  the  intensity 
of  incoming  las('r  radiation  and  modifying  th('  pr<'])roc('ssing  program  to  adjust  for 
this  change.  1  he  normalization  program  currently  in  \ise  selects  the  pc’ak  value 
in  each  data  fih'  and  calcnlat('s  a  scaling  factor  by  which  to  multiply  each  data 
jioint ( pix('!  N'alue)  such  that  the  new  peak  value  is  2o^>.  d’his  e.xpands  the  range  of 
\'alues  in  each  file  to  the  ma.ximum  but  does  not  adjust  the  ovc'rall  int('nsity  of  each 
file  to  a  common  reh’rence.  1  hns  tin'  normalization  routine'  shoidd  be  a])p('nd('d  to 
allow  cempensation  for  variable'  lase'r  inle'iisily  by  adjusting  all  files  to  a  common 
re'fe'i'enee'  belore'  pe'rfe)rming  the'  ])e'ak  normalization. 

rile-  oidy  sejtirce'  nse'd  in  this  t  he'sis  was  a  hOm\V  lle'-.\e  lase-r  with  an  eaitient 
wa\e'length  of  tid'd. Sunn  Sine-e-  .\d-''>’.'\(!  lasers  with  a  wave-length  of  l.lKi//;;/  aie' 
pre'vale'iit  in  military  a|)piicat ieeiis  fen-  targe-t  illumination,  the-  re-latieni  of  wavelength 
to  the'  lilx'r  diame-te-r  aiiel  libei-  h-ngth  shoedd  lx-  e'xamine'el. 

.\ne)t  he'i'  re'ce)nime'ndat  ion  invedve's  |)roce'ssing.  ( 'urre-nt  ly  t  he-  |)r(x-essing  scheme' 


iii\olv('s  creating  a  test  feature  vector  from  the  low  freriueiu  y  l\)uri<'r  components 
and  comparing  this  test  vector  to  previously  created  trunplate  feature  vectors  us¬ 
ing  Minkowski  distance  for  classification.  This  is  one  of  the  h-ast  eirecti\(‘  wa\s  to 
classify  data.  In  order  to  utilize'  information  from  all  frecpu'ucies,  a  ccurelation  he- 
twee'ii  t('st  and  template  vectors  can  he  performed  in  tlu'  frecpiency  domain.  1  hi' 
can  l)e  accomplished  by  modifying  existing  software  such  as  Cie-orge  LeTourneau  s 
I'k  r  ])rogram.  using  I.MSL  library  programs,  or  developing  a  working  vc'ision  ol  the 
rc'cc'utly  created  FFT  program  written  in  ‘Turbo  ("  for  the  PC.  ib  ])('rform  the  cor- 
relati(U),  both  the  real  ajid  imaginary  array  must  be'  accc’ssible'.  Fake-  the'  conpigate'  of 
the  imaginar\‘  component  ol  the-  template  array  and  the-n  |)e'rform  a  point  for  point 
comple'x  midtiplicat ion  betwee'ii  the  te'inplate  and  te'sl  arrays.  1  he'  multiplication 
which  \-ie'lds  the  highest  correlation  pe'ak  classifie's  the'  te'st  \’ector.  .Anothe'r  proci'ss- 
ing  sclie'me  involves  using  a  Neural  .Net  to  classify  the  test  \’ector.  Ib’cause  .Neural 
.Ne'ts  are  inhe'reiitly  slow,  care  must  be  taken  in  selecting  the'  inputs.  One'  ])ossible 
set  of  inputs  would  lie  tin'  4!)  low  free|uency  Fe>urie'r  e'ompone'ut s  currently  used  in 
the  Minkowski  distance  calculat ieui. 

.Many  data  reduction  and  proce'ssing  sche-me's  e'xist  and  the'  limited  re'se-arcli 
conducted  in  this  the'sis  will  hoi)efully  |)roviele'  ine-entive'  for  furthe'r  inve'st  igat  ion  in 
tills  are'a. 
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Appendix  A.  Spiricon  2250  Laser  Beam  Diagnostics  System 


This  appendix  was  written  to  assist  the  user  in  repeating  the  data  collection 
pc'i  formed  in  this  thesis.  Any  information  not  specifically  covered  here  can  he  found 
within  one  of  the  two  Spiricon  manuals  retained  by  7\FIT. 

Ihe  Spiricon  2250  Laser  Beam  Diagnostics  system  contains  a  processing  card 
which  uses  two  slots  in  an  IBM  compatible  personal  computer(PC).  Vf-rsion  5.0  Laser 
B('am  Diagnostic  software,  a  CID'l'ECB  512  x  512  Charge  Injection  Device(CID) 
camera,  a  real-time  display  monitor,  and  the  necessary  connecting  cabh's.  I  he  2250 
system  is  adverti.sed  to  uork  within  an  JB.M  .\d'.  IB.M  .580,  or  IB.M  cmnpatible.  .After 
much  experimentation  with  various  PC's,  it  was  discoverc'd  that  the  2250  aiipari'iitly 
only  works  on  an  IBM  AT.  'Lhis  may  be  due  to  an  addressing  problem  within  the 
/.-218‘.s  and  IBM  386  in  which  the  2250  system  failed  to  operatic 

Tlie  version  5.0  .software  is  contained  on  two  360Kbyte  iio])py  disks,  d'lu'  files 
were  copied  onto  the  hard  drivi'  into  the  directory  c:\spir.  Subdirectories  weri-  also 
created  to  handle  the  various  stages  of  the  data.  The  subdirectory  H .AW'D.A'F.A  con¬ 
tained  the  original  frames  of  data  received  from  the  CID  camera.  These  files  contain 
263.052  nnsigne<l  bytc's.  d'he  structure  of  these  fih's  will  be  addri'ssed  shortly.  The 
subdirectories  PROCD.Ad'.A  and  Fl'D.ATA  contained  piejirocc'ssed  data  and  Fourier 
transformfl'  1  )  data  ri'spect ively.  'I  hese  files  contain  65.536  nnsigiu'd  bytes.  Finallv 
the  subdiri'ct ories  I  Ic.MPL.Al  F,  and  I’l'iS'I  contain  t  lu'  h'atnre  vectors  composed  of 
th<'  low  frecjuency  fonrier  components  representing  th<'  sha|)e  of  tin'  original  images. 

1  liese  files  vary  in  lengt  li  depi'iiding  on  th('  mimbc'r  of  vi'ctors  contained  in  t  he  fill's, 
f  ill's  cmitaining  inferrmat  ion  relating  to  the  results  of  data  ])rocessing  are  also  stored 
in  t  hesi'  last  two  subdirectories.  1  he  fih's  contained  in  all  subdirectories  will  be 
examined  in  greati'r  detail  in  the  Ujrcoining  section  on  proci'ssing. 

lo  start  S[)iricon,  type  CQ3D  from  within  tin-  c:\spir  directory.  The  main 
mi'iiu  will  appear  on  the  comjiuter  monitor  within  a  f<’w  si'i-onds.  Several  sub  ineiins 


arc  available  from  the  main  menii.  Most  of  these  are  not  of  immediate  fonrern  and 
will  not  be  addressed  here.  Type  F7  to  open  the  setup  menu.  Id  allow  Spirieon  to 
access  fdes  in  the  subdirectories,  type  [ALT]F2  and  the  subdirectory  name  followed 
by  a  reinrn.  To  check  operation  of  the  camera,  type  F2  to  access  the  b('amlink  menu. 
Typing  F2  again  toggles  the  camera  on  and  off.  ddie  camera  outjjut  is  displave<i  on 
the  real-time  monitor.  Typing  FIO  automatically  turns  the  camera  off  and  grabs 
the  frame,  as  well  as  returning  to  the  setup  menu.  Typing  FIO  again  returns  to  the 
main  menu.  Online  help  is  available  from  within  any  menu  by  typing  FO. 

Now  tliat  the  camera  has  been  found  to  be  working  correclly  and  the  piupcj 
subdirectory  has  been  accessed,  the  functions  found  within  the  main  and  auxiliary 
menus  suffice  to  carry  out  the  remaining  operations  needed  for  this  thesis.  Fl  grabs 
a  iK'w  frame  of  data  and  displays  it  on  the  real-time  monitor.  S'everal  si’conds  later 
a  :?-D  plot  is  displayed  on  the  computer  monitor.  This  3-D  plot  can  be  modified  in  a 
varii'ty  of  ways  from  the  auxiliary  menu.  The  3-D  plot  is  a  128  x  128  representation 
of  the  512  X  512  camera  array.  The  slice  regulates  the  resolution  of  this  128  x  128 
ariay.  A  slice  of  }  provides  the  highest  re.solution  when  looking  at  the  ent ire(  i/'/m/r ) 
512  X  512  array.  A  slice  of  4  tells  Spirieon  to  re])resent  the  average  of  4  pixids  in  th('  x- 
direction  and  4  ])ixels  in  the  y-directioii  as  one  pixel  in  the  3-D  jilot.  Slice  values  can 
increase  by  powers  of  2  to  the  upp<‘r  limit  of  G4.  l^arge  slice  values  decrease  resoint  ion 
as  well  as  processing  time.  Pre.ssing  iVl  increa.ses  the  slice  value  while  jiressing  L  will 
decrease  the  slice.  There  are  two  cur.sors  on  the  plot  rejn’esented  by  crossc's  at  two 
diagonal  corners.  In  ord('r  to  achieve*  maximum  resolution(.s//cf  =  1),  the  user  must 
zoom  in  on  an  aica  no  larger  than  128  x  128.  Idiis  is  accomplished  In'  nun-ing  the* 
cursors  into  the  jilot  using  the  cursor  keys  on  the  keyboarel.  4  Ik*  llOMF  key  toggles 
b('t  wc'en  cursors.  4'he  cooialinates  of  <’ach  cursor  are*  displaye'd  at  the  bottom  of  the 
auxiliai  v  menu.  .Note  that  the*  cui’sors  move*  in  increiiK'iits  bast'd  on  tlu'  slice'  \'alue'. 
Once  t  he  cursors  are*  pexsit  ione'd.  press  W.  fliis  tells  Spiriceen  to  plot  only  the*  imii 
within  the'  enrseus.  N'eew  |)i'e'ss  L  as  ne'ce'ssary  tee  reeluce'  slice'  to  1.  Pressing  return 


will  cause  Spiricon  to  reprocess  the  image  according  to  the  new  parameters,  riiere 
now  exists  a  one  to  one  correspondence  between  the  pixels  on  the  3-D  plot  and  the 
pixels  on  the  camera  within  the  representative  area.  To  read  the  intensity  value  of 
a  particular  pixel,  just  move  the  cursor  to  that  pixel  and  read  the  value  of  Z  at  the 
bottom  of  the  auxiliary  menu.  The  eight  bit  intensity  range  of  0  to  255  is  represented 
in  Spiricon  by  values  from  0  to  2040.  Another  variable  is  the  magnification  of  the 
display.  If  the  display  is  too  tall  to  fit  on  the  screen,  the  user  can  reduce  the  height 
by  pressing  -  on  the  keyboard.  Typing  -f  as  necessary  will  increase  magnification. 
\5dues  of  magnification  vary  from  1  /2  to  8.  The  display  can  also  be  tilted  in  5  degree 
increments  from  -90  to  +90  degrees  and  rotated  in  15  degree  increments  from  -180 
to  +180  degrees  by  pressing  PgUp/PgDn  and  Ins/Del  respectively.  There  are  ten 
different  formats  for  the  3-D  plot.  Some  of  these  add  contour  lines  representing 
lc\el  curves  to  the  plot.  These  formats  are  changed  each  time  the  X  key  is  pressed. 
Once  the  desired  parameters  are  set,  they  can  be  saved  by  pressing  FT  to  access  the 
setup  menu.  Now  press  [AItjFl  to  save  the  configuration.  Press  FIO  to  return  to  the 
main  menu.  The  parameters  are  saved  in  the  file  SPIR.FIG  which  is  90C  bytes  long. 
Whenever  Spiricon  is  initially  loaded,  it  looks  for  the  file  SPIR.FIG  and  sets  the 
parameters  accordingly.  There  are  two  ways  to  send  the  plot  to  a  ])rint('r.  'fypiiig  P 
peiforms  a  screen  dump,  'flius  menus  are  printed  as  will  as  the  jilot.  To  print  the 
3-1)  plot  and  various  other  information  without  the  menus,  prc'ss  1''2.  To  save  the 
plot  to  disk  with  the  curixuit  parameters,  jiress  Ft,  enter  a  filename  of  your  choice, 
and  press  return.  The  user  is  now  prompted  for  a  record  number.  Spiricon  allows  as 
many  as  32,000  records  to  be  stored  under  each,  filenanu'.  'Fliis  is  ovi-rkill  sinci'  each 
record  is  263,053  l)3tes  long,  dype  1  and  hit  return.  If  record  1  alrcxuly  exists,  the 
user  can  overwrite  it  or  choose'  anotln'r  record  or  fih'uame.  The'  jilot  is  stored  with 
tlu'  parameters  current h'  in  use.  Once  the  file  is  stoix'd  according  to  tlu'  data  |)ath 
sclectc'd  ('arlier.  the'  main  menu  once'  again  appears  on  the  scre'e'ii.  lb  loael  another 
file-,  pre'ss  I'M.  1  he  use'r  is  gi\'en  a  list  of  files  in  t  lie*  elire'e  tori’  chosen  earlier.  Pre'ss 


It) 


the  up  or  down  cursor  keys  to  select  a  file  and  hit  return.  If  the  required  file  is  in 
another  directory  the  user  can  change  directories  by  pressing  D  and  typing  the  new 
directory  path.  This  is  an  alternate  way  of  changing  directory  path.  Files  can  also 
be  deleted  by  selecting  a  file  with  the  cursor  keys  and  pressing  the  DEL  key.  Once 
a  file  is  selected  for  loading,  type  return.  The  chosen  file  is  loaded,  displayed  on  the 
real-time  monitor,  sliced,  and  then  plotted  on  the  computer  monitor  according  to  the 
current  parameters.  Sometimes  it  is  desired  to  display  the  file  with  the  parameters 
set  as  they  were  when  the  file  was  stored.  This  is  achieved  by  typing  F7  from  the 
main  menu  to  access  the  setup  menu  before  trying  to  load  the  file.  Type  F7  again 
to  access  the  miscellaneous  menu.  Now  type  F6.  This  tells  Spiricon  to  load  files  and 
set  parameters  according  to  the  information  stored  in  the  files.  Type  FlO  twice  to 
return  to  the  main  menu.  To  quit  this  feature,  repeat  the  procedure. 

Two  anomalies  were  noted  while  using  the  Spiricon  system.  It  was  noted  that 
Spiricon  gave  an  error  message  when  trying  to  load  a  file  from  floppy  disk  although 
this  is  an  allowed  function.  Thus  any  data  files  from  another  system  converted  to 
Spiricon  format  had  to  be  copied  to  the  hard  drive  before  being  loaded.  The  second 
anomaly  concerns  a  user  created  file  called  BADPIX.  This  file  contains  coordinates 
of  bad  pixels  located  on  the  camera  sensor.  These  bad  pixels  are  very  sensitive 
and  register  at  the  maximum  value  when  little  or  no  light  is  present.  This  resulted 
in  the  warning  message  “Usable  range  of  camera  exceeded”.  This  message  occurs 
whenever  intensity  values  exceed  approximately  80%  of  peak.  Spiricon  is  supposed 
to  examine  this  file  and  ignore  these  bad  pixels  in  t  data  display.  This  file  was 
never  used  successfully  although  calls  were  made  to  Spiricon  Incorporated  and  an 
ui)dated  version  of  software  was  finally  sent. 

As  mentioned  earlier,  Spiricon  stores  the  data  received  from  the  camera  in  files 
2Gd.().'32  bytes  in  length.  The  first  2  bytes  contain  the  HEX  digits  .‘59  and  30.  Tlu'.se  are 
used  by  Spiricon  to  find  the  beginning  of  a  data  file.  The  next  262,144  bytes  are  the 
data  from  the  camera.  The  last  906  bytes  contain  the  display  parameters  which  were 


in  cfTect  when  the  file  was  stored.  Notice  this  is  equivalent  to  the  SlNR.FKl  file.  1  lu' 
data  is  written  to  the  real-time  monitor  starting  at  the  top  and  working  downward, 
however,  only  every  other  row  of  data  supplied  by  the  camera  is  displayed  on  the 
real-time  monitor.  All  columns  are  displayed.  Although  the  first  row  of  data(r>r2 
bytes)  in  the  file  represents  the  top  row  of  the  real-time  monitor,  it  is  listed  as  row 
512  on  the  3-D  plot  since  the  origin  is  in  the  lower  left  hand  corner  of  the  image.  ,411 
bytes  in  the  last  11  rows  in  the  data  file,  corresponding  to  the  first  1 1  rows  on  the 
3-D  plot,  are  always  zero. 

Spiricon  has  many  more  features  such  as  gain  adjustment,  curvefit,  and  capture 
length  that  are  not  addressed  here.  Spiricon  can  al.so  be  used  for  continuous  or  pulsed 


Appendix  B.  Program  LOG3X3.C 


The  program  logSxS.c  is  used  on  full  size  Spiricon  files.  This  program  cloes 
not  have  any  reduction  routines,  therefore,  the  output  remains  a  full  size  Spiricon 
file.  The  program  is  written  in  loops.  There  is  one  large  outside  loo])  containing  five 
subroutines.  These  subroutines  are  median  filter,  normalization,  thresholding,  peak 
value,  and  random  access.  The  last  two  subroutines  are  contained  in  individual  loops 
to  allow  repetition  without  running  through  the  entire  program  each  time.  The  flow 
chart  in  figure  16  demonstrates  this  process. 

Upon  entering  the  program  the  user  is  asked  if  the  data  is  to  be  median  filtered. 
If  the  user  answers  yes  then  he  is  prompted  for  a  data  path  and  filename  without 
the  three  letter  extension.  The  filename  in  this  case  cannot  have  an  ('xtimsion.  d'hc 
program  then  reads  FILENAME  and  the  data  is  median  filteri'd  whih'  the  peak 
intensity  and  the  sum  of  the  intensities  is  also  calculated.  I’he  median  filtered  data 
is  automatically  placed  in  FILENA.ME..MD3  while  the  intensity  data  is  displayed  on 
th('  monitor.  The  program  now  asks  the  user  if  he  wishes  to  normalize'  to  the  ])eak 
value.  If  yes,  then  the  program  automatically  reads  FILENAME..M  1)3  and  opi'iis 
output  file  FILENAME. NM3.  The  user  can  choose  the  value  he  wishes  to  s('t  the 
[)eak  value  (i.e.  normalize)  to,  log  compress  the  data  if  so  desired,  and  also  display 
a  histogram  of  the  data  on  the  monitor.  For  example,  if  the  jx-ak  value  of  the  data 
is  119  and  the  user  wishes  to  normalize  to  2.')5.  then  a  scaling  value  is  calculated 
which  converts  119  to  2.5.5.  The  scaling  value  is  then  multiplied  by  ('acli  data  point 
in  the  array.  Log  compression  calculates  the  logarithm  of  each  data  point  in  base' 
10.  calculates  a  different  scaling  factor,  and  then  normalizes,  d  his  oi)lion  efh’ct i vc'ly 
raises  the  low('r  intensities  such  that  the  data  is  not  s|)r('ad  over  a  large  range  of 
valiu's.  d  Ik'  normalization  loutiiK'  obtains  tlu'  peak  valiK'  from  the  iiu'dian  lilP'r 
routine'.  Once  tlu'  data  is  normalized,  the  data  is  se'gnu'iite'd  into  Ki  gioiips  of  Hi 


Figiii'c  IG.  Flowcliart.  For  I’rograin  lofi-lr.i.c 


v;ilu(“s.  This  information  is  disi)la.yed  on  the  monitor  upon  entering  the  thrc'sholding 
routine  and  is  used  to  help  choose  a  threshold  value. 

If  the  data  was  iiot  median  filtered,  the  program,  upon  entering  the  normaliza¬ 
tion  routine,  will  ask  for  an  input  filename  without  the  three  letter  extension  .  1  his 
time,  however,  the  program  appends  the  .MD3  extension  to  the  filename  which  was 
entered,  such  that  the  normalization  routine  still  reads  FILENAME. MD3.  It  then 
calculates  the  peak  value  since  it  was  not  calculated  in  the  median  filter  routine, 
d  h('  normalization  routine  then  runs  as  previously  described. 

Following  the  normalization  routine  is  the  thresholding  routine.  It  automat¬ 
ically  reads  FIIT1NAME.NM3  and  opens  output  file  FILENAME. TII3.  If  the  his¬ 
togram  option  was  chosen  in  the  normalization  routine,  it  is  displayed  on  tin'  screen 
at  this  time.  An  example  of  this  display  is  shown  in  figure  11  in  chapter  il.  The 
program  then  displays  the  peak  value  of  the  data  file  and  asks  the  user  to  enter  a 
threshold  vahu'.  All  values  below  threshold  are  set  to  zero.  The  user  can  also  decide 
to  binarize  the  data.  If  this  option  is  chosen,  then  all  the  values  above  the  threshold 
ai'e  s('t  to  205.  Otherwise  the  values  above  the  threshold  remain  as  they  are. 

If  the  normalization  routine  was  not  entered,  then  the  program  asks  fi)r  the 
directory  path  and  filename  without  extension.  It  automatically  appends  the  exten- 
sicjii  .NM.3,  such  that  it  still  r<‘ads  F1I>I'1N AME.NM3.  'hhe  i^rogram  then  asks  if  a 
histogram  of  the  flata  is  desired.  'The  program  then  calculates  the  maximum  inten¬ 
sity  and  segments  the  data.  It  tlnni  <lisplays  a  histogram  il  that  option  was  chosen, 
following  this  the  ])rogram  |)erforms  as  ch'seribed  in  the  previous  i)aragraph. 

rh('  next  routine  is  dedicated  to  finding  the  coordinates  of  the  peak  int<'nsity, 
valiK'  (jf  the  jx'ak  intensity,  and  the  sum  of  the  intensiti('s.  d'his  routine'  expects  a 
data  |)ath  and  filename  with  extension.  I’lins.  it  does  not  r('f|uire  a  spe'cilic  ('xtension 
since  the  user  enters  this  information.  Since  this  [trogram  lies  within  a  while  loop,  it 
can  b('  re  entered  immediately  upon  completion  and  gi\'en  a  new  filename  tcj  access. 


Following  this  routine  is  a  random  access  routine.  This  allows  the  user  to 
display  a  particular  area  of  data  within  a  fde  on  the  monitor.  The  user  once  again 
chooses  a  data  path  and  filename  with  extension.  The  user  also  choo.sc's  the  beginning 
row  and  column  and  the  ending  row  and  column.  There  are  512  columns  i)er  row  in 
the  data  file,  however,  the  display  on  the  screen  is  only  16  columns  wide.  Thus,  one 
row  from  the  data  file  requires  16  rows  on  the  screen.  This  routine  also  lies  within 
a  while  loop  and  can  be  accessed  as  many  times  as  necessary.  Upon  completion  of 
this  routine,  the  program  queries  the  user  for  changes.  If  the  user  wishes  to  make 
changes,  the  program  returns  to  the  beginning  and  starts  over.  If  changes  are  not 
re(|uired,  the  program  is  terminated. 

It  should  be  noted  that  the  routines  display  the  data  pointer  location  on  the 
screen  while  the  program  is  running.  In  some  ca.se.s  the  pointer  ma\'  indicate  that 
the  program  did  not  access  all  262, 141  bytes  of  data.  This  is  because  S])iricon  ])laces 
zeros  in  the  last  5632  bytes  of  the  file.  Thus,  some  of  the  calculations,  such  as  the 
intc'iisity  summation,  do  not  need  to  operate  beyond  the  first  501  rows  of  the  file. 

In  summary,  the  median  filter  routine  reads  a  filename  with  no  extension  while 
tin'  normalization  routine  reads  a  filename  with  the  .MD3  extension  and  the  thresh¬ 
olding  routine  reads  a  filename  with  the  .TH3  extension,  even  though  the  user  only 
(Miters  a  data  path  and  filename  with  no  extension  (i.e.  c:\spir\rawdata\myfile). 
riie  program  was  designed  in  this  manner  so  that  a  user  could  go  smoothly  through 
each  of  the  first  three  routines,  exit  the  program,  and  comparf'  the  results  of  each 
routine  since  each  routine  sends  iUs  output  data  to  a  different  fih'. 


/♦Program  L0G3X3.C*/ 


/♦Written  by  Capt  John  W.  Welker^/ 


/♦This  program  performs  the  pre-processing  on  a  Spiricon  image  file*/ 
/♦This  pre-processing  consists  of  a  median  filter  operation  which*/ 
/♦removes  spikes  in  the  data.  The  data  is  then  energy  normalized  on*/ 
/♦a  linear  scale.  The  data  may  be  log  compressed  to  increase  the*/ 
/♦lower  intensities.  A  thresholding  option  is  also  available  which*/ 
/♦will  set  all  bits  below  the  threshold  to  zero.  Bits  above*/ 
/♦threshold  can  be  set  to  255  or  left  unchanged.  An  option  to  find*/ 
/♦the  coordinates  of  the  peak  value  and  to  display  a  random  area  of*/ 
/♦data  is  also  available.  Recursion  means  that  the  median  value  is*/ 
/♦placed  into  the  middle  of  the  3x3  matrix  before  it  calculates*/ 
/♦the  next  median  value.*/ 

/♦COMPILER  DIRECTIVES*/ 

#include  <stdio.h> 

#include  <string.h> 

#include  <math.h> 

/♦DEFINE  MACROS*/ 

#define  inrow  104 
#define  incol  512 
#define  slice  5 

/♦DECLARE  POINTERS  TO  FILE  STRUCTURE*/ 

FILE  ♦ptrl,*ptr2; 

/♦START  MAIN  PROGRAM*/ 
void  mainO 
{ 

/♦DECLARE  LOCAL  VARIABLES*/ 

unsigned  char  filename [256] ,median_f ile [256] .array [inrow] [incol] ; 
unsigned  char  temp[9] ,random_file_access[256]  ; 
unsigned  char  norm_f ile[256] ,thresh_f ile[256] ,max_f ile[256] ; 
unsigned  record_num , row , col ,max=0 .count  1 . count 2 . counts , c . r , a=0 ; 


unsigned  length=0  ,peak=0  ,norni_maLX  ,threshold=0 ; 
unsigned  begin_row,begin_col ,end_row,end_col ,random_length=0 ; 
undigned  randoin_row=0 ,  random_col=0 ,  peakinax=0  ,peakrow=0  ,  peakcol=0 ; 
unsigned  long  sum=0  ,  of  f  set=0  ,num_byte=0  ,  slice_suin ; 
unsigned  long  norm_sum ,ptr_off set=0 ,peakbyte=0 ; 

unsigned  long  segment  1 , segment 2, segment 3, segment4, segments , segments; 
unsigned  long  segment? , segments ,segment9 , segmentlO , segment  11 ; 
unsigned  long  segmentl2,segmenti3,segmentl4,segmentl5,segmentl6; 
char  prograLm='y'  .median  .normal  .thresh,  random  .binary  .maxbyte  ; 
char  log_scale,histo,headerC]  =  {'\x39’,  '\x30’}; 
float  norm_scale=0 . 0 , thresh_scale=0 . 0 ; 

length=inrow*incol ; 

while (program== ' y ' ) 

{ 

printf ("\n\nDo  you  wish  to  smooth  data  with  a  median  filter, 
y  or  n?  \n")  ; 
scanfC"  y.c"  .ftmedian) ; 
if (median=='y' ) 

{ 

sum=0 ; 
max=0 ; 
num_byte=0; 
peakmax=0 ; 
peakrow=0 ; 
peakcol=0 ; 
peakbyte=0 ; 

/♦ENTER  FILES  AND  VARIABLE  VALUES  FROM  KEYBOARD*/ 

printf ("\nEnter  directory  path  and  file  name  of  input  file  without 

extension;  \n\n") ; 

scanf  ("'/,s"  .filename)  ; 

printf ("\nEnter  record  number  \n"); 

scanf  ("'/.u"  ,&record_num)  ; 

sprintf  (median_f  ile ,  "'/,s  .md3"  .filename)  ; 

/♦OPEN  INPUT  AND  OUTPUT  FILES*/ 

printf  ("VnOpening  input  file  */,s  \n"  ,  f  ilename)  ; 


if (  (ptrl=f open(f ilename, "rb"))==NULL) 

{ 

printf ("Can't  open  input  file  %s  \n" , filename) ; 
exit ()  ; 

} 

printf  ("Opening  output  file  ‘/,s  \n\n"  ,median_f  ile)  ; 
if  (  (ptr2=f open (mediaii_f ile ,  "wb")  )==NULL) 

{ 

printf  ("Can’t  open  output  file  */,s  \n"  ,median_f  ile)  ; 
exit  0  ; 

> 

/♦MOVE  POINTER  TO  DESIRED  RECORD  IN  FILE*/ 
of f set=(record_num-l)*263052+2 ; 
f seek(ptrl , offset , 0) ; 

f write (header , 1 ,2 ,ptr2) ; 


/♦BEGIN  MEDIAN  FILTER  PROCESSING^/ 

for (count  1=0 ; count 1< si i ce ; count 1++) 

{ 

if (countl==0)  /♦FIRST  PASS  ONLY*/ 

printf  ("reading  slice  #’/,u  \n"  ,  count  1  +  1)  ; 

f read (array , 1 , length ,ptrl) ;  /♦FIRST  READ  INTO  ARRAY*/ 

printf  ("ptrl  =  */,lu  \n\n"  ,f  tell  (ptrl)  )  ; 

r=0; 

f or (count2=0 ; count2<inrow-2 ; count2++)  /♦FIRST  MEDIAN*/ 

/♦FILTER  LOOP*/ 

{ 

c=0 ; 

for(count3=0 ; count3<incol-2 ;count3++) 

{ 

a=0 ; 

f or (row=r ; row<r+3 ; row++)  /♦SELECT  VALUES  FOR*/ 

{  /*3  X  3  ARRAY*/ 

for(col=c;col<c+3;col++)  /♦AND  PREPARE  FOR  SORT*/ 

{ 

temp [a] =array [row] [col] ; 
a+": 


YS 


} 


} 

sort (temp, a) ;  /*fORT  3X3  ARRAY*/ 

C++; 

arrayCr+l] [c]=temp[a/2] ;  /*PLACE  MEDIAN  VALUE  IN*/ 
}  /*CENTER  OF  3  X  3  ARRAY*/ 


r++; 

} 

f or(col=0;col<incol ;col++) 

array [0] [col] =array [l] [col] ; 
f or (row=0  ;row<inrow-2 ;row++) 

array [row] [0] =array [row] [l] ; 

array [row] [incol-l] =array [row] [incol-2] ; 

} 


/*SET  FIRST  ROW  EQUAL*/ 
/*T0  SECOND  ROW*/ 

/*RESET  FIRST  AND  LAST*/ 
/*C0LUMN*/ 


printf ("Calculating  first  slice  intensity  and  maximum  \n\n") ; 
slice_sum=0 ; 

for(row=0;row<inrow-2;row++)  /*CALCULATE  MAXIMUM  VALUE*/ 

{  /*AND  INTENSITY*/ 

f or (col=0 ; col<incol ;col++) 

{ 

sum+=array [row] [col] ; 
slice_sum+=array [row] [col] ; 
num_byte++ ; 

max  =  (max  >  array [row]  [col] )  ?  max  : array [row] [col] ; 
if(peakmax  !=  max) 

{ 


peakiaax=max ; 
peakbyte=num_byte ; 

} 

} 

} 

printf  ("The  sum  of  this  slice  is  V.lu  \n\n"  ,  slice_sum)  ; 
printf  ("Writing  slice  ft'/.u  \n"  ,countl+l)  ; 

fwrite(array,l,length-2*incol,ptr2) ;/*WRITE  ALL  BUT  LAST  TWO*/ 
printf  ("ptr2  =  ‘/.lu  \n\n"  , f  tell  (ptr2)  )  ; /*R0WS  TO  OUTPUT  FILE*/ 

} 

else  /*BEGIN  ALL  SUCCESSIVE*/ 

{  /*PASSES  THROUGH  FILTER*/ 

f or (row=0 ; row<2 ; row++) 

{ 

for(col=0;col<incol;col++)  /*LAST  TWO  ROWS  ARE  NOW*/ 


{  /*FIRST  TWO  ROWS*/ 

array [row] [col] =array [inrow-2+row] [col] ; 

> 

} 

printf  ("reading  slice  #*/,u  \n"  .count  1+1)  ; 
f read(&array [2] [0] , 1 ,length-2*incol ,ptrl) ; 

/*READ  IN  SECOND  SLICE  OF*/ 

/*DATA  STARTING  IN  THIRD*/ 

/*R0W  OF  ARRAY*/ 

printf  ("ptrl  =  */,lu  \n\n"  ,  f  tell  (ptrl)  )  ; 
r=0 ; 

for (count 2=0 ; count2< inrow-2 ; count 2++) 

{  /*SUCCESSIVE  MEDIAN  FILTER*/ 

c=0;  /*L00PS*/ 

for (count 3=0 ; count 3<incol-2 ; count3++) 

{ 

a=0 ; 

f or (row=r ; row<r+3 ; row++)  /*SELECT  VALUES  FOR  3  X  3*/ 

{  /*ARRAY  AND  PREPARE  FOR*/ 

for(col=c; col<c+3 ; col++)  /*S0RTING*/ 

{ 

temp [a] =array [row] [col] ; 
a++; 

> 

} 

sort (temp , a) ;  /*S0RT  3X3  ARRAY*/ 

C+  +  ; 

array [r+1] [c]=temp[a/2] ;  /*PLACE  MEDIAN  VALUE  INTO*/ 

}  /*CENTER  OF  3  X  3  ARRAY*/ 

r++; 

} 

f or (row=0 ; row<inrow-2 ; row++)  /*RESET  FIRST  AND  LAST*/ 

{  /*C0LUMNS*/ 

array [row] [0] =array [row] [l] ; 

array [row] [incol-l] =array [row] [incol-2] ; 

} 

if (countl<slice-l)  /*INTERMEDIATE  LOOPS*/ 

{ 

printf ("Calc  lating  intermediate  slice  intensity  and 
maximum  ',n\n")  ; 
si ice_sum=0 ; 


f or (row=0 ; row<inrow-2 ; row++) 

{ 

for(col=0;col<incol;col++)/*CALCULATE  MAXIMUM  VALUE*/ 

{  /*AND  INTENSITY*/ 

sum+=array [row] [col] ; 

slice_sum+=array [row] [col] ; 

num_byte++ ; 

max  = 

(max  >  array [row] [col] )  ^  max  :  array [row]  [col] ; 
if (peakmax  !=  nax) 

{ 

peakmax=maLX ; 
peakbyte=num_byte ; 

} 

} 

} 

printf("The  sum  of  this  slice  is  */.lu  \n\n"  ,  si ice_sum)  ; 
printf  ("Writing  slice  #'/.u  \n"  ,countl+l)  ; 
fwrite(array , 1 , length-2*incol ,ptr2) ; 

/*WRITE  TO  OUTPUT  FILE*/ 

/*DURING  INTERMEDIATE*/ 

/*PASSES  THROUGH  FILTER*/ 

printf  ("ptr2  =  */.lu  \n\n"  ,ftell(ptr2))  ; 

} 

else 

{ 

/*LAST  PASS  THROUGH  FILTER*/ 

/*RESET  FIRST  AND  LAST  COLUMN*/ 
array [inrow-2] [0] =array [inrow-2] [l] ; 
array [inrow-2] [incol- 1] =array [inrow-2] [incol-2]  ; 
f or (col=0 ; col<incol  ;col++)  /+RE5ET  LAST  ROW*/ 

array [inrow- 1] [col] =array [inrow-2] [col] ; 

printf ("Calculating  last  int  msity  and  maximum  \n\n")  ; 
slice_sum=0 ; 

f or(row=0 ; row<inrow ; row++)  /*CALCULATE  MAXIMUM  VALUE*/ 

{  /*AND  INTENSITY  IN  LAST*/ 

for(col=0; col< incol ; col++)  /*SLICE*/ 

{ 

sum+  =  array  [row]  [col] ; 
si ice_sum+=array  [row]  [col] ; 
num_byte++ ; 
max  = 


()1 


(max  >  array [row] [col] )  ?  max  :  array [row] [col] ; 
ifCpeakmax  !=  max) 

{ 

peakmax=max ; 
peakbyte=num_byte ; 

> 

> 

> 

printf("The  sum  of  this  slice  is  */,lu  \n\n"  ,  slice_sum)  ; 
printf  ("Writing  slice  #'/,u  \n"  , countl+1)  ; 
fwritev^array ,  1  .length, ptr2)  ;  /*WRITE  ALL  OF  LAST  SLICE*/ 
printf  ("ptr2  =  '/.lu  \n\n”  .ftell(ptr2))  ;  /*T0  OUTPUT  FILE*/ 

} 

} 

> 

if (length>906) 

{ 

f read (array , 1 ,906 ,ptrl)  ; 
f write (array , 1 ,906 ,ptr2)  ; 

> 

f close(ptrl) ;  /*CL0SE  FILES*/ 

fclose(ptr2) ; 

peakrow=(peakbyte/(long)incol) ; 
peakcol=(peakbyte  '/,  (long)  incol)-l ; 

printf  ("The  sum  of  the  first  V.lu  intensities  is  '/.lu  \n\n"  , 
num_byte , sum) ; 

printf  ("The  peak  intensity  of  ’/,u  after  filtering  occurs  at  \n",max 
printf  ("row  '/.u  and  column  '/.u  in  file  */,s  .  \n\n"  ,  peakrow  ,  peakcol , 
median_f ile) ; 

printf  ("This  corresponds  to  coordinates  (‘/,u,'/,u)  ".peakcol, 
512-peakrow)  ; 

printf ("on  Spiricon  plot.  \n"); 

} 


/*BEGINNING  NORMALIZATION  ROUTINE*/ 

printf  ("\n\nDo  you  wish  to  normalize,  y  or  n'^  \n")  ; 
scanf("  */,c"  ,&normal)  ; 
if (normal== ' y  0 
{ 

norm_sum=0 ; 


norm_max=0 ; 
num_byte=0 ; 
segment 1=0 ; 
segment2=0 ; 
segment3=0 ; 
segment4=0 ; 
segraent5=0 ; 
segment6=0 ; 
segment7=0 ; 
segment8=0 ; 
segraent9=0 ; 
segment  10=0 ; 
segment  11=0; 
segment  12=0 ; 
segment  13=0 ; 
segment  14=0 ; 
segment  15=0 ; 
segment  16=0 ; 
if (median ! = 'y ' ) 

{ 

max=0 ; 

printf ("\nEnter  directory  path  and  file  name  of  input  file 

without  extension:  \n\n") ; 

scanf  ("*/.s"  .filename)  ; 

sprintf  (median.f  ile ,  "'/,s  .md3"  .filename)  ; 

printf  ("\n0pening  input  file  */,s  \n\n"  ,median_f  ile)  ; 

if (  (pt rl=f open (median_f ile, "rb"))==NULL) 

{ 

printf  ("Can't  open  input  file  */,s  \n"  ,median_f  ile)  ; 
exitO  ; 

} 

f seek(ptrl ,2,0); 

printf ("\nCalculating  maximum  intensity  \n\n")  ; 

for (count  1=0 ; count l<s] ice ; count 1++) 

{ 

printf  ("reading  slice  #'/,u  \n"  ,countl  +  l)  ; 
f read (array , 1 ,length-2*incol ,ptrl) ; 
printf  ("ptrl='/,lu  \n\n"  ,ftell(ptrl))  ; 
f or(row=0 ; row<inrow-2 ; row++) 

{ 

for (col=0  ;col<incol ; col++) 

{ 


max=(max>array [row] [col] )  ?  max  : array  [row]  [col] ; 

> 

} 

} 

printf("The  maximum  intensity  in  file  */,s  is  */,u  \n" , 
median_f ile ,max) ; 
f close(ptrl) ; 

} 

sprintf  (norm_f  ile,  "’/,s  .nm3"  .filename)  ; 

printf ("\nWhat  number  do  you  wish  to  normalize  to  (2040  is 
peak)?  \n"): 
scanf  ("*/,u"  ,&peak) ; 

prinLi  ("\nHighest  intensity  will  be  */,u  \n",peak); 

printf("Do  you  wish  to  perform  a  log  compression  of  the  data, 

y  or  n?  \n") ; 

scanf  ("  y,c"  ,&log_scale)  ; 

if (log_scale== 'y ' ) 

norm_scale= (float )peak*255 . 0/2040 . 0/log( (double) max) ; 
else 

norm_scale=( (float) peak+0 . 5) *255 . 0/2040 .0/ (float) max ; 
printf  ("Scaling  factor  is  *'f  \n\n"  ,norm_scale)  ; 
printf ("Do  you  wish  to  see  a  histogram  of  the  data  ") ; 
printf ("for  thresholding  purposes,  y  or  n?  \n"); 
scanf  ("  y.c"  ,&histo)  ; 

printf ("\n0pening  input  file  y,s  for  normalization  \n", 
median_f ile) ; 

ptrl=f open (median_f ile, "rb") ; 

printf  ("Opening  output  file  '/.s  for  normalization\n\n"  ,norm_f  ile) 

ptr2=f open (norm_f ile, "wb") ; 

f seek(ptrl  ,2,0) ; 

fwrit e (header , 1 , 2,ptr2) ; 

for (count  1=0 ; count l<sl ice ; count 1++) 

{ 

printf  ("reading  slice  #y,u  \n"  ,countl+l)  ; 
f read(array , 1 , t  ength-2*incol ,ptrl) ; 
printf  ("ptrl=y,lu  \n\n"  ,f  tell  (ptrl) )  ; 

printf  ("Value  of  array[l][l]  before  scaling  is  '/.u  \n"  , 
array [1]  [1] )  ; 

f or (row=0 ; row<inrow-2 ; row++) 

{ 

for(col=0;  col<incol ; col++) 


(il 


{ 


if (log_scale== ’y ' ) 

{ 

arrayCrow] [col] =log((double)array [row] [col]+l) ; 

if(row==l  &&  col==l) 
printf ("Value  of  array [1]  [1]  after  log 
compression  is  '/,u  \n"  .array  [1]  [l]  )  ; 

} 

array [row] [col] *=norm_scale ; 
norm_sum+=array [row] [col] ; 
norm_max= 

(norm_max>array [row] [col] )  ?norm_max: array [row] [col] ; 
if (histo==' y ' ) 

{ 

if (array [row] [col]<16) 
segment 1++; 

else  if (array [row] [col] <32) 
segment2++ ; 

else  if (array [row] [col] <48) 
segment3++; 

else  if (array [row] [col] <64) 
segment4++; 

el39  if (array [row] [col] <80) 
segment5++; 

else  if (array [row] [col] <96) 
segment6++; 

else  if (array [row] [col]  <112) 
segment7++; 

else  if (array [row] [col]  <128) 
segment8++ ; 

else  if (array [row] [col] <144) 
segment9++; 

else  if (array [row] [col]  <160) 
segment 10++ ; 

else  if "ray [row] [col] <176) 
segment  11 ++; 

else  if (array [row] [col] <192) 
segment 12++ ; 

else  if (array [row] [col] <208) 
segment 13++ ; 

else  if (array [row] [col] <224) 
segmentl4++ ; 

else  if (array [row] [col] <240) 


segment 15++; 

else  if (array [row] [col] <=255) 
segment 16++; 
num_byte++ ; 

} 

} 

} 

printf  ("Value  of  array[i][l]  after  scaling  is  '/,u  \n\n" , 
array [1] [l] ) ; 

printf  ("Writing  slice  #y,u  \n"  .count  1+1)  ; 
fwrite(array , 1 ,length-2*incol ,ptr2) ; 
printf  ("ptr2=*/.lu  \n\n"  , f tell (ptr2) )  ; 

} 

printf("The  sum  of  the  intensities  after  normalization 
is  y,lu  \n"  ,norm_sum)  ; 

printf  ("The  normalized  maximum  is  y,u  \n"  ,norm_max)  ; 
f read(array , 1 , 1930 ,ptrl) ; 
f write(array , 1 , 1930 ,ptr2) ; 
f close(ptrl) ; 
fclose(ptr2) ; 

> 

/♦BEGINNING  THRESHOLDING  ROUTINE*/ 

printf ("\n\nDo  you  wish  to  set  a  threshold,  y  or  n?  \n") ; 
scanf("  y.c"  .ftthresh)  ; 
if (thresh=='y' ) 

{ 

if (normal ! = ' y ' ) 

{ 

max=0 ; 
num_byte=0 ; 
segment  1=0; 
segment2=0 ; 
segment3=0 ; 
segment4=0 ; 
segment5=0 ; 
segment6=0 ; 
segment7=0 ; 
segment8=0 ; 
segment9=0 ; 
segment  10=0 ; 


segment  11=0; 
segment  12=0 ; 
segment  13=0 ; 
segment  14=0 ; 
segment  15=0 ; 
segment  16=0 ; 

printf ("\nEnter  directory  path  and  file  najne  of  input  file 
without  extension:  \n\n") ; 
scanf  ("*/,s"  .filename)  ; 

?  orintf  (norm_f  ile ,  "Xs  .nm3"  ,f  ileneime)  ; 
printf  ("\n0pening  input  file  */.s  \n\n"  ,norm_f  ile)  ; 
if (  (ptrl=f open(norm_f ile , "rb") )==NULL) 

printf  ("Can  ’  t  open  input  file  */,s  \n"  ,norm_f  ile)  ; 
exit  ()  ; 

} 

f seekCptrl ,2,0); 

printf("\nDo  you  wish  to  see  a  histogram  of  the  data, 
y  or  n?  \n")  ; 

scanf  ("  */,c"  ,&histo)  ; 

printf ("\nCalculating  maximum  intensity  \n\n")  ; 
for (count  1=0 ; count 1< si i ce ; count 1++) 

{ 

printf  ("reading  slice  #'/,u  \n"  .count  1+1) ; 
fread(array , 1 ,length-2*incol,ptrl) ; 
printf  ("ptrl=y,lu  \n\n"  ,f  tell(ptrl)  )  ; 
f  or(row=0 ;  row<inrow-2  ;row+  •■) 

{ 

for(col=0 ;col<incol ; col++) 

{ 

max= (max>array [row] [col] )  ?  max  : array [row]  [col] ; 
if (histo=='y') 

if (array [row] [col]<16) 
segment 1++ ; 

else  if (array [row] [col] <32) 
segment2++ ; 

else  if (array [row] [col] <48) 
segment3++ ; 

else  if (array [row] [col] <64) 
segraent4++ ; 

else  if (array [row] [col] <80) 


(i? 


segment5++ ; 

else  if (array [row] [col] <96) 
segment6++ ; 

else  if (array[row] [col] <112) 
segment7++ ; 

else  if (array [row] [col] <128) 
seginent8++ ; 

else  if (array [row] [col] <144) 
segment9++ ; 

else  if (array [row] [ccl] <160) 
segment 10++; 

else  if (array [row] [col] <176) 
segment 11++; 

else  if (array [row] [col] <192) 
segment 12++; 

else  if (array [row] [col] <208) 
segment 13++; 

else  if (array [row] [col] <224) 
segment 14++; 

else  if (array [row] [col] <240) 
segmentl5++ ; 

else  if (array [row] [col] <=255) 
segment 16++; 
num_byte++ ; 

> 

} 

} 

} 

printf("The  maximum  intensity  in  file  '/,s  is  '/,u  \n" , 
norm_f ile ,max) ; 
peak=max*2040 . 0/255 . 0 ; 
fclose(ptrl) ; 

> 

if (histo== ’ y ' ) 

{ 

printf  ("\n\n\n\t\t\t\tHISTOGRAM  OF  DATA  \n’')  ; 

printf  (''\t\t  (Value  in  parentheses  is  value  on  spiricon 

plot)  \n\n\n"); 

printf ("'/,. 2f  \tpercent  of  the  values  lie  between 
0(0)  and  15(127).  \n" , (float)segmentl/ 

(f loat)num_byte*100 . 0) ; 

printf ("'/,. 2f  \tpercent  of  the  values  lie  between 


16(128)  and  31(255).  \n" , (f loat)segment2/ 

(f loat)num_byte*100 . 0) ; 

printf  ("y, .  2f  \tpercent  of  the  values  lie  between 
32(256)  and  47(382).  \n" , (float) segment3/ 

(f loat)num_byte*100 . 0)  ; 

printf  ("*/,. 2f  \tpercent  of  the  values  lie  between 
48(383)  and  63(510).  \n" , (float) segment4/ 

(f loat)num_byte*100 . 0) ; 

printf ("'/,. 2f  \tpercent  of  the  values  lie  between 
64(511)  and  79(637).  \n" ,  (float) seginent5/ 

(f loat)num_byte*100 . 0) ; 

printf  ("'/,. 2f  \tpercent  of  the  values  lie  between 
80(638)  and  95(765).  \n" , (float)segment6/ 

(f  loat)nuin_byte*100 . 0)  ; 

printf  ("'/,. 2f  \tpercent  of  the  values  lie  between 
96(766)  and  111(892).  \n" , (float) segment7/ 

(f loat)num_byte*100 . 0) ; 

printf  ("*/,.  2f  \tpercent  of  the  values  lie  between 
112(893)  and  127(1020).  \n" , (float) segments/ 

(f loat)num_byte*100 . 0) ; 

printf  ("*/, .2f  \tpercent  of  the  values  lie  between 
128(1021)  and  143(1147).  \n" , (float) segment9/ 

(f loat)num_byte*100 . 0) ; 

printf  ("'/,. 2f  \tpercent  of  the  values  lie  between 
144(1148)  and  159 ( 1275)  .  \n" , (float) segmentlO/ 

(f loat)num_byte* 100 . 0) ; 

printf  ("'/,.  2f  \tpercent  of  the  values  lie  between 
160(1276)  and  175(1402).  \n" , (float ) segment  1 1/ 

(f loat)num_byte*100 . 0) ; 

printf  ("'/,.  2f  \tpercent  of  the  values  lit-  >etween 
176(1403)  and  191(1530).  \n" , (float)segmentl2/ 

(f loat)num_byte*100 . 0) ; 

printf  (”'/,.  2f  \tpercent  of  the  values  lie  between 
192(1531)  and  207(1657).  \n" , (float) segment  13/ 

(float ) num_byte* 100 . 0) ; 

printf  ("'/,. 2f  \tpercent  of  the  values  lie  between 
208(1658)  and  223(1785).  \n" , (float ) segment  14/ 

(f loat)num_byte*100 .0) ; 

printf (''*/,. 2f  \tpercent  of  the  values  lie  between 
224(1786)  and  239(1912).  \n" , (float) segment  15/ 

(f loat)num_byte*100 . 0) ; 

printf ("*/,. 2f  \tpercent  of  the  values  lie  between 


240(1913)  and  255(2040).  \n" , (float) segment  16/ 
(float)num_byte*100.0) ; 

} 

printf (''\n\nEnter  threshold  value  ('/.u  is  peak)"^  \n",peak); 
scant  ("*/,u"  ,&threshold)  ; 

printf ("\nDo  you  wish  to  set  intensities  above  the  threshold  ") ; 
printf ("to  a  single  value,  y  or  n?  \n"); 
scant  ("  */,c"  ,&binary)  ; 
if (binary==’y ' ) 

{ 

printf ("XnEverything  below  y.u  will  be  zero  ".threshold); 
printf  ("and  everything  aoove  y.u  will  equal  2040  \n", 
threshold) ; 

} 

else 

printf  ("\nEverything  below  '/.u  will  be  zero  \n"  , threshold)  ; 
thresh_scale= (f loat )threshold*255 .0/2040.0; 
printf  ("Threshold  equals  ‘/.f  \n\n"  ,thresh_scale)  ; 
sprintf  (thresh.f ile , "'/.s . th3" .filename)  ; 

printf  ("Opening  input  file  '/.s  for  thresholding  \n"  ,r.- rr.:  ile)  ; 
ptrl=f open(norm_f ile , "rb") ; 

printf  ("Opening  output  file  ’/.s  for  thresholding  \n\:  , 

thresh.f ile) ; 

ptr2=f open (thresh_f ile , "wb") ; 

f seek(ptrl ,2,0)  ; 

fwritf (header, 1 ,2,ptr2) ; 

1  or  (c'-unt  1=0  ;  count  1<  si  ice;  count  1++) 

printf  ("reading  slice  #y,u  \n"  ,countl+l)  ; 
f read (array , 1 , length-2*incol ,ptrl) ; 
printf  ("ptrl=y,lu  \n\n"  ,ftell  (ptrl) ) ; 
f or (row=0 ; row<inrow-2 ; row++) 

{ 

f  or  (col=0  ;  col<iriCol  ;col++) 

{ 

if (binary== ' y ' ) 

{ 

if (array [row] [col] <thresh_scale) 
array [row] [col J  =0 ; 

else 

array [row]  [col]=255; 

} 


else 

if (array [row]  [col] <thresh_scale) 
array [row] [col] =0 ; 

} 

} 

} 

printf  ("Writing  slice  #y,u  \n"  , countl+l)  ; 
fwrite(array , 1 ,length-2*incol,ptr2) ; 
printf  ("ptr2='/,lu  \n\n"  ,ftell(ptr2))  ; 

} 

f read(array , 1 , 1930, ptrl) ; 
fwrite(array , 1 , 1930 ,ptr2) ; 
f close(ptrl) ; 
f close(ptr2) ; 

> 

/♦BEGINNING  MAXBYTE  AND  SUMMATION  ROUTINE*/ 

printf ("\n\nDo  you  wish  to  find  the  coordinates  of  the  peak 
intensity,  y  or  n?  \n") ; 
scanf("  y,c"  ,&maxbyte)  ; 
while (maxby t  e== ' y ’ ) 

{ 

max=0 ; 
sum=0 ; 
num_byte=0 ; 

printf ("\nEnter  directory  path  and  file  name  of  input  file 
:  \n\n"); 

scanf  ("'/.s"  ,inax_f  ile)  ; 

printf ("\nEnter  record  number  \n") ; 

scanf  ("‘/.u"  ,&record_num)  ; 

printf  ("XnOpening  input  file  '/.s  \n"  ,max_f  ile)  ; 
if (  (ptrl=f open (max_f ile ,"rb") )==NULL) 

{ 

printf  ("XnCan't  open  input  file  '/.s  \n"  ,max_f  ile)  ; 
exit  0  ; 

} 

/♦MOVE  POINTER  TO  DESIRED  RECORD  IN  FILE*/ 
off set =(record_num-l) *263052+2; 
f seek (ptrl , offset ,0) ; 


printf  ("\nCalculating  maximum  intent-ity  \n\n")  ; 
f or (count  1=0 ; count l<slice; count 1++) 

{ 

printf  ("reading  slice  #'/,u  \n"  ,countl  +  l)  ; 
f read(array , 1 , iength-2*incol ,ptrl) ; 
printf  ("ptrl='/,lu  \n\n"  ,f tell (ptrl)  )  ; 
f or(row=0 ;row<inrow-2 ; row++) 

{ 

f or (col=0 ; col<incol ;col++) 

{ 

sum+=array [row] [col] ; 
num_byte++ ; 

max=(max>array  [row]  [col]  )  max  :  array  [row]  [col]  ; 
if(peakmax  !=  max) 

{ 

peakmax=max ; 
peakbyte=num_byte ; 

> 

} 

} 

> 

fclose(ptrl) ; 
peak=max*2040 . 0/255 . 0  ; 
peakrow=(peakbyte/(long) incol) ; 
peakcol=  (peakbyte  '/,  (long)  incol) -1 ; 

printf  ("The  siim  of  the  first  '/,lu  intensities  is  ‘/.lu  \n\n", 
num_byte , sura) ; 

printf  ("The  peak  intensity  of  */,u  occurs  at  row  '/,u  and  \n",max, 
peakrow; ; 

printf  ("column  /,u  in  file  */.s.  \n\n\n"  .peakcol  ,max_f  ile)  ; 
printf  ("This  corresponds  to  a  peak  of  '/.u  at  ".peak); 
printf  ("coordinates  (’/.u./.u)  on  Spiricon  plot.  \n\n"  ,  peakcol , 
512-peakrow) ; 

printf ("Do  you  wish  to  find  the  peak  and  sum  in  another 
f lie ,  y  or  n?  \n" ) ; 
scanf("  */,c"  , &maxbyte)  ; 

> 

/=»^BEGINNING  RANDOM  ACCESS  ROUTINE*/ 

printf ("\n\nDo  you  wish  to  look  at  a  particular  aiea  of 
data ,  y  or  n'^  \n")  ; 


scanf("  Xc" ,&random) ; 
while (random== ’ y ' ) 

{ 

printf ("\nEnter  beginning  row  and  column  starting  at  0,0 

(i.e.  24,105  na"); 

scant  ("'/,u ,*/,u"  ,&begin.row ,&begi.'...col)  ; 

print!  ("The  beginning  row  and  column  is  ‘/,u,*/,u  \n\n"  , begin_row , 
begin_col)  ; 

printf ("Enter  ending  row  and  column  \n"); 
scant  ("*/,u,*/,u"  ,&end_row,&end_col)  ; 

printf  ("The  ending  row  and  column  is  '/,u,’/,u\n\n"  ,end_row,end_col)  ; 
ptr_off set =( long) begin. row* (long) incol+begin_col+2 ; 
r andom_length=( end. row*incol+end_col+ 2 ) -ptr .offset ; 
random_row=random. length/ incol ; 
random.  col=random.length‘/,incol ; 

printf  ("The  pointer  offset  =  */,lu  \n"  ,ptr.of  f  set)  ; 
printf  ("The  random  length  =  */.u  \n\n"  ,  random. length)  ; 
if (random. length<length) 

{ 

printf ( "Enter  directory  path  and  file  name  of  input  file 
:  \n\n"); 

scant  ("'/.s"  ,random_f  ile.access) ; 
printf ("\n0pening  input  file  /.s  for  random  access  \n\n", 
random.f ile.access) ; 

if (  (ptrl=f open (raindom.f ile.access, "rb") )==NULL) 

{ 

printf  ("Can' t  open  input  file  */,s  \n"  ,  random.!  ile.access)  ; 
exit  0  ; 

} 

printf  ("Below  is  the  16  column  array  of  values  from  file  '/.s 
\n\n" .random. file. access) ; 
f seek(ptrl ,ptr_offset ,0) ; 
tread (array , 1 , random.length ,ptrl) ; 
f or (row=0 ; row<random.row ; row++) 

{ 

for(cnl=0; CO 1< incol ;col++) 

{ 

if  (col'/,16==0) 

printf ("\n") : 
pr  intf  ("'/,5u"  ,  array  [row]  [col]  )  ; 

} 

1 


7;{ 


f oi(col=0 ; col<random_col ; col++) 

{ 

iJ(col'/.16==0) 
pnntf  ("\n")  ; 

printf  ("y.5u"  ,  array  [random_row]  [col]  )  ; 

} 

printf  (''\n\n'')  ; 
f close (ptrl ) ; 

printf  ("Do  you  wish  to  look  at  another  area,  y  or  n'i’  \n"); 
scanf("  y.c"  ,&random)  ; 

} 


printf  ("Maximum  number  of  allowed  rows  is  '/.d  and  the 
maximum  " , inrow) ; 

printf  ("number  of  allowed  columns  is  ‘/.d  \n",incol); 

} 

> 

printf  ("\nDo  you  wish  to  c  lange  something,  y  or  n"''  \n")  ; 
scant  i'  "  y.c"  ,&program)  ; 


;  (1  ist  .s  ino) 

1.  g  n  e  d  char  i  i  s  t  [  J  ; 
:gned  siz>'; 


/^SORTING  FUnCTIQN+/ 


■'U-'  -  ’ 

::  lit  =  >',  ;  c:  •, <s  1  ze-  1  ;  out  +  +  ) 


•  Ut  +  l  ;  irKsize;  in  +  +  ) 


1  s  f,  L  o  u  t  '  >  ]  1  s  t  [  1  n  ] 


tempi =1  i  St [ in]  ; 

1  1  s t [ ] n ]  =  1 1 st  i out] 
1  1  E  t  r out ] =temp 1  ; 

} 


Appendix  C.  Program  PREPROC.C 


This  program  is  similar  to  logSxS.c,  however,  it  is  (h'sigiied  to  ])repro(('ss  mul¬ 
tiple  files  without  user  intervention.  The  [)rogra!n  consists  of  a  reduction  routine, 
ni'-dian  filter,  normah.zation,  am!  i hreslioiding  oj)tions.  t'pon  enteiing  the  program, 
the  user  is  ask  if  the  data  files  are  test  or  t('m|)late  files.  J  he  only  difference  lies 
in  the  names  of  tlu'  input  and  output  files,  h'or  a  test  fih'.  tin'  input  lih  name  may 
he  dtstl  for  example,  while  a  temi)late  file  would  he  named  2d('gl.  1  he  output  liles 
are  a])pended  witfi  .‘ioG  (i.(‘.  2t st  1 .2oG).  The  2  r('presents  eitln'r  the  angle  of  ar¬ 
rival  (.\0.\)  for  the  template  file  or  the  numher  of  the  t('.st  fih'.  'The  1  represents  tin' 
h'ligth  of  the  fiber. 

Tlu'  program  next  prompts  the  user  for  the  init  ial  and  ending  ,\().\  for  tr-mplate 
files,  or  initial  and  ending  numlx'r  for  t<'st  files,  for  example',  if  tlu'  user  has  taken 
data  for  AOA  from  3°  to  15°  in  om'  de'gree  increments,  tlu'ii  h"  ente-i's  2)  and  1.5  for 
initial  and  ending  .AOA  r('sp('ct  iv<'ly.  He  tin'll  ente'rs  tin'  h'ligth  of  the  liln'r  (i.('.  2  for 
a  two  inch  filler)  at  the  next  prompt,  and  the  program  will  automat  icall\'  preproce'ss 
files  3d<'g2  through  15d('g2.  It  will  locate  and  j>r('proc<'ss  tln-sc'  liles  ai'cording  to  the 
inlormation  the  user  supplies.  Tin'  options  are'  log  compression,  a  histejgram  to  In'lp 
d'termine  threshold,  thre'sholel  value',  hinarizat  iein.  ami  input  ami  eaitput  elirccteery. 
1  lie'  preigram  edle'rs  an  i'lpiit  ami  out  put  elire'e  t e)ry  sugge'st  iejii.  aiiel  i  iio  u^er  e  an  ■'e'h'e  t 
tlii"  h\'  typing  'a',  eu  the'  use'r  e-an  e'ute'r  his  envn  input  ami  output  elin'elorw  Tin’ 
progra  III  use-s  a  S\\  hi  ('ll  si  at  e-me'iit  tee  eke  this  see  that  ot  In'r  eliie'et  eerie 's  ma\'  he'  ae  hh’il 
to  the  pieegi'ain  at  any  time’  just  hy  ajepe'mhng  tln'in  to  the'  ,S\\1  It'll  state’incni  in 
t  he  souree'  e'oeh'.  1  Ins  ’ 'ill  he'lj)  save'  typing  t  he- e'lit  ire- eliicel  oiy  path  e'ae  h  time’. 

In  a.elditieeii  to  I  In'  pre’|»re)e'e'sse'el  data  euitjeiit  file's  (whieh  are'  a  [epe'iie  Ici  1  with 
.2  it)  I  t  lie'l  l'  is  a  file'  l  a I li'el  1  f.M  1’  1 1  IS  I  .  |){)( '  lor  t e'lnplai  i'  I i h 's  c  u'  1  f.S  1  1 1  IS  1  , 1 )( )( ' 
loi  I  i.'st  h h 's  \\  1 1  n  h  IS  al  ways  phu'i'il  111  t  In'  1  i•..\l  1  ’  1  ..\  1  f.  oi'  1  f.S  1  el  1 1 1  'i  t  or\’  l  e’Spi'c- 
ti'.'l\.  11ns  lilc  contains  tin'  jn-ak  \aha'  lii'loie'  ami  altei  ina  i  na  1  i/'a  i  na  i  as  well  as 


it's  (•()()r(liiiat('.s,  tlic  sum  of  iho  iiilt  nsitics,  and  a  histogram  for  cacli  data  file.  I  hc 
1)0  ginning  of  this  fih'  indicatc's  what  program  oj)tions  worn  soh'ctisl. 

In  summary,  thr  program  reduces  each  512  x  512  file  to  a  250  x  250  file.  It  also 
('xeentes  multiph'  fih's  autotnat  n  ally,  ddie  fih's  must  Ix'  in  oik'  d('gree  increments  for 
template  fdes  and  one  numher  increments  for  test  files.  For  example,  using  a  one  inch 
filxT,  the  user  may  take  test  data  at  3°,  19°,  and  12°,  howe\(‘r.  he  w<)u!d  name  thes<> 
files  Otst,  1,  Itstl,  and  2tstl  respectively.  When  all  processing  is  completed  (i.e.  the 
exec’ution  of  Sdft,  t(  mpldh  .c,  and  (listaiicr.c),  tin*  program  ilisidiir:  .c  will  classify  the 
test  fih's,  hop('fully  corix'ctly.  d  h<'  |)rogram  also  riMjuires  expanded  memoiA'  to  run 
sinc('  it  reads  and  writes  to  t('mporary  fih's  in  drive  1).  Fach  file  takes  approximatidy 
1  minute  and  2(J  sc'conds  (o  pre|)i'ocess. 


/♦Program  PREPROC.C*/ 


/♦Written  by  Capt  John  W.  Welker^/ 


/♦This  prograia  converts  a  512  x  512  spiricon  file  into  a  256  x  256^/ 
/♦file  by  grabbing  only  the  center  256  x  256  pixels  in  the  spiricon*/ 
/♦file.  ♦/ 

/♦This  program  performs  the  pre-processing  on  a  Spiricon  image  file-^/ 
/♦This  pre-processing  consists  of  a  median  filter  operation  which^/ 
/♦removes  spikes  in  the  data.  The  data  is  then  normalized  to  the  ♦/ 
/♦maximum  of  255.  The  data  may  be  log  compressed  to  increase  the^/ 
/♦lower  intensities.  A  histogram  of  the  data  can  then  be  examined*/ 
/♦to  determine  the  best  threshold.  Data  below  the  threshold  is*/ 

/♦set  to  zero.  Data  above  the  threshold  is  either  set  to  255.*/ 

/♦or  left  unchanged.  Recursion  means  that  the  median  value  is*/ 
/♦placed  into  the  middle  of  the  3x3  matrix  before  it  calculates*/ 
/♦the  next  median  value.*/ 


/♦COMPILER  DIRECTIVES*/ 
ftinclude  <stdio.h> 

#include  <string.h> 

#include  <math.h> 

/♦DEFINE  MACROS*/ 

#define  inrow  64 
#define  incol  512 

#define  slice  4 

#define  outrow  66 
#define  outcol  256 

/♦DECLARE  POINTERS  TO  FILE  STRUCTURE*/ 

FILE  *ptr_in , ♦terap_ptr , *templ_ptr,*ptr_out , ♦histogram_ptr ; 

/♦START  MAIN  PROGRAM*/ 
void  mainO 


{ 


/^DECLARE  LOCAL  VARIABLES*/ 

unsigned  char  f ilename.in [256] ,f ilename_out [256] , temp_f ile [256] ; 
unsigned  char  templ_f ile[256] ,directory_in [256] ,directory_out [256] ; 
unsigned  char  hi stogram.f ile [256] .temp [9] ; 

unsigned  char  array_in [inrow] [incol] ,array_out [outrow] [outcol] ; 
unsigned  long  Iength_in=0,length=0,sujii=0,num_byte=0,peakbyte=0; 
unsigned  long  slice_sum=0 ,norm_sum=0 ; 

unsigned  long  segment  1 , segment 2 .segments, segment4 , segments , segments ; 
unsigned  long  segment/, segments, segment9, segment  10 .segment  11 ; 
unsigned  long  segment  12 , segment 13, segment  14, segment  15 , segment  16 ; 
unsigned  row , col , count , count  1 , count 2 .count 3 , c , r ,a=0 ,max=0 ,peakmax=0 ; 
unsigned  peakrow=0 ,peakcol=0 ,begin_degree=0 , end_degree=0 ; 
unsigned  f iber_length=0 ,begin_test=0 , end_test=0 ,norm_max=0 ; 
unsigned  threshold=0; 

char  test , log_ scale , hist o .thresh .binary ; 
float  norm_scale=0 . 0 ; 

length_in=inrow*incol ; 
length=outrow*outcol ; 

strcpy  (temp_f  ile,  "d :  Wtemp")  ; 
strcpy  ( tempi _file ,  "d :  Wtempl")  ; 


printf ("\n\nAre  these  test  files,  y  or  n'i’  "); 
scanfC  y.c"  ,&test)  ; 
printf ("\n\n") ; 
if (test=='y' ) 

printf ("What  is  the  number  of  your  initial  test  file?  ") ; 
scanf  ("'/,u"  ,&begin_test)  ; 
printf ("\n\n") ; 

printf  ("What  is  the  number  of  your  last  test  file'?’  ")  ; 

scanf  ("'/.u"  ,&end_test)  ; 

printf ("\n\n" ) ; 

begin_degree=begin_test ; 

end_degree=end_test ; 

strcpy (histogram_f lie , "c ; \\spir\\test\\testhist . doc" ) ; 

> 

else 


TS 


{ 

printf("Your  first  input  file  corresponds  to  what  angle  of 

incidence?  "); 

scanf  ("'/,u"  ,&begin_degree)  ; 

printf ("\n\n") ; 

printf("Your  last  input  file  corresponds  to  what  angle  of 
incidence?  ") ; 

scanf  ("'/.u"  ,&end_degree)  ; 
printf ("\n\n") ; 

strcpy  Chistograjn.f  ile,  "c .  WspirWteraplateWtemphist  .doc")  ; 

} 

printf ("How  long  is  your  fiber?  ") ; 
scanf  ("*/,u"  ,&fiber_length)  ; 
printf ("\n\n") ; 

printf ("Do  you  wish  to  perform  a  log  compression  of  the 

data,  y  or  n?  ") ; 

scanf  ("  y.c"  ,&log_scale)  ; 

printf ("\n\n") ; 

if (log_scale== 'y ' ) 

fprintf (histograra_ptr,"A  log  compression  was  performed  on  this 
data.  \n\n"); 

printf ("Do  you  wish  to  see  a  histogram  of  the  data  ") ; 
printf ("for  thresholding  purposes,  y  or  n?  ") ; 
scanf  ("  */,c"  ,&histo) ; 
printf ("\n\n") ; 

printf  ("Do  you  wish  to  set  a  threshold,  y  or  n'r’  "); 
scanf  ("  */,c"  ,&thresh)  ; 
printf ("\n\n") ; 
if (thresh== ' y ' ) 

{ 

printf ("Enter  threshold  value  (255  is  peak).  ") ; 
scanf  ("'/.u"  ,&threshold)  ; 
printf ("\n\n") ; 

fprintf  (histogram_ptr  ."Threshold  was  set  at  ’/,u  \n\n"  , threshold) 
printf ("Do  you  wish  to  set  intensities  above  the  threshold  "); 
printf  ("to  a  single  value,  y  or  n"^  \n"); 
scanf  ("  y.c"  ,&binary)  ; 
if (binary=='y ') 

{ 

printf  ("Everything  below  '/.u  will  be  zero  ",  threshold)  ; 
printf  ("and  everything  above  or  equal  to  '/.u  will  equal 
255  \n\n" , threshold) ; 


f printf  (histograjn_ptr ,  "This  data  was  binarized  \n\n"); 

} 

else 

printf  ("Everything  below  ‘/,u  will  be  zero  \n"  .threshold)  ; 

} 

/♦ENTER  FILES  AND  VARIABLE  VALUES  FROM  KEYBOARD*/ 
printf ("\nEnter  directory  path  of  input  file  or;  \n\n")  ; 
printf ("Enter  'a'  for  c : \\spir\\rawdata\\  "); 

scanf  ("'/.s"  ,directory_in) ; 
printf ("\n\n") ; 
switch  (*directory_in) 

{ 

case  ' a’ ; 

strcpy ( directory. in, "c ; \\spir\\rawdata\\")  ; 
break ; 

} 

printf ("\nEnter  directory  path  of  output  file  or;  \n\n")  ; 

printf ("Enter  ’a'  for  c ; \\spir\\procdata\\  "); 

scanf  ("'/.s"  ,directory_out)  ; 

printf ("\n\n") ; 

switch  (*directory_out) 

case  ' a ' ; 

strcpy (directory.out , "c ; \\spir\\procdata\\") ; 
break ; 

> 

printf  ("Opening  output  file  */s  \n\n"  ,histogram_f  ile)  ; 
if (  (histogram_ptr=f open (histogram. file ,"wb") )==NULL) 

{ 

printf  ("Can't  open  output  file  V.s  \n"  .histogram.!  ile)  ; 
exit  0  ; 

} 

f or (count=begin_degree ; count<=end_degree ; count++) 

{ 

if (test== ' y ' ) 

{ 

sprintf  (f  ilename.in  ,  "*/,sy,utst*/,u"  .directory. in,  count ,  f  ibcr.length)  ; 


SO 


sprintf  (f  ilename_out ,  "'/.sXutstXu .  256"  .directory^out ,  count , 
f iber_length) ; 

} 

else 

{ 

sprintf  (f  ilenaine_in , "*/.s‘/,udeg*/,u"  ,directory_in , count ,  f  iber_length)  ; 
sprintf  (f  ilenaine_out ,  "*/,s%udegXu .  256"  .directory^out ,  count , 
f iber_length) ; 

} 

/♦OPEN  INPUT  AND  OUTPUT  FILES*/ 

printf  ("\nOpening  input  file  */,s  \n"  ,f  ilenajne^in)  ; 
if (  (ptr_in=f open (filename. in, "rb") )==NULL) 

{ 

printf  ("Can't  open  input  file  ‘/.s  \n"  ,f  ilename.in)  ; 
exit ()  ; 

} 

printf  ("Opening  output  file  */,s  \n\n"  ,temp_f  ile)  ; 
if (  (temp_ptr=f open (temp.f ile , "wb") )==NULL) 

{ 

printf  ("Can't  open  output  file  */,s  \n"  ,temp_f  ile)  ; 
exit  0  ; 

> 

fseek(ptr_in, 65538,0) ; 

for (count  1=0 ; count 1< si ice ; count 1++) 

{ 

printf  ("reading  slice  #'/,u  \n\n"  ,countl  +  l)  ; 
f read(array_in, 1 ,length_in,ptr_in) ; 
f or (row=0 ; row< inrow ; row++) 

{ 

for(col=0; col<incol ; col++) 

{ 

if(col>127  &&  col<384) 

{ 

array.out [row] [col-128] =array_in[row] [col] ; 

} 

> 

} 

printf  ("\n\nWriting  slice  #*/,u  \n\n"  .count  1  +  1)  ; 
f write (array.out , 1 , length-2+outcol , temp.ptr) ; 

} 

f closa(ptr_in) ; 


SI 


fclose(temp_ptr) ; 


/♦BEGIN  MEDIAN  FILTER  PROCESSING+/ 

/♦OPEN  INPUT  AND  OUTPUT  FILES^/ 

printf  ("\nOpening  input  file  '/,s  \n"  ,temp_f  ile)  ; 

if (  (temp_ptr=f open(temp_f ile, "rb") )==NULL) 

printf  ("Can’t  open  input  file  */,s  \n"  ,  temp_f  ile)  ; 
exit  ()  ; 

} 

printf  ("Opening  output  file  '/.s  \n\n"  ,templ_f  ile)  ; 
if (  (templ_ptr=f open (tempi _f ile , "wb"))==NULL) 

printf  ("Can’t  open  output  file  '/,s  \n"  ,templ_f  ile)  ; 
exit  0  ; 

} 


sum=0 ; 
max=0 ; 
num_byte=0 ; 
peakmax=0 ; 
peakbyte=0 ; 
peakrow=0 ; 
peakcol=0 ; 

for (count  1=0 ; count 1< si ice ; count 1++) 

{ 

if (countl==0)  /♦FIRST  PASS  ONLY* 

{ 

printf  ("reading  slice  #*/,u  \n"  ,countl  +  l)  ; 

/♦FIRST  READ  INTO  ARRAY*/ 
fread(array_out , 1 , length-2+outcol , temp_ptr) ; 
printf  ("temp_ptr  =  */.lu  \n\n"  ,  f  tell  (temp_ptr)  )  ; 
r=0 ; 

for (count 2=0 ; count 2<out row-4 ; count 2++) 

/♦FIRST  MEDIAN  FILTER  LOOP*/ 

{ 

c=0  ; 

for (count 3=0 ; count3<outcol-2 ; count3++) 


S2 


{ 

a=0 ; 

for(row=r;row<r+3;row++)  /^SELECT  VALUES  FOR*/ 

{  /*3  X  3  ARRAY*/ 

for(col=c;col<c+3;col++)/*AND  PREPARE  FOR  SORT*/ 

{ 

temp [a] =array_out [row] [col] ; 
a++; 

} 

} 

sort (temp, a) ;  /*SORT  3  X  3  ARRAY*/ 

C++; 

array_out[r+l] [c]=temp[a/2] ;/*PLACE  MEDIAN  VALUE  IN*/ 

}  /*CENTER  OF  3  X  3  ARRAY*/ 

r++; 

} 

for(col=0;col<outcol;col++)  /*SET  FIRST  ROW  EQUAL*/ 

array_out [0] [col] =array_out [1] [col] ;  /*T0  SECOND  ROW*/ 

f or (row=0 ; row< out row-4 ; row++) 

{  /*RESET  FIRST  AND  LAST*/ 

array_out [row] [0] =array_out [row] [l] ;  /*C0LUMN*/ 
array_out [row] [outcol-l]=array  out  [row]  [outcol-2] ; 

> 

printf ("Calculating  first  slice  intensity  and  maximum  \n\n"); 
slice_sum=0 ; 

for(row=0;row<outrow-4;row++)  /*CALCULATE  MAXIMUM  VALUE*/ 

{  /*AND  INTENSITY*/ 

for(col=0; col<outcol ; col++) 

{ 

sum+=array_out [row] [col] ; 
si ice_sum+=array_out [row] [col]  ; 
num_byte++ ; 
max  = 

(max  >  array_out  [row]  [col] )  '?’max  :  array_out  [row]  [col]  ; 
if(peakmax  !=  max) 

{ 

peakmax=max ; 
peakbyte=num_byte ; 

} 

} 

} 

printf  ("The  sum  of  this  slice  is  V.lu  \n\n"  ,  slice_sum)  ; 


S3 


printf  ("Writing  slice  #y,u  \n"  .count  1+1)  ; 
f write (array_out , 1 , length-4*outcol ,templ_ptr) ; 

/♦WRITE  ALL  BUT  LAST  TWO*/ 

/♦ROWS  TO  OUTPUT  FILE*/ 

printf  ("templ_ptr  =  */,lu  \n\n"  ,ftell(tenipl_ptr) )  ; 

} 

else  /*BEGIN  ALL  SUCCESSIVE*/ 

{  /*PASSES  THROUGH  FILTER*/ 

f or (row=0 ; row<2 ; row++) 

{ 

f or (col=0 ; col<outcol ; C0I++)  /*LAST  TWO  ROWS  ARE  NOW*/ 

{  /*FIRST  TWO  ROWS*/ 

array_out [row] [col] =array_out [outrow-4+row] [col]  ; 

} 

} 

printf  ("reading  slice  #'/,u  \n"  ,countl+l)  ; 
f read(&ari ay_out [2] [0] , 1 , length-2*outcol ,temp_ptr) ; 

/♦READ  IN  SECOND  SLICE  OF*/ 

/♦DATA  STARTING  IN  THIRD*/ 

/♦ROW  OF  ARRAY*/ 

printf  ("temp.ptr  =  '/.lu  \n\n"  ,ftell(temp_ptr) )  ; 
r=0; 

for (count 2=0 ; count 2 < out row- 2 ; count 2++) 

{  /♦SUCCESSIVE  MEDIAN  FILTER*/ 

c=0;  /*L00PS*/ 

for (count 3=0 ; count3<outcol-2 ; count 3++) 

{ 

a=0 ; 

for(row=r;row<r+3;row++)  /*SELECT  VALUES  FOR  3  X  3*/ 

{  /*ARRAY  AND  PREPARE  FOR*/ 

for (col =c ; col<c+3; C0I++)  /* SORTING*/ 

{ 

temp [a] =array_out [row] [col] ; 
a++; 

} 

> 

sort (temp , a) ;  /*S0RT  3X3  ARRAY*/ 

C++; 

/♦PLACE  MEDIAN  VALUE  INTO*/ 

/♦CENTER  OF  3  X  3  ARRAY*/ 

array_out [r+l]  [c] =temp [a/2]  ; 


} 


r++; 

> 

f or (row=0 ; row<outrow-2 ; row++)  /♦RESET  FIRST  AND  LAST*/ 

{  /*COLUMNS*/ 

array_out [row] [0] =array_out [row] [l] ; 
array_out [row] [outcol-l] =array_out [row] [outcol-2] ; 

} 

if (countl<slice-l)  /*INTERMEDIATE  LOOPS*/ 

{ 

printf ("Calculating  intermediate  slice  intensity  and 
maximum  \n\n") ; 
slice_sum=0 ; 

f or (row=0 ; row<outrow-2 ; row++) 

{ 

f or (col=0 ; col<outcol ; col++) / *CALCULATE  MAXIMUM*/ 

{  /*VALUE  AND  INTENSITY*/ 

sum+=array_out [row] [col] ; 
slice_sum+=array_out [row] [col] ; 
num_byte++ ; 
max  = 

(max  >  array_out [row] [col] )  ?max;array_out [row] [col] ; 
if(peakmax  !=  max) 

{ 

peakmax=max ; 
peakbyte=num_byte ; 

} 

} 

•» 

/ 

printf  ("The  sum  of  th-is  slice  is  '/,lu  \n\n"  ,  slice_sum)  ; 
printf  ("Writing  slice  #*/.u  \n"  ,countl  +  l)  ; 
fwrite(array_out ,  1  ,length-2*outcol  .tempi _ptr'' 

/♦WRITE  TO  OUTPUT  FILE*/ 

/♦DURING  INTERMEDIATE*/ 

/♦PASSES  THROUGH  FILTER*/ 

printf  ("tempi_ptr  =  '/.lu  \n\n"  , ftell  (templ_ptr) )  ; 

> 

else 

{  /*LAST  PASS  THROUGH  IILTER*/ 

/♦RESET  FIRST  AND  LAST  COLUMN*/ 
array_out [outrow-2] [0] =array_out [outrow-2] [l] ; 
array_out [outrow-2]  [outcol-l]  = 


Hr> 


array_out  [outrow-2] Coutcol-2] ; 

f or (col=0 ; col<outcol ; C0I++)  /*RESET  LAST  ROW*/ 

array_out  [outrow-l] [col] =array_out [outrow-2]  [col]  ; 

printf ("Calculating  last  intensity  and  maximum  \n\n") ; 
slice_sum=0 ; 

for(row=0;row<outrow;row++)  /+CALCULATE  MAXIMUM  VALUE*/ 

{  /*AND  IN'^ENSITY  IN  LAST*/ 

f or (col=0 ; col<outcol ; col++)  /♦SLICE*/ 

{ 

sam+=array_out  [row]  [col]  ; 
slice_sum+=array_out [row] [col] ; 
num_byte++ ; 
max  = 

(max  >  array_out [row] [col] )  ?max : array_out [row] [col] ; 
if (peakmax  !=  max) 

{ 

peakmax=max ; 
peakbyte=num_byte ; 

} 

> 

} 

printf  ("The  sum  of  this  slice  is  V.lu  \n\n"  ,  slice,  sum)  ; 
printf  ("Writing  slice  #*/,u  \n'' .count  1+1)  ; 
fwrite(array_out , 1 , length, tempi. ptr) ; 

/♦WRITE  ALL  OF  LAST  SLICE*/ 

/*T0  OUTPUT  FILE*/ 

printf  ("templ_ptr  =  */lu  \n\n"  .ftell  (templ_ptr)  )  ; 

} 

} 

> 

f close(temp_ptr) ;  /*CL0SE  FILES*/ 

f close (templ_ptr) ; 

peakrow= (peakbyte/( long) out col) ; 

peakcol=  (peakbyte  '/,  (long) outcol) -1 ; 

printf  ("The  sum  of  the  first  7.1u  intensities  is  '/,lu  \n\n"  , 
num_byte , sum) ; 

f  printf  (histogram_ptr , ''\n\n\n\n\nThe  sum  of  the  first  '/.lu  intensities 
after  filtering  is  '/,lu  \n\n"  ,num_byte,  sum)  ; 

f  printf  (histogram_ptr ,  "The  peak  intensity  of  '/,u  after  filtering 
occurs  at  \n",max); 

f  printf  (histogram_ptr ,  "row  ‘/,u  and  column  '/.u  in  file  '/,s.  \n\n"  , 
peakrow .peakcol ,f ilename.out) ; 


/♦BEGINNING  »'QRMALIZATION  ROUTINE*/ 

norm_sum=0 ; 
norm_max=0 ; 
nuin_byte=0 ; 
segment  1=0 ; 
segment 2=C ; 
segment3=0 ; 
segment4=0 ; 
segment5=0 ; 
segment6=0 ; 
segment7=0 ; 
segment8=0 ; 
segment9=0 ; 
segment  10=0 ; 
segmentll=0; 
segmentl2=0; 
segmentl3=0 ; 
segmentl4=0 ; 
segmentl5=0; 
segment  16=0; 
if (log_scnle== ’ y ' ) 
norm_scale=255 . 5 /log ( (double) max) ; 
else 

norm_scale=255 . 5/ Cf loat)max ; 

printf  ("Scaling  factor  is  |',f  \n\n”  ,norra_scale)  ; 

printf  ("\n0pening  input  file  */,s  for  normalizationXn"  ,templ_f  ile 

if (  (templ_ptr=f open (tempi _f lie , "rb") )==NULL) 

{ 

printf  ("Can' t  open  input  file  '/.s  \n"  ,templ_f  ile)  ; 
exitO  ; 

} 

printf  ("Opening  output  file  */.s  \n\n"  ,f  ilename_out)  ; 
if (  (ptr_out=fopen(f ilename_out , "wb") )=-NULL) 

{ 

printf  ("Can  '  t  open  input  file  */,3  \n"  ,f  ilename_out)  ; 
exit  0  ; 

} 

for(countl  =  0;countl<slicc;  count  1+-'-) 

{ 
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printf  ("reading  slice  #*/,u  \n"  ,countl+l)  ; 

fread(array_out , 1 ,length-2*outcol ,templ_ptr) ; 

printf  (”templ_ptr='/,lu  \n\n" ,ftell(tenipl_ptr) ) ; 

printf  ("Value  of  array.out  [1]  [1]  before  scaling  is  '/.u  \n", 

array_out  [1]  [1] ) ; 

f or(row=0;row<outrow-2 ; row++) 

-C 

for(col=0; col<outcol ; col++) 

{ 

if (log_scale==^y ^ ) 

{ 

array_out [row] [col] = 

(double)norm_scale*(log( (double) array_out [row] [col] ) ) ; 
if(row==l  &&  col==l) 

printf ("Value  of  array_out  [1] [1]  after  log 
compression  is  */,u  \n"  ,  array_out  [l]  [1]  )  ; 

> 

else 

array_out [row] [col] *=norm_scale ; 
if(row==l  &&  col==l) 

printf ("Value  of  array.out [l] [1]  after  scaling 
is  y,u  \n\n"  ,  array.out  [1]  [l]  ) ; 
norm_sum+=array_out [row] [col] ; 
norm_max= 

(norm_max>array_out  [row]  [col]  )  ?norin_max :  array.out  [row]  [col]  ; 
num_byte++ ; 
if (histo=='y') 

{ 

if (array_cut [row] [col]<16) 
segment 1++ ; 

else  if (array_out [row] [col] <32) 
segment2++ ; 

else  if (array_out [row] [col] <48) 
segment3++ ; 

else  if (array_out [row] [col] <64) 
segment4++; 

else  if (array_out [row] [col] <80) 
segment5++ ; 

else  if (array_out [row] [col] <96) 
segment6++ ; 

else  if (array_out[row] [col]<112) 
segraent7++ ; 
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else  if (array_out [row] [col] <128) 
segmeiit8++ ; 

else  if (array_out [row] [col] <144) 
segment9++ ; 

else  if (array_out [row] [col] <160) 
segment 10++ ; 

else  if (array_out [row] [col] <176) 
segment 11++; 

else  if (array_out [row] [col] <192) 
segment 12++ ; 

else  if (array_out [row] [col] <208) 
segment 13++ ; 

else  if (array_out [row] [col] <224) 
segment 14++ ; 

else  if (array_out [row] [col] <240) 
segment 15++ ; 

else  if (array_out [row] [col]<=255) 
segment 16++ ; 

} 

if (thresh=='y' ) 

{ 

if (binary=='y ' ) 

{ 

if (array.out [row] [col] <threshold) 
array_out  [row]  [col] =0 ; 

else 

array_out [row] [col] =255 ; 

} 

else 

{ 

if (array_out  [row] [col] <threshold) 
array_out [row] [col] =0 ; 

> 

} 

} 

} 

printf  ("Writing  slice  tf'/.u  \n"  ,countl  +  l)  ; 
fwrite(array_out , 1 , length-2*outcol ,ptr_out) ; 
printf  ("ptr_out=‘/,lu  \n\n" ,ftell(ptr_out) )  ; 

} 

f printf (histogram_ptr , "The  sum  of  the  first  '/.lu  intensities 
after  normalization  is  '/.lu  \n"  ,num_byte,norm_sum)  ; 
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fprintf (histogram_ptr , "The  normalized  maximum  in 
file  y,s  is  */,u  \n\n"  ,  f  ilename_out  ,norm_max)  ; 
if (histo== 'y ' ) 


{ 

fprintf (histogram_ptr , "\n\t\t\t\tHISTOGRAM  OF  DATA  \n"); 
fprintf (histogram_ptr, "\t\t (Value  in  parentneses  is  value 
on  spiricon  plot)  \n\n\n'')  ; 

fprintf  (histogram„ptr 2f  \tpercent  of  the  values  lie 
between  0(0)  and  15(127).  Xn” , 

(f 1 oat) segment  1/ (float )num_ byte* 100 . 0) ; 

fprintf (histogram_ptr 2f  \tpercent  of  the  values  lie 

between  16(128)  and  31(255).  \n" , 

(float) segment2/ (f loat)num_byte*100 . 0) ; 

fprintf  (histogram_ptr,  "*/,.2f  \tpercent  of  the  values  lie 

between  32(256)  and  47(382).  \n" , 

(float) segment3/ (float)num_byte*100.0) ; 

fprintf (histogram_ptr,''*/,.2f  \tpe'"cent  of  the  values  lie 

between  48(383)  and  63(510).  \n" , 

(f loat) segment4/ (f loat)num_byte*100 . 0) ; 

fprintf (histc'gram_ptr,"*/,.2f  \tpercent  of  the  values  lie 

between  64(511)  and  79(637).  Xn", 

(f 1 oat) segment 5/ (float )num_byte* 100 .0) ; 

fprintf  (histogram.ptr, ''*/,.2f  Xtpercent  of  the  ''alues  lie 

between  80(638)  and  95(765).  Xn" , 

(float) segment 6/ (f loat)num_byte*100 . 0) ; 

fprintf  (histogram_ptr 2f  Xtpercent  of  the  values  lie 

between  96(766)  and  111(892).  Xn" , 

(float) segment7/ (f loat)num_byte*100 .0) ; 

fprintf (histogram_ptr,"y,.2f  Xtpercent  of  the  values  lie 

between  112(893)  and  127(1020).  Xn" , 

(float) segment8/ (float)num_byte*100 .0) ; 

fprintf (histogram_ptr,"y,.2f  Xtpercent  of  the  values  lie 

between  128(1021)  and  143(1147).  Xn" , 

(float) segment9/ (f loat)num_byte*100 . 0) ; 

fprintf (histogram_ptr,"y,.2f  Xtpercent  of  the  values  lie 

between  144(1148)  and  159(1275).  Xn" , 

(float) segmentlO/ (f loat)num_byte*100 .0) ; 

fprintf (histograun_ptr,"y,.2f  Xtpercent  of  the  values  lie 

between  160(1276)  and  175(1402).  Xn" , 

(float) segment  11/ (float )num_byte* 100 .0) ; 

fprintf (histogram_ptr,"y,.2f  Xtpercent  of  the  values  lie 

between  176(1403)  and  191(1530).  Xn" , 
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(float)  segmentl2/ (f  loat)nm_byte*100 .0) ; 

fprintf  (histogram_ptr,"*/, .2f  \tpercent  of  the  values  lie 

between  192(1531)  and  207(1657).  \n" , 

(float) segmentl3/ (f loat)num_byte*100 . 0) ; 

fprintf (histograiii_ptr,"'/..2f  \tpercent  of  the  values  lie 

between  208(1658)  and  223(1785).  \n" , 

(float) segment  14/ (float )num_byte* 100 .0) ; 

fprintf (histogram_ptr, 2f  \tpercent  of  the  values  lie 

between  224(1786)  and  239(1912).  \n" , 

(float) segmentl5/ (f loat)num_byte*100 .0) ; 

fprintf (histogram_ptr,"'/,.2f  \tpercent  of  tV''^  values  lie 

between  240(1913)  and  255(2040).  \n" , 

(float) segment  16/ (f loat)num_byte*100 .0) ; 

> 

f close(templ_ptr) ; 
f close(ptr_out) ; 

} 

f close(histogram_ptr) ; 

> 

sort (list .size) 
unsigned  char  list[]; 

unsigned  size;  /♦SORTING  FUNCTION* 

{ 

int  out , in, tempi ; 
for(out=0;out<size-l ;out++) 

{ 

f or (in=out+l ; in<size; in++) 

{ 

if (list [out] >list  [in] ) 

{ 

templ=list  [in] ; 
list [in] =list [out] ; 
list [out] =templ ; 

} 

} 

} 

return ; 

} 


Appendix  D.  Program  TEMPLATE. C 


Thi?  program  is  used  to  extract  the  low  freciiieury  Fourier  componenls  (  1!) 
numbers)  and  create  a  fde  containing  feature  vectors,  each  of  which  represents  the 
sliape  of  the  output  intensity  pattern  for  each  degree  of  incidence'.  I'hc  ini)uts  to 
this  program  are  files  which  were  previously  Fourier  translormed  and  have  the  f  1 
pix'fix  (i.e.  F'TlSdeg.'F 2.5(1).  The  FT  indicates  that  the  file  was  j)reviously  Fourier 
transformed.  The  15  represents  a  1-5°  angle  of  orrival  (AOA).  deg  is  indicat i\e  of  a 
te'inplate  file,  the  3  is  for  a  3  inch  length  of  fiber,  and  the  25G  shows  that  the  file  was 
pre[)rocessed  with  the  program  pnproc.c.  More  information  on  input  iidorinat  ion, 
fih'  tiames,  directory  jtatlis,  and  preprocessing  can  be  found  in  appendix  (’. 

I'he  program  places  the  output  data  (i.e.  featnia'  vectors)  in  the  fih'  TF.M- 
PL.Al’F.DAT  or  TEST.DAT  in  directory  TEMPLATlf  or  FE.s  r  n'spect ively.  E<uh 
row  in  this  file  represents  a  particular  AOA  and  contains  19  columtis.  'I'his  file  is 
strictly  used  for  further  processing  by  the  program  di.clance.c.  Atiotlu'r  outirut  fih'. 
1  EMPLAd  E.DOC’  or  TEST. DOC’,  contains  the  same  itiformation  but  can  Ix'  vic'wed 
on  the  screen  or  .sent  to  a  printer.  d'he.se  files  are  also  contaiiu'd  in  tlu'  diix'ctories 
'I  FMPLATE  and  I'ES'l’  respectively. 

Iti  summary,  the  outirut  files  am  variable  h'ugth.  I'ln'  file  ri'.MPL.Al  E.D.Vf 
cu'  TEST.DAT  is  always  19  columns  wide,  but  tin'  nnmlx'r  of  rows  d('pends  on  tin' 
numbf'r  of  files  proc('ss('d.  This  it. formation  is  neech’d  for  the  ])rogram  ilisln nc(  .r. 


/♦Program  TEMPLATE. C*/ 


/“^Written  by  Capt  John  W.  Welker*/ 


/♦This  program  places  the  fundamental  and  first  three  harmonics*/ 
/♦from  one  or  several  256  x  256  fourier  transform  files  into  a  49*/ 
/♦column  array  where  each  row  represents  a  49  component  feature*/ 
/♦vector.  The  files  must  be  named*/ 

/♦  [ft  (degrees  number)deg(f  iber  length). 256]  Excimple:  ft  lOdegS  .  256  .  ♦/ 
/♦This  is  the  fourier  transform  file  representing*/ 

/♦an  angle  of  incidence  of  10  degrees  and  a  fiber  length  of  3*/ 
/♦inches.  The  number  of  your  initial,  test  file  refers  to  the*/ 
/♦number  preceeding  tst  in  the  first  test  file*/ 

/♦filename[ftOtst (fiber  length). 256]  and  then*/ 

/♦  nurber  them  consecutively*/ 

/*  You  have  your  choice  of  input  and  output  directory*/ 

/♦COMPILER  DIRECTIVES*/ 

#include  <stdio.h> 

#include  <string.h> 

#include  <math.h> 

/♦DEFINE  MACROS*/ 
tfdefine  inrow  7 
#define  incol  256 

/♦DECLARE  POINTERS  TO  FILE  STRUCTURE*/ 

FILE  *ptr_in , *template_data_ptr , *template_doc_ptr ; 

/♦START  MAIN  PROGRAM*/ 
void  mainO 
{ 

/♦DECLARE  LOCAL  VARIABLES*/ 

unsigned  char  f ilename_in[256] ,filename_data_out [256] ; 
unsigned  r-har  f  ilename_doc_out  [256]  ,  vector  [49]  ; 
unsigned  char  directory_in[256] ,directory_out [256] ; 


unsigned  char  array_in [inrow] [incol] ; 
unsigned  long  length_in=0; 
unsigned  row, col .count ,a=0; 

unsigned  begin_degree=0 , end_degree=0 ,begin_test=0 , end_test=0 ; 
unsigned  f iber_length=0 ; 
char  test; 

length_in=inrow*incol ; 

/+ENTER  FILES  AND  VARIABLE  VALUES  FPOM  KEYBOARD*/ 
printf ("\n\nAre  these  test  files,  y  or  n?  ") ; 
scanf("  */c",&test); 
printf ("\n\n") ; 
if (test== 'y ' ) 

{ 

printf  ("What  is  the  number  of  your  initial  test  file'?’  ")  ; 

scanf  ("*/,u"  ,&begin_test)  ; 
printf ("\n\n") ; 

printf ("What  is  the  number  of  your  last  test  file?  ") ; 

scanf  ("*/,u"  ,&end_test)  ; 

printf ("\n\n") ; 

begin_degree=begin_test ; 

end_degree=end_test ; 

} 

else 

{ 

printf ("Your  first  fourier  transform  input  file  corresponds  to"); 
printf ("  what  angle  of  incidence?  ") ; 
scanf  ("’/,u"  ,&begin_degree)  ; 
printf ("\n\n") ; 

printf ("Your  last  fourier  transform  input  file  corresponds  to 
what") ; 

printf ("  angle  of  incidence?  "); 
scanf  ("*/,u"  ,&end_degree)  ; 
printf ("\n\n") ; 

} 

printf ("How  long  is  your  fiber?  ") ; 
scanf  ("'/,u"  ,&fiber_length)  ; 
printf ("\n\n") ; 

printf ("\nEnter  directory  path  of  input  file  or:  \n\n"); 
printf ("Enter  'a'  for  c : \\spir\\ftdata\\  ")  ; 
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scanf  ("*/,s"  ,directory_in)  ; 
printf ("\n\n") ; 
switch  (*directory_in) 

case  ’a’ : 

strcpy (directory. in, "c ; \\spir\\ftdata\\") ; 
break; 

} 

if (test== 'y ' ) 

{ 

printf ("\nEnter  directory  path  of  output  file  or:  \n\n")  ; 

printf ("Enter  'a'  for  c ; \\spir\\test\\  "); 

scanf  ('"/.s"  ,directory_out)  ; 

printf ("\n\n") ; 

switch  (♦directory.out) 

{ 

case  ' a' ; 

St rcpy  (direct ory.out ,  "c  :  WspirWtestW")  ; 
break; 

} 

sprintf  (f  ilename_data_out ,  "*/,stest .  dat"  ,directory_out)  ; 
sprintf  (f  ilename.doc.out  ,‘'/(stest  .doc"  ,directory_out) ; 

} 

else 

i 

printf ("\nEnter  directory  path  of  output  file  or:  \n\n") ; 

printf ("Enter  'a'  for  c : \\spir\\template\\  ") ; 

scanf  ("'/,s"  ,directory_out)  ; 

printf ("\n\n") ; 

switch  (*directory_out) 

{ 

case  ’ a' : 

strcpy (direct ory.out ,  "c :  WspirWtemplateW" )  : 
break ; 

} 

sprintf  (filename_data_out ,  "'/.stemplate .  dat"  .directory. out)  ; 
sprintf  (f  ilename.doc.out  ,"*/,stemplate.  doc"  ,  directory  .out)  ; 

} 

/♦OPEN  OUTPUT  FILES*/ 

printf ("Opening  output  file  %s  \n\n" .filename. data.out) ; 
if (  (template. data_ptr=fopen(filename_data_out , "wb") )==NULL) 


{ 

printf  ("Can't  open  output  file  %s  \n"  ,  f  ilenaine_data_out 
exit  ()  ; 

} 

printf  ("Opening  output  file  */,s  \n\n"  ,f  ilenajtie_doc_out)  ; 
if  (  (teinplate_doc_ptr=f  open(f  ilenaine_doc_out ,  "wb")  )==NULL) 

{ 

printf  ("Can’t  open  output  file  %s  \n"  ,f  ilenaine_doc_out) 
exitO  ; 

} 

for (count =begin_degree ; count <=end_degree ; count ++) 
a=0 ; 

if (test== 'y ’ ) 

{ 

f  printf  (template_doc_ptr,  "Test  vector  */.u  :  \t",  count); 

sprintf  (filename,  in ,  "5(sf  t'/(utstj(u .  256" ,  directory,  in  .count , 
fiber. length) ; 

> 

else 

{ 

fprintf (template. doc.ptr, "Template  vector  '/.u:  \t", count); 
sprint  f  (filename,  in,  "y,sft'/,udeg’/,u .  256"  .directory,  in  .count , 
fiber. length) ; 

} 

/♦OPEN  INPUT  FILE*/ 

printf ("\n\nOpening  input  file  Xs  \n\n" , f ilename.in) ; 
if (  (ptr.in=f open (filename. in, "rb") )==NULL) 

{ 

printf  ("Can’t  open  input  file  */,s  \n"  .filename. in)  ; 
exit  0  ; 

} 

I  seek (pti. in, 32000 ,0) ; 
f read(array.in, 1 , length. in ,ptr. in) ; 
f or (row=0 ; row< inrow ; row++) 

{ 

f or(col=0 ;col<incol ; col++) 

{ 

if(col>124  &&  col<132) 

{ 

vector [a] =array.in [row] [col] ; 


printf  ("'/,3u  ",  vector  [a]  )  ; 

f printf  (template_doc_ptr ,  "’/.u,"  .vectorLa] )  ; 
if (a==12| la==25| la==38) 

fprintf (template_doc_ptr ,"\n\t\t\t") ; 

a++ ; 

} 

} 

printf ("\n") ; 

} 

fprintf ( tempi at e_doc_ pt r , "\n\n") ; 
printf  ("\nWriting  vector  #y.u  \n\n‘’ .couiii.)  , 
fwrite(vector,l,a, template_data_ptr) ; 
fclose(ptr_in) ; 

} 

f close (template  ,data_ptr) ; 
f close (template_doc_ptr) ; 

> 


Aj.'jienclix  E. 


Profjraiii  DISTAS'CEX 


I  in's  proifiam  calcnlati's  distaiu'cs  Ix'tvvccti  Icinplal*'  xccioi''  in  lln'  lilc  I  1,.\! 
I’l.A  1  la.DA  I  t ()  \  ('rif \'  adixpiat <'  scparat  ion.  it  <  an  alsc.)  calculat  c  i  lie  ani  i".  lie 

twia'ii  lli<'  test  \'('ctois  in  till'  I  I'.S  r.D.AT  and  tin-  tcinpla'ic  \iTtoi''.  In  tlii^  (.i-c 
till'  ininimimi  <listanc<'  rlassilit's  ('ach  lest  vc'ctor.  I'Ik'sc  dislanci's  aic  plaii-d  iii  tin' 

li!r  ri:.\Il’DISl'.l)()('  or  ri'.STDlS  !.!)()(’  in  tin-  1  I'Al  IM.  A  H',  .a  1  IS  I  dn,  , 
r«  'I  x'ct  i\Tly.  I  !',S  1  I )  !S  I  d )( )( ’  alscj  cont  ains  the  c-'t  iinat  cd  a  n  rdn  ol  a  1 1 1  \  a  i  ( i  a  ■,  i  .  i ' 
-'dn  atioii)  for  <'a('li  test  \('(  toi'. 

I  ho  proiarani  oidy  nocils  to  know  which  distances  ai’e  re(|uired  li.e.  di'iaiici  s 
hci'aiTii  tem[)late  \'ector''  or  distances  hetwceii  test  and  I'Miiplate  ■.■.•ci  i  ,i  >  i  ami  hew- 
iiiriiiy  \'ectors  (i.c.  nuinlx'i'  ol  I'ows)  are  containeil  m  each  Idle.  1  he  imei'  can  also 
ci loose  t  !ie  input  directory,  which  should  he  eit  her  l  l  ’.M  !M,.\ T  1',  or  I  l.S  1  and  1  I  ..M 
IM  A  ( C. 


•ts 


/♦Program  DISTANCE. C*/ 


/♦Written  by  Capt  John  W.  Welker*/ 


/♦This  program  calculates  distances  between  feature  vectors.*/ 

/♦You  have  the  option  of  calculating  distances  between  the  template*/ 
/♦vectors  in  the  file  TEMPLATE.DAT  to  determine  if  your  template*/ 
/♦vectors  are  adequately  separated  in  your  feature  space,  or  you*/ 
/♦can  calculate  the  distances  between  test  vectors  in  the*/ 

Z+file  TEST.DAT  and  the  template  vectors  in  the  file  TEMPLATE.DAT.*/ 

/♦COMPILER  DIRECTIVES*/ 

#include  <stdio.h> 

#include  <string.h> 

#include  <math.h> 

/♦DECLARE  POINTERS  TO  FILE  STRUCTURE"/ 

FILE  *test_ptr , *template_ptr , ♦distance_ptr ; 

/♦START  MAIN  PROGRAM*/ 
void  mainO 
{ 

/♦DECLARE  LOCAL  VARIABLES*/ 

unsigned  char  test _f ilename [256] ,template_filename [256] ; 
unsigned  char  distance.! ilename [256] ; 

unsigned  char  test.directory [256] ,template_directory [256] ; 
unsigned  char  test_vector[49] ,template_vector[49] ; 
unsigned  countl,count2,a,min_vector=0; 
unsigned  num_test_rows  =  '^‘ .  num.template.rows  ; 
double  sum=0 . 0 , distance=0 . 0 ,min_distance=0 . 0 ; 
char  test ; 

/♦ENTER  FILES  AND  VARIABLE  VALUES  FROM  KEYBOARD*/ 

printf ("\n\nAre  you  comparing  test  vectors  to  the  template 

vectors,  y  or  n"^  ")  ; 

scanfC"  y.c"  , &test )  ; 

printf ("\n\n")  ; 


printf("How  many  vectors  are  in  yo'’.r  template  file?  "); 
scanf  C'/.u”  ,&num_template_rows)  ; 
printf ("\n\n") ; 
if (test== 'y ' ) 

{ 

printf  (."How  many  vectors  are  in  your  test  file?  "); 
scanf  ("'/,u"  ,&num_test_rows)  ; 
printf ("\n\n") ; 

} 

else 

num_te3t_rows=num_templo.te_rows ; 

} 

printf ("\nEnter  directory  path  of  template  file  or;  \n\n"); 

printf ("Enter  ’a’  for  c ; \\spir\\template\\  ") ; 

scanf  ("*/,s"  ,template_directory)  ; 

printf ("\n\n") ; 

switch  (*tempiate_directory) 

{ 

case  ' a' ; 

strcpy (template.dircctory ,"c :\\spir\\template\\") ; 
break; 

} 

sprintf  (template_f  ilename ,  "'/.stemplate .  dat"  ,  template_di rectory)  ; 
/♦OPEN  TEMPLATE  FILE*/ 

printf  ("Opening  template  file  /.s  \n\n"  ,template_f  ilename)  ; 
if  (  (template_ptr=f  open(temp''  ate_f  ilename ,  "rb")  )==NULL) 

{ 

printf  ("Can't  open  output  file  '/,s  \n"  ,template_f  ilename)  ; 
exit  0  ; 

} 

if (test== 'y ' ) 

{ 

printf ("\nEnter  directory  path  of  test  file  or;  \n\n")  ; 

printf ("Enter  'a'  for  c : \\spir\\test\\  ") ; 

scanf  ("*/,s"  ,test_directory)  ; 

printf ("\n\n") ; 

switch  (*test_directory) 


{ 

case  ' a’ : 

strcpy (tPst_directory , "c :\\spir\\test\\") ; 
break; 

> 

sprint  f  (test  .filename ,  "V.stest  .dat"  ,  test,  directory)  ; 
sprintf  (distance.f  ilename,"*/,stestdist  .doc"  , test.directory)  ; 

/♦OPEN  TEST  FILE+/ 

printf  ("Opening  test  file  */,s  \n\n"  .test. filename)  ; 
if  (  (test_ptr=f  open  (test. filenajne,  "rb")  )==NULL) 

{ 

printf  ("Can't  open  input  file  '/,s  \n"  .test.f  ilename)  ; 
exitO  : 

} 

> 

else 

{ 

sprintf  (test.f  ilename,  "'/.stemplate .  dat"  ,  tempi  at  e.di  rectory) ; 
sprintf (distance.f ilename,  '"/stempdist . doc" .template.directory) 

/♦OPEN  TEST  FILE^/ 

printf  ("Opening  test  file  */,s  \n\n"  .test.f  ilename)  ; 
if  (  (test.ptr=f  open  (test.f  ilenajme,  "rb") )  ==NULL) 

printf  ("Can't  open  input  file  y.s  \n"  ,  test.f  ilename)  ; 
exit  0  ; 

} 

} 

/♦OPEN  DISTANCE  FILE^/ 

printf  ("Opening  distance  f:ile  */,s  \n\n"  ,  distance.f  ilename)  ; 
if (  (distance_ptr=f open (distance.f ilename, "wb") )==NULL) 

{ 

printf  ("Can't  open  input  file  */,s  \n"  .distance.f  ilename)  ; 
nxitO  ; 

} 

for (count  1=0 ; count l<num_test. rows ; count 1++) 

{ 
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f read (test .vector , 1 ,49 , test.ptr) ; 
f seek (template.ptr ,0,0); 
min_distance=5. OEIOO ; 

f or (count2=0 ; count2<num_template_rows ; count2++) 

{ 

f read(template_vector , 1 ,49 ,template_ptr) ; 
sura=0 . 0 ; 
distance=0 . 0 ; 
f or (a=0 ; a<49 ; a++) 

{ 

sum+= 

pow ( ( (double) test .vector [a] - (double) template. vector [a]  ) , 2 . 0) ; 

} 

distance=sqrt (sum) ; 
if (distance ! =0 .0) 

min. distance  =  (min. distance  <  distance)  ?  min. distance  ;  distance 
if (min.distance==distance) 
min_vector=count2 : 

} 

if (test== 'y ' ) 

{ 

printf("Test  vector  '/.u  arrived  at  */,u  degrees  \n\n'', 
countl ,min_vector) ; 

fprintf  (distance.ptr , "Test  vector  V.u  arrived  at  */,u 
degrees  ",  count  1 , rain. vector) ; 

fprintf (distance.ptr ," (Distance  =  Xf)  \n\n" , min. distance) ; 

> 

else 

{ 

printf("Pass  '/,u  \n"  ,countl+l)  ; 

fprintf  (distance.ptr ,  "The  minimum  distance  of  ’/,f  occurs  ", 
min. distance) ; 

fprintf  (distance.ptr ,  "between  vector  '/,u  and  vector  '/,u  \n\n"  , 
countl , min. vector) ; 

} 

> 

f close (test. ptr) ; 
f close(template.ptr) ; 
f close (distance.ptr) ; 

} 


!()2 


Appendix  F.  Program  2562SPIR .C 


Tills  program  is  used  to  expand  a  256  x  256  data  file  to  a  512  x  512  SpiiT  oti 
fd('  for  viewing  by  the  Sjiirieon  system.  Tlu*  nsi'r  onl\'  needs  to  enter  a  dirertery 
patli  and  filename  without  extension.  The  program  automatically  appends  ,256  to 
the'  end  of  the  input  filename,  'fhe  outjuit  file  from  this  ])rogram  is  a|)pen(led  with 
a  ,512  to  indicate  that  it  is  a  512  x  512  file.  It  expands  the  file  by  zero  filling,  inst'its 
the  two  byte  header,  and  then  appends  the  SPIR.FIG  file  to  the  end  so  that  Spiricon 
will  recognize  it.  The  program  expands  the  file  within  a  few  sc-conds. 


/♦Program  2562SPIR.C*/ 


/♦Written  by  Capt  John  W.  Welker^/ 


/♦This  program  converts  a  256  x  256  file  previously  reduced  from  a+/ 
/♦512  X  512  spiricon  file  back  into  a  512  x  512  spiricon  file.^/ 


#include  <stdio.h> 

#include  <string.h> 

#include  <math.h> 

#define  inrow  64 
#define  incol  256 

#define  slice  8 

#define  outrow  64 
#define  outcol  512 

FILE  ♦ptrl , ♦ptr2 , ♦ptrS ; 

void  mainO 

{ 

unsigned  char  filename [256] ,filename_in [256] ,f ilename_out  [256] ; 

unsigned  char  array_in [inrow] [incol] ,array_out [outrow] [outcol]; 

unsigned  long  length_in=0 ,length_out=0 ; 

unsigned  row , col , count  1 ; 

char  header []  =  {’\x39’  ,  '\x30’}; 

length_in=inrow^incol ; 
length_out=outrow^outcol ; 

printf ("\n\nEnter  directory  path  and  filename  of  input  file  without 
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extension;  \n\n") ; 
scant  ("'/,s"  .filename)  ; 

sprint  f  (filename,  in ,  "*/,s  .  256"  .filename)  ; 
sprintf  (f  ilename_out .  "‘/,s  .  512"  .filename)  ; 
printf  ("Opening  input  file  '/,s  \n"  ,f  ilename_in)  ; 
ptrl=f open (f ilename_in . "rb") ; 

printf  ("Opening  output  file  ‘/,s  \n\n"  ,f  ilename_out)  ; 

ptr2=f open(f ilename_out . "wb") ; 

fwrite(header. 1 .2 .ptr2) ; 

for (count  1=0 ; count 1< si i ce ; count 1++) 

{ 

printf  ("reading  slice  #'/.u  \n\n"  ,countl+l) ; 
if(countl>l  &&  countl<6) 
f read(array_in. 1 .length_in.ptrl) ; 
f or (row=0 ; row<outrow ; row++) 

{ 

f  or(col=0  .•col<outcol  ;col++) 

{ 

if(countl>l  &&  countl<6) 

{ 

if(col>127  &&  col<384) 

array.out [row] [col] =array_inCrow] [col-128] ; 
else 

array.out [row] [col] =0 ; 

} 

else 

array_out [row] [col]=0; 

} 

} 

printf  ("Writing  slice  #y,u  \n\n"  .  countl  +  1) ; 
fwrite(array_out , 1 ,length_out ,ptr2) ; 

> 

printf ("Opening  input  file  c : \\spir\\spir . f ig  \n" .filename) ; 
ptr3=fopen("c:\\spir\\spir .f ig" ."rb") ; 
f read (array _out . 1 .906 .ptr3) ; 
f write (array_out . 1 . 906 .ptr2) ; 
f close(ptrl) ; 
f close(ptr2) ; 
fclose(ptr3) ; 

} 


Dibliography 


1.  Cole,  Capt  Geno.  Angk  of  An'ivnl  Ijelecfor.  Patent  Apf)lirat  ion:  APIS  lot).  Air 
Force  Weapons  Laboratory.  Kirtland  AFB  NM,  1988. 

2.  Gaskill,  Jack  D.  Linear  Fourier  'l’rane:fonus,  anti  OpUc.-t.  .New  ^’ork: 

John  Wiley  iC  Sons.  1978. 

9.  Ginsburg,  Arthur  P.  \  isual  Information  Processing  Bast li  on  Spaiial  L'/lh  rs 
Constrained  hij  Biological  Data.  Doctoral  Diss('rtat ion.  A.MRL  1  H  -  78  129, 

Wright-  Patterson  ,\FB  Oil,  December  1978. 

1.  Hecht,  Eugene.  Optics.  Heading  MA:  Addison-Wesley  Publishing  (^oinpany. 
1987. 

•').  Holl,  Herbert  B.  .Angle  of  Arriral  Meter.  Patent  Application;  S.N.  919128.  I’n- 
knovvn.  1982. 

f).  Kabrisky,  Matthew  E.  Class  lecture  in  EE.NG  621,  i’attern  IC'cognition  II. 
School  of  Engineering,  Air  Force  Institute  of  Technology  (AF).  Wright - 
Patterson  AFB  OH,  April  1989. 

7.  Keiser.  Gerd.  Optical  Fiber  Communications.  .New  ^’ork:  McGraw  Hill  Book 
Company,  1983. 

8.  Letourneau,  Capt  George  B.  .Aspect  .Angle  Inrariant  .Sgnllitlic  .Aptr- 
ture  Radar  Pattern  Recognition  Using  the  Hough  L'ran.^form.  .MS  thesis. 
AEIT/GE/ENG/89D.  School  of  Engineering,  Air  Force  Institute  of  Jechnol- 
ogy  (AU),  Wright-Patterson  AFB  OH.  December  P989.  (Dl'lC  uumlx-r  not 
available  at  this  time). 

9.  Optical  Sensor  With  High  Directional  Resolution.  SBIR  Phase  1  Report.  Topic 
AF88-184.  Interactive  Intelligent  Imagery  (’orporat  icni,  Nh-nlo  Park  C.\.  .March 
1989. 

10.  Parker.  Jack  H.  Jr.  Common  Opto -Fleet ronic  Imser  Diliclion  Sgsttm 
(COLDS).  Final  Report,  July  1985  -  August  1985.  AFWAL  J'R  86  -  K)18, 
Wright  Patterson  AFB  OH,  December  1986  (AD  -  B1  19557). 


South  Carolina  in  November  1978  and  entered  the  University  of  South  Carolina 
in  January  1980.  In  December  1983,  he  enrolled  in  the  College  Senior  Engineering 
Program(CSEP)  of  the  United  States  Air  Force.  After  receiving  a  Bachelor  of  Science 
degree  in  Electrical  Engineering  in  December  1984,  he  attended  Officer  Training 
School  and  was  commissioned  a  Second  Lieutenant,  5  April  1985.  Following  a  three 


year  assignment  at  Edwards  Air  Force  Base,  where  he  performed  Defensive  Avionics 
testing  for  the  Bl-B,  he  entered  the  Electro-optics  program  at  the  Air  Force  Institute 
Technology(AFIT). 


UNCLASSIFIED 


SECURITY  cla5s;p‘Ca:iotj  of  this  ^age 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


h'orm  Approved 
CMB  No.  0704-0138 


la.  REPORT  SECURITY  CLASSiFiCATlON 

UNCLASSIFIED 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

2b.  DECLASSlFlCATION/OOWNGRAOiNG  SCHEDULE 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER{S) 

afit/ge/ei:g/89d  -  57 

6a.  NAME  OF  PEHFO.RMiNG  ORGANIZATION 

6b  OFFICE  SYMBOL 

(If  applicable) 

School  of  Engineering 

AFIT/ENP 

6c.  ADDRESS  {City,  State,  and  ZIP  Code) 

Ai.r  Force  Institute  of  Technology  (AFIT) 

(AU)  VJPPFB,  CH  45433  -  6583 

'  3a.  N.'.ME  Hie  e.  :^Jn|^,^3  SPONSORING 

3b  OFFICE  SYMBOL 

ORGANIZATION 

(If  applicable) 

Air  Force  Weapons  Lab  (WL) 

3c,  ADDRESS  (Oty,  State,  and  ZIP  Code) 

iTL Atget 

Kirtland  AFB,  37117  -  6008 

3  DISTRIBUTION /AVAILABILITY  O'  REPORT 

Approved  for  public  release; 
distribution  unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUM3ER;S) 


b  AonREsS{Oty,  State,  ^nd  ZIP  Code) 


’0.  SOURCE  OF  'FUNDING  NUMBERS 


PROGRAM  I  PROJECT 

ELEMENT  NO  I  NO 


WORK  UNIT 
ACCESSION  NO 


11.  TITLE  (Include  Security  CUiSification) 

ANGLE  OF  ARRIVAL  DETECTION  THROUGH  ANALYSIS  OF  OPTICAL  FIBER  INTENSITY  PATTERNS (UNCLASSIFIED) 


12.  PERSONAL  AUTHOR<S) 

John  W.  WeLker,  Captain ,  USAF 


13a,  TYPE  OF  REPORT  1 3b.  TIME  COVERED 

MS  thesis  FROM _  TO _ 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  15.  PAGE  COUNT 

1989,  Deceirber  107 


17.  COSATI  CODES 


FIELD  GROUP  SUB-GROUP 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Angle  of  Arrival,  Optical  Fibers,  Pattern  Recognition 

12  09  I  I 

20  06  Ol 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number)  Thesis  Advisor:  Maj  S.  Rocers ,  PhD 

Associate  Professor,  Dept  of  Electrical 
and  Ccnputer  Engineering 


This  thesis  examined  the  feasibility  of  using  optical  fibers  for  angle  of  arrival  (AOA) 
detection.  Specifically,  a  uniform  anplitude  plane  wave  was  transmitted  through  a  lirm 
diameter,  multimode,  step  -  index,  plastic  fiber.  Fiber  lengths  of  lin,  2in,  and  3in  were 
investigated.  The  output  intensity  pattern  from  each  optical  fiber  was  detected  by  a  512  X 
512  resolution,  charge  injection  devic:e  (CID)  camera.  The  data  array  fron  the  camera  was 
processed  with  a  median  filter,  and  normalized  to  the  peak  intensity.  A  fast  Fourier 
transform  converted  the  data  to  the  frequency  dcmain  where  the  fundamental  and  first  three 
harmonic  frequencies  were  extracted.  These  49  numbers  represented  components  of  a  feature 
vector  representing  the  pattern  related  to  a  specific  ADA.  Twenty  -  six  template  vectors, 
each  representing  a  1  degree  increment  from  0  degrees  to  25  degrees,  were  created.  Test 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

^UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT  □  DTIC  USERS  LTvCLASSIFIED  _ 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b.  TELEPHONE  (Indude  Area  Code)  22c  OFFICE  SYMBOL 

Steven  K.  Pogers,  Asscciate  Professor  (513)  255  -  9266  AFIT/ENG  _ 


DD  Form  1473,  JUN  86  Previous  editions  are  obsolete.  _ SECURITY  CLASSIFICATION  OF  this  page 

UNCLASSIFIED 


19.  vectors  were  coipared  to  template  vectors  using  Euclidean  distance.  The  mirdmum 
distance  classified  the  test  vector  AOA. 


